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ABSTRACT 

This report explains the fundamental theory and equations 
used in writing a set of software applications which predict 
antenna radiation parameters. Each application predicts the 
radiation characteristics of a particular type of antenna over 
a planar surface which serves as a model of either earth or 
seawater. The radiation parameter predictions are based 
solely on an antenna's physical dimensions, the properties of 
the underlying surface, and electromagnetic theory. Existing 
electric field equations provide the basis for radiation 
parameter predictions, and the accuracy of the predicted 
radiation parameters is totally dependent upon the extent to 
which the electric field equations used realistically model 
the actual radiated electromagnetic fields of the antennas. 

In addition to a review of applicable electromagnetic 
field theory, this report is also intended to be a user's 
guide for the corresponding computer applications. The 
appendices contain computer hardcopies of sample calculations 
for several antenna types and remarks regarding the conformity 
of predicted radiation parameters to expectations. Radiation 
parameters computed thus far are consistent with expectations 
based on other computational programs and empirical 


measurements. 
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I. BACKGROUND AND PURPOSE 

This report and associated computer software applications 
are submitted in fulfillment of the thesis requirements for 
the degree of Master of Science, Aeronautical Engineering from 
the Naval Postgraduate School in Monterey, CA. The thesis 
text explains the fundamental theory and equations used in 
writing a set of software applications which predict an 
antenna’s radiation parameters. Each application predicts the 
radiation characteristics of a given type of antenna over a 
planar surface which serves as a model of either earth or 
seawater. The thesis requirement was generated by a request 
from the Naval Maritime Intelligence Center (NAVMARINTCEN) . 
At the request of NAVMARINTCEN, Mathcad mathematical software, 
DOS version 3.1, was used to write the applications. The 
NAVMARINTCEN request was extensive enough to be tasked as two 
separate thesis topics. This thesis fulfills the first half 
of the NAVMARINTCEN requirements. Mathcad applications are 
complete for the following types of antennas: 
Vertical Monopole Antenna 
Elevated Vertical Dipole Antenna 
Elevated Horizontal Dipole Antenna 
Arbitrarily Oriented Dipole Antenna 
Inverted L Antenna 
Long-Wire Antenna 
Terminated Sloping V Antenna 
Side-Loaded Vertical Half Rhombic Antenna 
Terminated Sloping or Horizontal Rhombic Antenna 
10. Sloping Double Rhomboid Antenna 
11. Vertically Polarized Log-Periodic Dipole Array 


12. Horizontally Polarized Log-Periodic Dipole Array 
13. Horizontal Yagi-Uda Array 


DONAN PUNE- 


d 


II. INTRODUCTION 

New American Heritage Dictionary defines antenna as 'an 
apparatus for sending and receiving electromagnetic waves'. 
An electromagnetic wave 1S a time-varying, self-propagating, 
interrelated combination of electric and magnetic fields. 
Antennas radiate electromagnetic fields as a function of their 
time-varying surface charge density and surface current 
density distributions. There are many different types of 
antennas, each with its own set of radiation characteristics 
or radiation parameters. These radiation characteristics are 
related to an antenna's radiated electromagnetic field 
distribution, and they determine the useful applications for 
a particular antenna design. 

Current personal computers with high-speed processors can 
quickly and accurately predict antenna radiation parameters 
from electromagnetic theory. The Mathcad routines described 
in this report predict an antenna's radiation parameters based 
solely on its physical dimensions and electromagnetic theory. 
Accurately predicted radiation characteristics provide 
intelligence analysts with a reliable estimate of an antenna's 
capabilities, limitations, and vulnerabilities. Chapters III 
and IV are a review of applicable electromagnetic theory, and 
the remaining chapters describe the calculations of the 
associated Mathcad computer code, and are intended to be a 
user's guide for the associated Mathcad applications. 
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III. FUNDAMENTAL ELECTROMAGNETIC THEORY 
A. MAXWELL'S EQUATIONS AND HELMHOLTZ'S EQUATIONS 
Maxwell's equations are the basis for electromagnetic 
field calculations. Their differential form is given by [Ref 


O 321-325] 


= — OH s - Y 
NOS E ba e Ta) UE : E) 
VXH-J«e$E Siang) V-H-0 (3.1d) 


Equation 3.1a is Faraday's Law, 3.1b is Ampere's Circuital 
Law, 3.1c is Gauss's Law, and 3.1d postulates the inexistence 
of magnetic charge. 

Maxwell's equations dictate the relationship between the 
electric and magnetic field intensities (E and H) and the 
Charge/cuüurrent density source] dTstraibutions (p and J) (ài.e., 
the antenna). Maxwell's equations are often applied in a 
source-free, current-free region, and an e'* time dependence 
is assumed. The time-harmonic, free space Maxwell's equations 


are [Ref 1: p. 340] 


VXEs--jopg,H (3.2a) Dae, = 0 sens. 2C) 
VXH= jwe,E (3.2b) V:-H=0 (3.2d) 
where 

] eae 1 2 Farads 

¢,=free spac Any. = 10° | AZ 

: es ett "— 
— 1| Henrys 
unoocscruasccascnmedbsdbtysudm - 10°) |———-— 
Fo d E ^ eS 


Rather than solve Maxwell's equations directly, the 
method of vector potentials is often used where the scalar 
potential (V) and magnetic vector potential (A) can be 
obtained from the nonhomogeneous Helmholtz equations 

Vv + B2V = 2 (3.3) VA + BA = -p,J (3.4) 

0 
where Hs(i/nug)(V X A) [Ref 1: pp. 338-340]. 

The solutions to Helmholtz's equations are the retarded 
scalar potential and the retarded vector potential. The value 
of V and A at some distance (R) from the source depends on the 
source's charge/current density at an earlier time (t-R/c), 
where car is equal to the speed of light in free space. 
The delay (R/c) is the time required for electromagnetic waves 
to propagate through a distance, R, in free space from source 
to observation point. The solutions to the nonhomogeneous 


Helmholtz equations are [Ref 1: pp.338-340] 





-jB|R-&'| 
R)= 1 f PEZ (SO) 
4TE, y! | RR | 
-jp|R-E"| 
T ICAA (376) 
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where the wavenumber (p) is EAT in free space, R is the 
observation vector, and R' is the source vector. The wave- 
number is most often expressed as w/c = 2m/X, where A is the 
wavelength. The Mathcad applications use the expression 27/A 


to calculate f. 


The Mathcad applications described herein assume that the 
antenna under analysis is radiating in free space above a 
defined ground plane. The free space Maxwell's equations and 
nonhomogeneous Helmholtz equations (using the free space 
wavenumber) provide the theoretical basis for deriving the 
equations for an antenna's radiated electromagnetic fields. 
The radiated E and H fields can be calculated directly from 
their source's charge and current densities using Maxwell's 
equations. However, aS previously mentioned, it is often 
Simpler to find the magnetic vector potential using 
Helmholtz's equations. Once the vector potential is 
determined, the electric and magnetic fields are calculated 
from H-(1/ug)(V X A) and E=(1/jwe,) (V X H) [Ref 1: pp. 338- 
341]. Existing theoretical equations for the radiated 
electric fields of many antenna types have been derived in 
this manner. These equations are used in the Mathcad 
applications to predict an antenna's radiated electric field 
distribution. Other radiation parameters are then calculated 
on the basis of the predicted electric field distribution. 
B. ELECTROMAGNETIC FIELD REGIONS 

The space surrounding an antenna is divided into three 
regions: the reactive near-field, radiating near-field, and 
far-field [Ref 2: pp. 22-24]. The reactive near-field 
occupies the space immediately surrounding the antenna out to 
a radius of about 0.62- (D'/A)5, where À is the wavelength and 


DASS maximum wdxmension of the “antenna. Reactive 


electromagnetic fields dominate in this region I Tea UT 
near-field occupies the space from the boundary of the 
reactive field out to a radius of about 2-(D'/A). Radiation 
fields dominate in this region, and radiation patterns are 
often a function of both radial and angular coordinates. The 
far-field region occupies all space outside the radiating 
near-field which meets two additional far-field requirements: 
27TR/X >> 1 and R > 5D. Radiation fields also dominate in the 
far-field, but field components are primarily transverse, and 
radiation patterns are normally a function of directional 
variables only. The Mathcad antenna applications are valid 
only in the far-field region of a given antenna. 
C. RADIATION PATTERNS 

A radiation pattern is a three-dimensional representation 
of an antenna's radiated electromagnetic field distribution or 
power distribution at a fixed distance from the antenna. 
Because it is difficult to depict three-dimensional images, 
the patterns are most often plotted in a defined vertical or 
horizontal plane. The Mathcad antenna applications plot two- 
dimensional radiation patterns in polar coordinates depicting 
the far-field electric field distribution. Patterns are 
plotted in the $20 and $-27/2 vertical planes, and also in a 
horizontal plane at an elevation angle set by an index from 
the applicable 'Elevation Angle Index Table' (Table 3.1 or 
3.2). Field magnitudes are normalized with respect to the 
maximum radiated electric field intensity. The magnitude to 


which each pattern is normalized is displayed below its plot. 
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TABLE 3.2 
ELEVATION ANGLE INDEX TABLE 
FOR LOG-PERIODIC AND YAGI-UDA ARRAYS 
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BeMMSthcad “applications plot the far-field radiation 
patterns for an antenna's space wave, surface wave, and total 
radiated electric field. The space wave is comprised of the 
electromagnetic waves which propagate directly from antenna to 
the field point and electromagnetic waves reflected from the 
surface below the antenna. The surface wave is composed of 
electromagnetic waves ducted along the surface, much like 
waveguide propagation. The total field is composed of the 
Space wave and surface wave combined. 

Space wave patterns typically vary with directional 
spherical coordinates, @ and 6$, but they are predominately 
invariant with distance from the antenna. The surface wave 
attenuates exponentially with distance, and far-field surface 
wave patterns are a function of distance from the antenna. 
Total field patterns are still predominately invariant with 
distance, because the maximum surface wave intensity is 
usually many orders of magnitude less than the maximum space 
wave intensity. Even at distances which just meet far-field 
requirements, the surface wave contribution to the total field 
1s barely significant m parameter values of interest. 

D. POWER DENSITY, POYNTING VECTOR AND RADIATION RESISTANCE 

The instantaneous real power flux density of an electro- 


magnetic wave is called the Poynting vector (@), 


O(x.y.z) = EX H = WRe(E(x,y,z) e??') x 8e(H(x,y,z)e?^?* (3.7) 


== RelE X H*) « = RelE X He?jet] 


The time invariant term on the right hand side of equation 3.7 
is the time average power flux density or average Poynting 


vector. Far-field E and H fields are related by the simple 


expression |E|=(1/y))-|H|, where  mqvnu,/ey-120*" is the 
intrinsic impedance of free space. The average Poynting 


vector solely in terms of the radiated electric field 
intensity is given by [Ref 1: pp. 382-385]. 


Pu (x.y, 2)=2Rele xH)=2Re Ex E = E (3.8) 
2 2 o No 


The total power flux density is composed of both the 
Poynting vector and the imaginary part of ¥%-(E X H’) which 
represents the reactive power of the radiated electromagnetic 
fields. Far-field power flux density is predominately real 
and is approximately equal to the average Poynting vector. An 
antenna's total average radiated power is calculated by 
integrating the poynting vector over a Gaussian surface which 


surrounds the antenna in the far-field [Ref 1: pp. 382-386]. 


1 " 
PEE f, Pag ds = >$ RAE X H'): ds (3.9) 


The Mathcad applications use equations 3.8 and 3.9 for power 
calculations by integrating over the hemisphere which encloses 
the antenna with the ground plane as the lower boundary. 
Radiation resistance relates an antenna's total average 
radiated power to its peak input current. If a DC current 
equal to the rms input current were applied to a resistance 


equal to the radiation resistance, the dissipated power would 
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equal the antenna's total average radiated power [Ref 3: pp. 
47-48]. For the Mathcad antenna applications, the radiation 
resistance is defined in terms of the total average radiated 


power and input current by 











2 rad 2 rad ORELO 


A sinusoidal input current (or voltage response across the 
input terminals for Yagi-Uda Arrays) with a maximum of unity 
is assumed in the Mathcad applications. 
E. DIRECTIONAL GAIN, DIRECTIVITY AND POWER GAIN 

An isotropic radiator is a hypothetical antenna which 
radiates electromagnetic energy equally in all directions. 
Its radiation pattern in any plane passing through the source 
is a circle. All real antennas are directional radiators 
which radiate or receive electromagnetic energy more 
effectively in some directions compared to others. Some 
antennas radiate an omnidirectional pattern which is 
directional in elevation, but non-directional in azimuth. 

uc ATA DA IS the trato OF an antenna’s Poynting 
vector ina given direction, divided by the power flux density 
of an isotropic source with equal total average radiated 
power. Directivity (D is the maximum value of the directive 
gain, and it is the primary measure of an antenna's 
directional properties. The power flux density of an 
isotropic radiator is equal to the total average radiated 


power divided by the surface area of a sphere at a given 


p 


distance from the antenna (47R’). The expression for 


directive gain is [Ref 2: pp. 29-31] 


_ ©(8,) _ 4nR*0 (8,6) 
Dg (8,6) == = 


liso rad 


D 


The Mathcad applications calculate only directivity using 


9 e max) E anR*® On Os.) 
c P 


liso rad 


D,= (ONES 


Power gain accounts for an antenna's radiation efficiency 
in addition to its directional properties. Power gain (G,) is 
defined as 4TR* times the ratio of the Poynting vector 
magnitude in a given direction to the total input power to the 
antenna. Total input power includes losses due to ohmic power 
dissipation and impedance mismatches between transmission line 
and antenna. The ratio of total average radiated power to the 
total input power is the radiation efficiency, ).,- The 
relationship between power gain, directive gain, and radiation 


efficiency is given by [Ref 2: pp. 42-45] 





P 
P (3.13) 
G4(0,$) -p,, D, (8, $) (3.14) 
G¿=Gg, (8, $) =PerrD, (9,9) =PoreDo (3.15) 


Since radiation efficiency and power gain account for 
ohmic losses and impedance mismatches, these parameters are 
difficult to predict based on antenna dimensions alone. For 


these reasons, the Mathcad applications do not predict power 
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gain and radiation efficiency. Gain calculations are limited 
bO directivity since it depends only on total average radiated 
power and the radiated electromagnetic field distribution. 
The maximum Effective Isotropic Radiated Power (EIRP) is also 
Saleculated, from EIRE- D [Ref 37 p. 62]. 
F. HALF-POWER BEAMWIDTH 

Radiation lobes are portions of the radiation pattern 
where a local maximum field intensity is bounded on either 
side by relative minimum intensities. The major lobe, or 
main beam, is the lobe in the direction of maximum radiated 
power flux density. Split beam antennas have more than one 
main beam. A minor lobe is any lobe except the main beam. 
Beamwidth refers to the half-power beamwidth of an antenna's 
main beam (or beams), and it can be measured in the vertical 
Giuminori zontal . For antennas with a defined beamwidth, the 
power intensity drops off on either side of the point of 
maximum radiated average power density until it is one half of 
the maximum value. The angle between the radials which 
intersect the half power points is defined as the half-power 
beamwidth. Main beams are not necessarily symmetrical since 
the radial which intersects the point of maximum power density 
does not necessarily bisect the half-power radials. 

The Mathcad applications do not calculate half-power 
beamwidth directly, but radiation patterns provide a good 
estimate. Since the radiation patterns depict the electric 


ElevemmcastributiLon smormalezed to “Unity, a magnitude of 


ds 


approximately 0.7 represents the half power-points of the 
pattern. If two radials are drawn from the origin through the 
half-power points, the angle subtended by the two radials is 
the antenna's half-power beamwidth. 
G. POLARIZATION AND BANDWIDTH 

Polarization describes the time varying properties of an 
antenna's radiated electric field at a set spatial coordinate 
as electromagnetic waves propagate outward from the antenna. 
Signals are usually transmitted with linear, elliptical, or 
circular polarization. Linear polarization is where the field 
vector at the spatial coordinate is directed along a fixed 
linear path. An elliptically polarized signal is one in which 
the tip of the electric field vector traces an ellipse around 
the fixed coordinate as the wave propagates outward. Circular 
polarization is a special case of elliptical polarization 
where the magnitude of the electric field remains constant as 
it rotates about the spatial coordinate. Clockwise rotation 
of the electric field vector is right-hand polarization, and 
counter-clockwise rotation is left-hand polarization. The 
Mathcad applications do not quantitatively predict the 
polarization at a selected spatial coordinate because of the 
computational intensity of such a calculation. Only the 
predominant type of polarization expected is addressed briefly 
in the introductory remarks of each application. 

The bandwidth of an antenna is the range of frequencies 


over which its radiation characteristics meet or exceed a 
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specified acceptable performance standard. Because radiation 
characteristics of dissimilar antennas vary differently with 
frequency, there are no set rules or guidelines which define 
bandwidth for all antennas in general. Bandwidth standards 
vary depending on antenna design and intended application. 
üEcstore, the Mathcad applications do Nok predict operational 
bandwidth. Instead, multiple frequencies can be input for a 
given configuration to determine the bandwidth for whatever 
standards are used to define acceptable performance. 
H. EFFECTIVE LENGTH AND EFFECTIVE AREA 

The effective area, or effective aperture, gives the 
power delivered to the antenna load when multiplied by the 
incident average power flux density. It is defined by [Ref 2: 


pp 59-63] 


E (3.16) 


inc 
The maximum effective area is related to directivity for all 
antennas by [Ref 4: pp. 46-47] 


_ 42D, 


EET 
TOEF ( ) 





This is a theoretical maximum effective area since directivity 
does not account for radiation efficiency. Like directive 
gain and power gain, actual effective area and theoretical 


effective area are related by radiation efficiency 


À*G, 


Bert” Ta Pett nax (3.18) 





LS 


The maximum theoretical effective area is a good approximation 
of actual effective area for frequencies within the bandwidth 
of a well designed antenna with a high radiation efficiency. 

Effective length is defined as the ratio of the voltage 
induced across the antenna terminals by the incident electric 


field and is given by [Ref 4: pp. 40-42] 


V. 
1 E das 
aes [EN | EA 


Effective length is related to effective area by 


R 
TANZ Ses ras (SON) 
0 


Substituting equation 3.17 into equation 3.20, we get the 
maximum theoretical effective length in terms of wavelength, 
directivity, radiation resistance and intrinsic impedance of 


free space as 


K D 
4TNo 


l __=2 (321) 


- 
Finally, the actual effective length is related to the 
theoretical effective length by 

occ Door max (3.22) 


Effective area and effective length are actually a 
function of directional coordinates with respect to the 
source. The equations given here are for maximum values of 


effective area and effective length along the radial of 
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maximum directive gain (directivity). While effective area is 
a parameter more generally applied to all antennas, the 
Mathcad applications use equations 3.17 and 3.21 to calculate 
both a theoretical maximum effective area and effective 
length: Penance S adaton efficiency for ay given 
frequency is known, the actual effective area and length can 
be calculated by multiplying the theoretical values by the 
racdgwson etficieney. 

The next chapter on the dipole antenna explains in detail 
the derivation of the radiated electric field equations for a 
vertical dipole antenna. The elevated vertical dipole Mathcad 
application uses this equation to obtain the electric field 
distribution and total average radiated power (P.,,) based on 
an assumed Sinusoidal current distribution on the antenna with 
a maximum of unity. Total average radiated power is then used 
to predict the radiation parameters discussed in the preceding 
paragraphs. In follow-on chapters, only the final electric 
field equations are presented because the derivation is 
Similar for all types of antennas. They are primarily summary 
chapters which describe the calculations of each individual 
application, but any new or important concepts not presented 


in previous chapters will be discussed as required. 
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IV. THE VERTICAL DIPOLE ANTENNA 

A. THE ELEMENTAL DIPOLE 

The current element, “or” Clemencal a1 vel, ls err 
building block for calculating antenna radiation parameters. 
The elemental dipole is a hypothetical differential length 
thin-wire antenna with a sinusoidally oscillating surface 
current. The antenna length (1) is many orders of magnitude 
smaller than the wavelength (A), so the instantaneous current 
is considered uniform over the length of the dipole. Figure 
4.1 depicts the elemental dipole oriented along the z axis in 


a rectangular coordinate system with its center at the origin. 





FIGURE 4.1: Elemental dipole spatial orientation. 
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The surface current is given by 


1(t)=I cos (wt) =Re(Ie?“*) (4.1) 


and the magnetic vector potential, A, is found from equation 


3.6 to be 











-jB& Tal -JBR 
ai Je- dl'=a, +? - | (4.2) 
4AtJc’ R 4 T R 


Primed coordinates refer to source points, R is the distance 
Ao IMSS gin "ro the fierd point, B=Wu é> is the free space 
wavenumber, and the magnitude of the current density (J) is 
Givens, 1 [Ret 2: “pp.” 100-102]: The magnetic vector 


potential in spherical coordinates is 








] i Ew mun | 
A,-A, cos (0) e c CN cos (0) (4.3a) 
M u P | 
Ag7-A,sin(0)- n = sin (0) (4.3b) 
A¿=0 (4.30) 


The magnetic field intensity, H, is found from 





1 1109 za 
== La e S oe 4.4 
es a 22 (aa) E ( ) 
Word i -jBR 
H, =-—— S n (0) - —————— je" "s 
ah JPR (FBR)? 








and the electric field intensity, E, is found from Ampere's 














Law by 

E- Li Y x ma wL SO IS 

jwe, jwe,| *Rsin(0) 00 ^? Ron i 
Tal 2 1 i P 
E zc c lets A (4383) 
Rag Nob e (GBR): 
ran oe dit l E 
Ea = - — 2sin(6) | —— = +———_ + —_——__ Je JR (4.7 b) 
oa nop (JBR) (GBR)? GBA? 
Hp=Hg= Ey, =0 (4.70) 


At far-field distances, the only significant components are 
the first terms inside the parenthesis on the right hand side 
of equations 4.5 and 4.7b. The far-field E and H fields are 
functions of a single angular coordinate (0) for a vertically 
orientated element [Ref 1: pp. 602-607]. 

Retaining only those terms that vary as 1/R, we get the 


far-field electric and magnetic field intensities for the 


current element as 








-jBR 

s, NM [EI (4.8) 
-jpR 

Ec JnoBsin (0) (4.9) 


Note that the far-field relationship |H|-*,|]E| holds true for 
the elemental dipole's radiated fields in the far-field 
region. 


The elemental dipole's radiated average Poynting vector 


(power density) is [Ref 1: 607-612] 
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_i re 


avg 


| | ss n.B?sin? (0) (4.10) 


The Poynting vector integrated over a Gaussian surface in the 
far-field yields the total average radiated power 
(Iddl) i 
Prad=$ Pavg AS = -n Pf Ji sin?(8)]sin(0)d0d$ (4.11) 


(Id1)? E | (Idi)? 4 
P -—ÁÀ E ERE ua E RE 0 I ms ces 
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D 
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The directive gain and directivity are found from the 


power density and total radiated power as 


AnR*O amsin (0) 3... 
E avg _ 2 2 
D,(0,9) Ea SEO sin (0) (4.12) 
3 
Dy=D,,,, (8,6) =D 3 0 0)=1.5 or 10 log, 1.5-1.76 dB CURES 


Radiation resistance is found directly from P,, and equation 


2.10 as [Ref 1: pp. 607-612]. 
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B. THE ELECTRICALLY SHORT DIPOLE 

The derivation of the electric field for an electrically 
short dipole is the next step in deriving the field equations 
for a practical finite length dipole. The electrically short 
dipole is a center fed thin-wire antenna oriented as in Figure 
4.1 with length greater than 0.02A, but less than 0.10A. The 


current distribution is estimated as triangular where 


aI, (2-52, 0sz/«d 
De a ze eee 2 E (4.16) 
a.I, zal “5 $2's0 


where I, 1s the current maximum ae e2eie@ewe sea. 


The dotted line in Figure 
4.2 is a more exact 
representation Of the 
actual current distri- 
bution, but the triangular 
estimate 1S a reasonable 


approximation for short 





FIGURE 4.2: Current distribution dipoles. 


for an electrically short dipole. 
The electrically 


short dipole has a magnetic vector potential given by 
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Because dipole length is small in comparison to the far-field 
distance, R is considered constant throughout the integration 
over the length of the dipole. Note that the magnitude of the 
magnetic vector potential for the electrically short dipole is 
half that of the elemental dipole for the same current input. 
The equations for the radiated electromagnetic fields are 
derived in the same way as for the elemental dipole, and the 
final expressions for the electric and magnetic field 


distributions of the electrically short dipole are given by 


eo. lowe, 
DAS MES S . 
e sin(0) (4.18) 
¿PI _le3PR 


O 


Direc a nftend directivity forthe electrically short 
dipole are the same as those of the elemental: dipole.The 
magnitudes of the radiated power density, total radiated 
power, and radiation resistance for the short dipole are given 


Diane mea: PP- 109-112] 


2 
ep Lo; eB sin! (0) (1.209 
2 
ra n p? (4-21) 
R,.4-20m? (=) (4.22) 
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C. THE FINITE LENGTH DIPOLE 

Finite length dipoles are center fed thin-wire antennas 
with negligible diameter (diameter << à and radius << length). 
Dipoles with a length greater than 0.10A can no longer use the 
assumption of triangular current distribution with acceptable 
accuracy The finite length dipole must account for the 
sinusoidal current variation along the length of the antenna. 
The sinusoidal surface current distribution for a finite 


length dipole oriented as depicted in Figure 4.1 is 


a.I, sin[B(h-z^], 0sz/«d 
I(x 2) IE 7 (4.23) 
a,I, sin(B(h«z/], “5 $2's0 


were I 19S the poak Orren OUI 


Nlu 


Considering the finite length dipole as an aggregate of 
elemental dipoles, we get the expression for the electric far- 
field from équation 4.7b as (Ref 2: pp. 118120] 


Nob 
atti 





&-[ Ad&-) sin(0) f T(x’, y’, z’) ezaz A 

2 

For purposes of the radiated electric field's phase 
angle, the radial field distance (R) cannot be considered 
constant while integrating over the length of the dipole. The 
variable r (which stands for |R-R' |) in the exponential 
(phase) term is approximated as [R-z'cos(0)], where R is the 
distance from the origin to the field pomt fama 2 Ene 
source coordinate of the vertical dipole along the z axis 
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E ADA. The 
approximation is depicted 
Ti —Pagures49398 and" is 


acsumatesnsn the far field 





—- z'cos (0) Ron ccc o Nu For 


Far Field Approximation purposes of calculating 


r=R-z'cos (8) , | | 
electric field magnitude, 





FIGURE 4.3: Phase approximation the variation of r over 


for distance variable, r. 
the Tength of thei dipole 


is negligible, and the distance from the origin to the field 
point (spherical coordinate R) is used in the denominator of 
equation 4.24. Post biting eguationa 23 for- I(x! yir Z) 
into equation 4.24 and [R-z'cos(8)] for r in the exponential 


term of equation 4.24 and integrating, we get 


Y nope 7^ 
Ll RR 


sin(6)/f Tsinlp (h+z‘) Je tIbz’cos (8) qz! (eS) 
-h 
+ ft sinlp (h-z/)] e*JBz cost?) qz! 
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- nN pe 75 


sin(0): 
¿mE 


/,sosinlp (h«z^) [cos [Bz’cos (@) ]+jsin[Bz‘cos (8) ] |dz’ 
a 
+f "rsin[f (h-z’) [cos [Bz'cos (8) ] +jsin[fz'cos (0) | Jaz" 
0 
Since the current distribution is an even function, the 
odd term of j:sin(f-z'-cos(0)) integrates to zero over the 


interval from -h to h. Since the cos(8-z'-cos(@)) term is 


even, the integral can be written as 


25 


= nope o 


EN ZEE sin(0) f "sin[B (h-2^]eos [Bz/cos (6) ) dz' (4.26) 


After integrating, the final expressions for radiated electric 
and magnetic field distributions of a finite vertical dipole 
are given by [Ref 2: pp. 118-120] 


-jBR, uan l 
mja ooe [-sesiBacos “6 | - cos (Bh) | (4.27) 


2TR Sin O" 


E TI 
m = ) 
No 2TR 


AA (4.28) 


sin(0) 





The coefficient outside the brackets is called the element 
factor, and the expression inside the brackets is the space 
fcc one 

The electric field radiation pattern is a function of 
elevation (0) and wavelength (A). For a fixed wavelength, the 
radiation pattern for any vertical plane through the antenna 
axis is obtained by plotting the magnitude of equation 4.27 in 
polar coordinates as Ü varies from zero to 2m. The pattern is 
normalized by multiplying the electric field magnitude by a 
factor such that the maximum result is one. Since electric 
field distribution is independent of $, the pattern for -ms6@s0 
is symmetrical with that for OsÜsm. Figure 4.4 depicts the 
normalized radiation patterns of four finite dipoles with 
successive lengths of A/4, A/2, 3A/4, and A. These lengths 
simplify calculation of the electric field patterns since fh 
is m/4, 1/2, 371/4, and mw for each successive dipole. If 


dipole length is greater than one wave-length, the patterns 
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become more complex and the 
number of lobes increases. 

The radiated power 
density and total radiated 
power are found just as 
they were for the elemental 
dipole and the electrically 
short dipole. The final 


expressions for radiated 





power density and total 


FIGURE 4.4: Vertical plane 


patterns for dipoles of length radiated power are [Ref 2: 
ME Ec 


1=3A/4 2s o E ; and l=A ------ ; DEM 120-124] 
2 
1 : Nolo [ cos [g hcos (8) ] -cos (B4) 1? 
0 ve H” o= 4.2 
ag? Re(E X H")=a, in LU (4.29) 
Z 2T (TX T A 
eer i EIU co (4.30) 


_ — [cos[Bhcos (8) 1-cos (BA) 1* 48 
0 


4T sin (80) 


The Mathcad applications evaluate equation 4.30 to obtain 
total radiated power for a Sinusoidal input current with a 
Maximum of unity. An expression for radiation resistance is 
also derived by Balanis [Ref 2: p. 124]. However, the Mathcad 
applications calculate radiation resistance from total average 
radiated power and equation 2.10. 

When the dipole is elevated above the origin along the z 


axis, an additional phase shift term of exp[-jfH,cos(0)] must 
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be included in equations 4.27 and 4.28, where H, is the 


elevation of the feed at the center of the dipole; hence, 





PE rc ee cos[pAcos (0)] - cos(BA) (4.31) 

: 2mR sin(0) | 

a eee I 0 2 moma cosíBAcos (0) ]-cosQ n (4.32) 
ETE 2TR sin (8) 








D. GROUND PLANE EFFECTS ON THE RADIATED SPACE WAVE 

The electric field equations derived in Chapters 3A-3C 
are free space equations which do not consider the presence of 
a ground plane. Since antenna orientation is described with 
respect to the ground plane, orientation descriptions are not 
really valid until a ground plane is introduced. The free 
space electric field distribution for a dipole oriented along 
the x axis or y axis is identical to the distribution for the 
same dipole along the z axis. The radiation pattern is 
displaced to correspond to the new position of the dipole, but 
the radiated electric field distribution relative to the 
antenna is unaffected. 

The free space electric field equations no longer provide 
an accurate description of an antenna's radiated electric 
field distribution above a ground plane. The free space 
equations only provide the direct wave component of the space 
wave which propagates directly from the antenna to the field 
points above ground (--7/2« 0« m/2). They are not valid for 


10| > 7/2 (points below the ground plane). 
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A complete description of the space wave above a ground 
plane must include the reflected wave component. The 
reflected component consists of electromagnetic waves radiated 
at elevation angles corresponding to |6| > 7/2 in the free 
space equations. These waves strike the ground plane and are 
absorbed, partially reflected, or completely reflected before 
reaching a field point above the ground plane. Trhestotal 
space wave is obtained by adding the direct wave and reflected 
wave components for field points defined by -7/2 < 0 « m/2. 
The distribution of the reflected waves vary relative to the 
antenna as its spatial orientation changes, so the electric 
field distribution (radiation pattern) is also affected when 
the antenna is repositioned with respect to the ground. 

The reflected wave component is determined using 
reflection coefficients and image theory. Reflection 
coefficients determine how much of an incident electric field 
is reflected by a given ground plane, and they are a function 
of @ and the complex index of refraction. The complex index 
of refraction (n) in terms of frequency, conductivity (0), and 


relative dielectric constant (¢€,) of the ground plane is 


The vertical reflection coefficient (T,) 


2 EN 2-a4n2 
p= 227608 (8,) -/n?-sin? (6,) — 
pocos r a Sin (0) 
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is used to find the reflected wave contribution to the @ 
component of the space wave. The horizontals shoei on 


Coeficiente 


r,- cos (0%) -/n*-sin?(0,) i SES 
cos (0,) «/n?-sin?(0,) 
is used to find the reflected wave contribution to the 6$ 
component of the space wave [Ref 5: pp. 629-633]. 

Antenna image theory accounts for phase differences 
between the direct and reflected wave components. Phase 
differences are caused by the different distances travelled by 
the two components to a given field point. The image is a 
'mirror image' located below the ground plane at a distance 
equal to the antenna height above ground. The distance from 
the image to a given field point is the distance travelled by 
the reflected wave to that field point. This distance is 
approximated by R+H,cos(@), similar in derivation to the 
approximation illustrated in Figure 4.4 for the direct wave 
distance of R-H,cos(@). The distances travelled by the two 
components differ by 2H,cos(@), but the direct wave expression 
already includes the term exp[+j8H,cos(@)], so the reflected 
wave expression includes the term exp[-jG6H,cos(@)] to account 
for any resulting phase differences. 

In equations 4.34 and 4.35, the variable 0, is the angle 
between the z axis and the radial from antenna image to field 
point. If the height of the antenna is small in comparison to 


the far-field distance, 0, is approximated as 0 (the angle 


T 
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SS Wecweeie Zea xis alee tie radial from origin to field point). 
The Mathcad applications use this approximation to find the 
Petlecwlon coerticienes tor 632 discrete values of 0 (316 for 
MMC EpIEcsmens)erHotHe interval -7/2 <0 < 1/2 at a 
specified distance (R) from the antenna. 

The direct wave component of the space wave is calculated 
for discrete values of 0 with the free space electric field 
equation and the phase terms derived from the antenna height. 
The reflected wave components are then calculated for each 6 
by multiplying the direct wave components by their 
corresponding reflection coefficients and applicable phase 
terms. The reflected wave components are then added to the 
direct wave components to obtain the total radiated space wave 
at each @. The discrete results for the space wave are 
normalized with respect to the value with maximum magnitude, 
and the space wave radiation pattern is depicted by plotting 
the normalized field magnitudes for -m/2 « 0 « m/2. 

The Mathcad applications plot radiation patterns assuming 
a planar surface below the antenna. The planar earth model 
provides accurate results for (1) antenna heights within 
several wavelengths of the surface and (2) wavelengths much 
smaller than the radius of the earth. The antennas addressed 
in this report are usually mounted well within several 
wavelengths of the ground, so the planar assumption in not a 
problem for the majority of calculations. If a configuration 


Pomencelitereee where the antenna is more than four or five 
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wavelengths above the ground, the Mathcad applications still 
provide a fair estimate of the radiation parameters. However, 
if highly accurate results are required, a computer program 
which accounts for the curvature of the earth should be used. 

The final expression for the radiated space wave of a 
vertical dipole antenna is [Ref 6: p. 166] 
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This is the equation used in the Elevated Vertical Dipole 
application to calculate the space wave radiation pattern and 
radiated electric field distribution. 
E. THE SURFACE WAVE 

The radiated electric field of an antenna above a ground 
plane includes a surface wave component in addition to the 
Space wave. The surface wave radiation pattern is a function 
of the distance from the antenna and the relative dielectric 
constant and conductivity of the ground plane. Since the 
surface wave attenuates exponentially with increasing range, 
the surface wave radiation pattern varies with distance from 
the antenna. While the space wave can be calculated at 
different distances from the antenna, its radiation pattern is 


predominately invariant with distance from the antenna. 
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EGuUaeLOns f£er the radiated electric fields of a vertical 
antenna above planar earth were derived by Sommerfeld in 1909 
[Ref 7], but the final form of these equations is complicated 
and difficult to evaluate. Norton simplified these equations 
to a form suitable for engineering applications [Ref 8]. 
Norton's expression for the surface wave of an elemental 


vertical dipole above a ground plane is [Ref 5: p. 644] 


-JP [R+H,cos (0) ] 
, (4.399 


jsoprat aT r2 E 


&a-un «in (9) [1325509 u y1-u?sin*(80) 


Its derivation is not explained here, but references 6 and 7 
give a detailed derivation. In equation 4.37, R is distance 
comte Eleld polne, I and k are unit vectors perpendicular and 
parallel to the antenna axis respectively, u is the inverse of 
the complex index of refraction (n), and F is the surface wave 
attenuation factor. The function F introduces a surface wave 
attenuation dependent upon distance, frequency, and tne 
paca Ome demeoarelmavemewhleh Ehe wave is travelling. 


The attenuation factor is given by [Ref 5: pp. 644-651] 
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where P, is the complex numerical distance for vertical 
polarization and the erfc term is the complementary error 
function defined by 


erfc(j/P,) si (4.41) 


Mathcad is incapable of evaluating the complementary error 
function for complex numbers, so an asymptotic expansion and 
infinite series approximation [Ref 9] was written into the 
Mathcad code to evaluate the erfc terms. 

When the magnitude of j: VP. is less than 2.18, the 


Mathcad applications use 





-x? 
erf(x+jy) =er£(x) += |(1-cos (2xy)) + j sin(2xy) (Ae) 
oo - +7? 
eae Yin E alo DU, Y RIOS Vel ya) where 
n mai n*+4x? 


£,(X, y) =2x-2x cosh(ny) cos (2xy) +n sinh(ny) sin(2xy) 
g, (x,y) =2x cosh(ny)sin(2xy) +n sinh (ny) cos (2xy) 


le(x,y)| = 107 |erf(x+Jy) | 


an infinite series approximation which evaluates erf(jVP,), 
where 'erf' is the standard error functiomes(Ret 9p 127 


The solution to erfc(jVP,) is then found by 


e 


erfc(x*jy)-1-erf(x-jy) (4.43) 
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When the magnitude of j-VP. iS greater than 2.18, the 


e 


Mathcad code uses 


IC GE C T EE v osse Ty een (4.44) 
m=1 ee)” 
cy EE] a 
m=1 m Pau 


(2-00, largz| < a 


an asymptotic expansion which solves for j-VaP, -erfc(jVP. ) 
[Ref 9 p. 298]. Either of these solutions can be substituted 
directly into equation 4.38 to find the surface wave 
attenuation factor. 
A magnitude of 2.18 is used as the transition point in the 
Mathcad applications because it provides the smoothest 
transition in both magnitude and phase between the series 
approximation and asymptotic expansion results. 

For antennas which radiate an electric field with a ¢ 
component, the surface wave attenuation factor for horizontal 


polarization (F,) must be found using 


F,-1-j/nP,e "erfc(j/P,) (4.45) 
and E UBI coo br os (a TO reg dg (0) (4546) 
2sin?(0,) 


where P, is the complex numerical distance for horizontal 
polarization. The variable 0, in equations 4.40 and 4.46 for 
P, and B,, respectively is approximated as 0 (just as it was 


for the reflection coefficients in equations 4.34 and 4.35). 
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The numerical distances (P, and Pa) and the surface wave 
attenuation factors (F, and F,) are evaluated for 0?s0s90? to 
determine surface wave radiation patterns. Their value is of 
most interest, however, directly at the surface (0-90?) where 
the magnitude of F is called the ground wave attenuation 
factor. Numerical distance is proportional to distance from 
the antenna and the square of the frequency, while it is 
approximately inversely proportional to ground conductivity. 
When the magnitude of the numerical distance is less than 
about 4.5, empirical formulas show that the ground wave 
attenuation factor varies exponentially with numerical 
distance. For numerical distance magnitudes greater than 4.5, 
empirical formulas show that the ground wave attenuation 
factor is inversely proportional to the numerical distance. 
Numerical distance magnitudes are typically much greater than 
4.5 3n the far 11era This implies that the ground wave 
(0290?) field strength varies inversely with the square of the 
distance from the antenna since numerical distance is directly 
proportional to distance from the antenna and equation 4.37 
shows an explicit 1/R dependence in addition to the 
attenuation factor. 

The Mathcad applications determine the complex numerical 
distances for the same discrete values of 0 described in 
Chapter 4.D. After determining the complementary error 
function for each numerical distance result, the surface wave 


attenuation factor is calculated for each discrete 4. The 


Bo 


surface wave 1s evaluated for each @ using the equations for 
Bieweteoceric ficld distribution, and the discrete results of 
the surface wave are normalized with respect to the value with 
maximum magnitude. The surface wave pattern is depicted by 
PloeeligGethia nOormalizedymagnitudes for -7/2 < 0 « m/2. The 
numerical distances (P, and Pa) at the surface (#=90°) are also 
listed with the predicted radiation parameters. 

With all of the terms defined, equation 4.37 is 
integrated over the length of a generic finite length vertical 
dipole to obtain the radiated surface wave distribution of an 


elevated vertical dipole [Ref 6 p.166] 





- JP [R+H¿cos (8) ] » 
prop EE c ol (4.47) 
(1 poss sinë (0) ALTO cos (0) 
n 





The total radiated far-field electric field distribution 
of an elevated vertical dipole above a planar earth can now be 


expressed by combining equations 4.36 and 4.37 to write 


e JP eos Tr eos[BA cos(0)] -cos (BA) 
ee |: 4.48 
o7 8 R sin (8) | 
; : enn 2 
o csi H (1-T.) ROO sin? (0) x Eos (0) 
n 











The first two terms inside the brackets represent the space 
Wave and the third term represents the surface wave. 
Chapter IV has defined the necessary terms and provided 


a d@taztled derivation of the total electric field distribution 
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for an elevated vertical dipole antenna. The result is a 
radiated electric field distribution with a @ component only. 
Most antenna types radiate an electric field with both 0 and 
$ components. While a detailed derivation of an electric 
field $ component was not presented here, the necessary 
horizontal polarization terms associated with deriving 6$ 
components were discussed. When these horizontal terms are 
substituted for their vertical counterparts, the method used 
to derive the electric field equation for the $ component of 
a radiated space wave or surface wave is identical to that 
used in deriving the 0 component for the elevated vertical 
dipole. 

The remaining chapters of this thesis explain the antenna 
configuration and the calculations of each individual Mathcad 
application. Only the final expressions for the radiated far- 
field electric field distributions are provided since the 
derivation of each follows the same procedure as explained in 
chis chapter. There is a lot of text repetition among the 
remaining chapters because each is written as a stand alone 
users guide for its associated Mathcad application. Any 
important concepts not presented in previous chapters are 


discussed as required. 
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V. THE ELEVATED VERTICAL DIPOLE 
The orientation of the elevated vertical dipole is 
depne cedin HT Gure Stieewicre 1941s the length of the dipole, 
Hy is the height of the feed above ground, R is the radial 
coordinate, 0 is the elevation coordinate, and ¢ is the 


azimata coordinate. 
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FIGURE 5.1: Spatial orientation of the elevated vertical 
dipole antenna for its corresponding Mathcad application. 
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The Mathcad application for the elevated vertical dipole 


requires the following inputs: 


jo EE length of the dipole (meters) 

ES MINE. Operational frequency (Hertz) 

Hodie antenna feed height above ground plane (meters) 
Re eee distance from antenna (meters) 

uM relative dielectric constant of ground plane 

ger UT conductivity of ground plane 


The length of the dipole is used to determine the frequencies 
which correspond to a A/4, A/2, 3A/4, or à length dipole. The 
user inputs the frequency and the radiation parameters 
discussed in Chapter 3 are calculated. Feed height (Hy) must 
be at least equal to the half-length of the dipole (h-1/2), 
and the distance from the antenna (R) must meet the far-field 
requirements of Chapter 3.C. 

The equation for the radiated electric field distribution 


of the elevated vertical dipole was derived in Chapter 4 as 


i. e JPI2-eos(9)) | cos [Bh cos (09) -cos (Bh) 
E1j601,—,———| e (5.1) 
-j2PH,¿cos (8) -j2BH,cos (8) a n?^-sin? (0) 
1+P,e PODE e SA 
n 


The first two terms inside the brackets represent the space 
wave and the third term represents the surface wave. I, is 
one since a sinusoidal current input with a maximum of unity 
is assumed. 

The requested inputs are used to calculate variables at 
the top of page 41 for 632 discrete values of 0 which are 


equally incremented from -7/2 to m/2. 
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para ao half-length of the dipole (meters) 


Nie. nie eee wavelength of the operational frequency (meters) 
o], NUM free space wavenumber for operational frequency 
DEMENS a distance from antenna to field point (meters) 

ce oe distance from antenna image to field point (meters) 
DEA ens index of refraction 

A toos complex numerical distance (vertical polarization) 
DNE Us vertical reflection coefficient 

NER us vertical surface wave attenuation factor 


The calculated variables are used to evaluate the 
radiated far-field space wave and surface wave for each 
discrete @. The space wave and surface wave results are 
combined for corresponding values of @ to obtain the total 
radiated electric field distribution. The space wave, surface 
wave, and total electric field results are then normalized 
with respect to the maximum field intensity of each, and the 
normalized magnitudes are plotted for each value of @ to 
depict the radiation patterns. The patterns are valid for any 
vertical plane containing the antenna axis, because the 
radiated electric field does not vary with changes in 6$. 

Equations 3.8 and 3.9 are used to integrate equation 5.1 
over the hemispherical Gaussian surface above the ground plane 
at a fixed radius (R) from the antenna to find total average 
radiated power (Pa). With the discrete values of the electric 
field and total average radiated power determined, the Mathcad 
application predicts the following radiation characteristics 


from the equations in Chapter 3: 


Ro ae a radiation resistance (Ohms) 

De w-. directivity 

BREUI effective radiated isotropic power (Watts) 

/ COS maximum theoretical effective area (square meters) 
— >.. maximum theoretical effective length (meters) 

PP e a numerical distance (vertical polarization, 0-290?) 


Angle... -elevation angle of maximum directive gain (degrees) 


41 


As an example, the Mathcad vertical dipole application 


was executed with the following inputs: 


dipole length 3 meters 
frequency 100-10? Hertz 
height of antenna feed 1.5 meters 
distance from the antenna 3000 meters 
relative dielectric constant 15 
ground conduc iv a SO 


Figures 5.2 and 5.3 depict the radiation patterns for the 


space wave and surface wave for this example. 
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FIGURE 5.2: Vertical dipole FIGURE 5.3: Vertical dipole 
space wave radiation pattern. surface wave pattern. 


The following radiation parameters were predicted by the 


Mathcad application for a sinusoidal current input of one Amp: 


Total power radiated 31.233 Watts 
Radiation resistance 62.467 Ohms 
Direceivaty 6.034 

Effective isotropic radiated power 188.470 Watts 
Maximum theoretical effective area 4.322 sq meters 
Maximum theoretical effective length 1.692 meters 
Numerical distance (vertical) 19 Seet7 9 

Elevation directi TR So degress 
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These results are consistent with expectations for this 
particular configuration: Appendix A contains computer 
hardcopies of additional example calculations for the vertical 
dipole and compares predicted radiation parameters to those 


expected based on previous calculations. 
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VI. THE VERTICAL MONOPOLE 
The orientation of the vertical monopole is depicted in 
Figure 6.1, where h is the length of the monopole, the antenna 
feed is at ground level, R is the radial coordinate, 0 is the 


elevation coordinate, and ġġ is the azimuth coordinate. 


e FIELD POINT 
(R, 0, Qœ) 


x axis 


Vertical Monopole Orientation 





FIGURE 6.1: Spatial orientation of the vertical monopole 
antenna for its corresponding Mathcad application. 
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The Mathcad application for the vertical monopole 


requires the following inputs: 


n" ee length of the monopole (meters) 

MM cr operational frequency (Hertz) 

E a ió distance from antenna (meters) 

E oss relative dielectric constant of ground plane 
ee cn Conductivity of ground plane 


The length of the monopole is used to determine frequencies 
which correspond to a A/8, A/4, 3/8, or A/2 long monopole. 
The user inputs the frequency for which the radiation 
parameters discussed in Chapter 3 are calculated. Distance 
from the antenna (R) must meet the far-field requirements of 
Ghepter 3 .C: 

The monopole's electric field distribution is derived in 
the same way as the vertical dipole equation was derived in 


Chapter 4 and is given by [Ref 6: p.152] 








^ eJPR| A«JB A-jB 
AO sin(0) ^ "sin(0) n 
TENES LA-JBS sin? (9) -Y2^-sin*(6) os (8) 
= Sano) n? 


UL sin (ph (sin? (8 pies nal Neos e| n?-sin?(0) | 
n 


nessuno) 
A=cos [Bh cos (8) } -cos (BA) 


B=sin[Bh cos (8) ] -cos (8) sin( BA) 


The first two terms inside the brackets represent the space 
wave and the second two terms represent the surface wave. The 
vertical monopole's far-field electric field has a @ component 
only , just as the vertical dipole. I, iS one since a 


Sinusoidal current input with a maximum of unity is assumed. 
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The requested inputs are used to calculate the following 
variables for 632 discrete values of 0 which are equally 


incremented from -T/2 to m/2: 


A E E wavelength of the operational frequency (meters) 
Bo R free space wavenumber for operational frequency 
D. 5 IG: index of refraction 


complex numerical distance (vertical polarization) 
vertical reflection coefficient 
vertical surface wave attenuation factor 

The calculated variables are used to evaluate the 
radiated far-field space wave and surface wave for each 
discrete l. The space wave and surface wave results are 
combined for corresponding values of @ to obtain the total 
radiated electric field distribution. The space wave, surface 
wave, and total electric field results are then normalized 
with respect to the maximum field intensity of each, and the 
normalized magnitudes are plotted for each value of @ to 
depict the radiation patterns. The patterns are valid for any 
vertical plane containing the antenna axis, because the 
radiated electric field does not vary with changes in q. 

Equations 3.8 and 3.9 are used to integrate equation 6.1 
over the hemispherical Gaussian surface above the ground plane 
at a fixed radius (R) from the antenna to find the total 
average radiated power (Paa). With the discrete values of the 
electric field and total average radiated power determined, 
the Mathcad application predicts the radiation characteristics 
listed at the top of the next page from the equations in 


Chapter 3. 
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Radiation characteristics determined by the vertical 


monopole application are: 


oe radiation resistance (Ohms) 

DOM" E. directivity 

¡A M. effective radiated isotropic power (Watts) 
D-—. maximum theoretical effective area (square meters) 
AA maximum theoretical effective length (meters) 


P, ..... numerical distance (vertical polarization, 60-90?) 
Angle, .elevation angle of maximum directive gain (degrees) 


As an example, the Mathcad monopole application was 
executed with the following inputs: 


monopole length 0.5 meters 
frequency 150-10? Hertz 
distance from the antenna 3000 meters 
relative dielectric constant T5 
ground conductivity 5-10? 


Figures 6.2 and 6.3 depict the radiation patterns for the 


Space wave and surface wave for this example. 





FIGURE 60.2: Monopole space FIGURE 6.3: Monopole surface 
wave radiation pattern. wave radiation pattern. 


The following radiation parameters were predicted by the 


Mathcad application for a sinusoidal current input of one Amp: 
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Total power radiated 5.607 Watts 


Radiation resistance 11.215 ORMS 
Directivity 3.246 

Effective isotropic radiated power 18.203 Watts 
Maximum theoretical effective area 1.033 sq meters 
Maximum theoretical effective length 0.351 meters 
Numerical distance (vertical) 293.016 

Elevation of directivity 26.815 degrees 


These results are consistent with expectations for this 
Daktdcular VeowmiguLrarren- Appendix B contains computer 
hardcopies of additional example calculations for the vertical 
monopole and compares predicted radiation parameters to those 


expected based on previous calculations. 
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VII. THE HORIZONTAL DIPOLE 
ii eoh levitate Cmmom=e te Get zZeoncal dipole is depicted in 
Figure 7.1, where 1 is length of the dipole, H, is the height 
of the feed above ground, R is the radial coordinate, Ü is the 


elevation coordinate, and $ is the azimuth coordinate. 


e Eom SIS ee © Nii. 
(R, O, o) 


Beruomitalsbapole Orientation 





FIGURE 7.1: Spatial orientation of the horizontal dipole 
antenna for its corresponding Mathcad application. 
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The Mathcad application for the horizontal dipole requires 


the following inputs: 


jeu aT X length of the dipole (meters) 

Hee. 2. . antenna feed height above ground plane (meters) 
ae eee operational frequency (Hertz) 

RON MN E distance from antenna (meters) 

C. UEM relative dielectric constant of ground plane 
oom i conductivity of ground plane 

A NE elevation angle index (from Table 3.1) 


The length of the dipole is used to determine frequencies 
which correspond to a A/4, A/2, 3A/4, or A length dipole. The 
user inputs the frequency for which the radiation parameters 
discussed in Chapter 3 are computed. Feed height, H,, must be 
greater than or equal to zero, and distance from the antenna 
(R) must meet the far-field requirements of Chapter 3.C. The 
elevation angle index (d) sets the elevation coordinate, 60, 
for which a horizontal radiation pattern is determined. Table 
3.1 lists possible indices and their corresponding elevation 
angles from .285° to 88.8579 In increments Aor mE Donn e 
Indices between those listed can be used to interpolate a 
better approximation of a desired elevation. 

The horizontal dipole's radiated electric field is a 
combination of 6 and $ components. In addition, unlike the 
vertical dipole, its radiation pattern varies with changes in 
$. The total radiated electric field is the vector sum of the 
ô and $ components. The electric field for the horizontal 
dipole is obtained in a manner analogous to that used for the 
vertical dipole discussed in Chapter 4. The horizontal 


dipole's § and $ field components are given by [Ref 6: p. 144] 
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E,2-j60I,cos (o) [sos tbh sini8i conti] -cos (BA | (7.2) 


1-sin? (8) cos? (@) 


-jB [R-H, cos (8) ] " 
= cos (0) (1-T,e polo 


R 
n ES) eee ten pS (8) [sin (8) E asii (0) e (B) | 
n 11 


ms 


g,=J607,sin(o)( SS (Bh sin (6) cos (9) “cos (Bh) | 


1-sin? (6) cos? () 
-jP [R-H,cos (8) ] - 

e E [1 Te j2PH¿cos (8) i (1-D,) Fe jode d 

The first two terms inside the brackets of each equation are 

the space wave components and the third terms are the surface 

wave components. I, is taken to be one, since a sinusoidal 

current input with a maximum of unity is assumed. 

The requested inputs are used to calculate the following 
variables using a constant $ of $20 and $-7/2 for 632 discrete 
values of 0 which are equally incremented from -7/2 to m/2: 

[p ee half-length of the dipole (meters) 
a v wavelength of the operational frequency (meters) 
eer free space wavenumber for operational frequency 
P emes distance from antenna to field point (meters) 
R distance from antenna image to field point (meters) 
A index of refraction 
P complex numerical distance (vertical polarization) 
P complex numerical distance (horizontal polarization) 
DTP A vertical reflection coefficient 
IE eo. horizontal reflection coefficient 
F vertical surface wave attenuation factor 
F horizontal surface wave attenuation factor 

The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 for $0 and ¢=7/2. The total space wave and surface wave 


distributions are determined from vector addition of 


x 


corresponding ô and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of @ to obtain the total radiated electric field 
distribution for the $-0 and $-7/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 0 to depict the radiation patterns. The Mathcad 
horizontal dipole application computes the space wave, surface 
wave, and total electric field radiation patterns and 
radiation parameters in the $-0 and $-7/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0, as $ varies from O to 2m in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the 6 and $ electric field 
components to the total average radiated power (Paa), equations 
3.8 and 3.9 are used to integrate equations 7.1 and 7.2 over 
the hemispherical Gaussian surface above the ground plane at 
a fixed radius (R) from the antenna. With the discrete values 
of the electric field and total average radiated power, the 
Mathcad application predicts the following radiation 


characteristics from the equations in Chapter 3: 
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R radiation resistance (Ohms) 

AE Cite civil 

EIR Poa effective radiated isotropic power (Watts) 

A ee maximum theoretical effective area (square meters) 

SR Ls maximum theoretical effective length (meters) 

PME es numerical distance (vertical polarization, 0290?) 

PE numerical distance (horizontal polarization, 0-90?) 

Angle,, .elevation angle of maximum directive gain (degrees) 
The directivity (Dj, effective isotropic radiated power 
(EIRP), maximum theoretical effective area (Anax), maximum 


theoretical effective length (1,,), and the elevation angle of 
maximum directive gain (Angle,,) are all determined for both 
the $20 and ¢=7/2 vertical planes. 

As an example, the Mathcad horizontal dipole application 


was executed with the following inputs: 


dipole length 1.0 meter 
frequency 150-10° Hertz 
height of antenna feed 1.0 meter 
distance from the antenna 3000 meters 
relative dielectric constant Ilo 

ground conductivity 1o 
elevation angle index EE 


The following radiation parameters were predicted by the 


Mathcad application for a sinusoidal current input of one Amp: 


Total power radiated (Watts) 26.116 
Radiation resistance (Ohms) 52.233 
Numerical distance (vertical) 292.419 
Numerical distance (horizontal) 66215.4 

=0 Si 2 
Directivity 253550 gH 215 
EIRP (Watts) pau 2 188.435 
Max eff area (sq meters) 0.417 Dc 
Max eff length (meters) 0.481 jS 
Angle,, (degrees) 45.071 28.241 


Figures 7.2 through 7.7 are the space wave and surface 
wave radiation patterns in the $-0 and $-7/2 vertical planes 


and the designated horizontal plane for this example. 
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FIGURE 7.2: Horizontal 


dipole space wave radiation 
pattern for $20 vertical 
plane. 


o NA, APA 1, an | 
FIGURE 7.4: Horizontal 
dipole space wave radiation 
pattern for g¢=m7/2 vertical 
plane. 
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PIGURE 73: 


dipole surface wave radiation 
pattern for $z0 vertical 
plane. 
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dipole surface wave radiation 
pattern for dqQ=r/2 vertical 
plane. 
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FIGURE 7.6: FIGURE 7.7: Horizontal 


dipole space wave radiation dipole surface wave radiation 
pattern for horizontal plane. pattern for horizontal plane. 


These results are consistent with expectations for this 
particular configuration. Appendix C contains computer 
hardcopies of additional example calculations for the 
horizontal dipole and compares predicted radiation parameters 


to those expected based on previous calculations. 
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VIII. THE ARBITRARILY ORIENTED DIPOLE 
The orientation of the arbitrarily snventedma1 ole mas 
depicted in Figure 8.1, where 1 is length ofitie dipole, BT- 
the feed height above ground, « is the angle of the dipole 
axis with the horizontal, R is the radial coordinate, 0 is the 


elevation coordinate, and $ is the azimuth coordinate. 


FIELD POINT 
* (R,0,0) 


x X axis 
Arbitrarily Oriented Dipole Orientation 





FIGURE 8.1: Spatial orientation of the arbitrarily 
oriented dipole antenna for its Mathcad application. 
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tHe Mathead application for the arbitrarily oriented 


dipole requires the following inputs: 


DE coh length of the dipole (meters) 

Oei TES antenna feed height above ground plane (meters) 
GUESS riis E angle of antenna axis with the horizontal (degrees) 
LC MM operational frequency (Hertz) 

RE n c E distance from antenna (meters) 

cu... relative dielectric constant of ground plane 

DNE D soe conductivity of ground plane 

OU NEN E elevation angle index (from Table 3.1) 


The length of the dipole is used to determine frequencies 
which correspond to a A/4, A/2, 3A/4, or À length dipole. The 
user inputs the frequency for which the radiation parameters 
discussed in Chapter 3 are computed. Feed height, H,, must be 
greater than or equal to (1/2):sin(o), and the distance from 
the antenna (R) must meet the far-field requirements of 
Chapter 3.C. The elevation angle index (d) sets the elevation 
coordinate, @, for which a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857° in 
increments of about 2°. Indices between those listed can be 
used to interpolate a better approximation of a desired 
elevation. 

The radiated electric field for the arbitrarily oriented 
dipole is a combination of § and $ components, and its 
radiation pattern varies with changes in q. The total 
electric field is the vector sum of the 6 and 9 components. 
The arbitrarily oriented dipole's radiated electric field is 
obtained in a manner analogous to that used for the vertical 


dipole discussed in Chapter 4, and is given by 
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E,--j601, 5 e E (8.1) 
A, | | | 
————— ——— (cos (a)sin($)cos(0)-sin(«)sin(9)) 
1-cos?(w) 
-p,— — 22... (eos (a) sin (6) cos (8) -sin(a) a 
1-cos?(w) 
PA c enne! A (8)- ¥2*-sin’ (8) e e| 
1-sin?(w) n? 
(eos (a) sino) LEE sin(a) sinto] 
n 
, e Pk j8H,cos (8) 
E,7 j601, e^t cos (a)cos (o) (8.2) 
A, A, -jJ2BH,cos (9) T | 
$$ 4 ———À— 2 (gr qs 
I 1-cos? (p’) j Cat nl Em) 


where A,=cos [Bh cos(p)] -cos (Bh) 
cos (y) =cos (8) sin(a) +sin(8)cos(a)sin(@) 
A,=cos [Bh cos (p’) ] -cos (Bh) 


cos (’) =sin(@) cos (a) sin() -cos (8) sin(a) 


Equation 8.1 is the 6 component (vertically polarized), and 
equation 8.2 is the ¢ (horizontally polarized) component [Ref 
6: p. 174-175]. The first two terms inside the brackets of 
each equation represent the space wave components and the 
third terms represent the surface wave components. Io is 
taken to be one, since a sinusoidal current input with a 


maximum of unity is assumed. 
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The requested inputs are used to calculate the following 
variables using a constant $ of $-0 and $-T/2 for 632 discrete 


values of 0 which are equally incremented from -7/2 to m/2: 


Di ÓN half-length of the dipole (meters) 

A on. wavelength of the operational frequency (meters) 
D. free space wavenumber for operational frequency 

Ro f PM distance from antenna to field point (meters) 
RE i. distance from antenna image to field point (meters) 
ic NR index of refraction 


P complex numerical distance (vertical polarization) 
P complex numerical distance (horizontal polarization) 
in vertical reflection coefficient 
E MEET horizontal reflection coefficient 
F vertical surface wave attenuation factor 
F horizontal surface wave attenuation factor 

The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 for $20 and $-7/2. The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding § and 4 components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of 0 to obtain the total radiated electric field 
distribution for the $20 and ¿=1/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete @ to depict the radiation patterns. The arbitrarily 
oriented dipole Mathcad application computes the space wave, 
surface wave, and total electric field radiation patterns and 
radiation parameters in the $-0 and ¢=n/2 vertical planes. 


The variables corresponding to the selected elevation 


angle index (d) are used to evaluate the radiated electric 
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field components for the space wave and surface wave at the 
fixed elevation angle (0, as $ varies from O to 2m in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those in the vertical planes. 

To find the contributions of the @ and $ electric field 
components to the total average radiated power, P.,, equations 
3.8 and 3.9 are used to integrate equations 8.1 and 8.2 over 
the hemispherical Gaussian surface above the ground plane at 
a fixed radius (R) from the antenna. The sum of the two 
integrals is the antenna's P.j. With the discrete values of 
the electric field and total average radiated power 
determined, the Mathcad application predicts the following 


radiation characteristics from the equations in Chapter 3: 


eames” «ae radiation resistance (Ohms) 

Dar T directivity 

Bgl Pee effective radiated isotropic power (Watts) 

a eee maximum theoretical effective area (square meters) 

les ceo cus maximum theoretical effective length (meters) 

rn S numerical distance (vertical polarization, 60-90?) 

Paes numerical distance (horizontal polarization, 0-909?) 

Angle,, .elevation angle of maximum directive gain (degrees) 
The directivity (D), effective isotropic radiated power 
(EIRP), maximum theoretical effective area (A,,,), maximum 


theoretical effective length (lax), and the elevation angle of 
maximum directive gain (Angle,,) are all determined for both 


the $20 and $-7/2 vertical planes. 
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IX. THE INVERTED L ANTENNA 
The orientation of the inverted L antenna is depicted in 
Figure 9.1, where H is length of the vertical antenna segment, 
l is the length of the horizontal antenna segment, the feed is 
at ground level, R is the radial coordinate, @ is the 


elevation coordinate, and $ is the azimuth coordinate. 
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FIGURE 9.1: Spatial orientation of the inverted L antenna 
for its corresponding Mathcad application. 
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The Mathcad application for the inverted L antenna 


requires the following inputs: 


| A length of vertical antenna segment (meters) 

Jo dU esos length of horizontal antenna segment (meters) 
BE EIE Operational frequency (Hertz) 

Re Ud E distance from antenna (meters) 

CA a a relative dielectric constant of ground plane 
QT. EN conductivity of ground plane 

Cj M elevation angle index (from Table 3.1) 


The combined length of the vertical and horizontal antenna 
segments (H+l) is used to determine frequencies which 
correspond to a A/8, -A/4, 93A/8, 40505472 ond verces E 
antenna. The user inputs the frequency for which the 
radiation parameters discussed in Chapter 3 are computed. 
Distance from the antenna (R) must meet the far-field 
requirements of Chapter 3.C. The elevation angle index (d) 
sets the elevation coordinate, 0, for which a horizontal 
radiation pattern is determined. Table 3.1 lists possible 
indices and their corresponding elevation angles from .285? to 
88:857? in increments of abone k Indices between those 
listed can be used to interpolate a better approximation of a 
desired elevation. 

The inverted L’s radiated electric field is a varying 
combination of Ü and $ components, and its radiation pattern 
varies with changes in $. The total radiated electric field 
is the vector sum of the ÓÜ and $ components. The radiated 
electric field for the inverted L antenna is obtained in a 
manner analogous to that used for the vertical dipole 


discussed in Chapter 4. The radiated fields of each antenna 
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segment are considered separately, and superposition is used 
to determine the total radiated electric field. The vertical 
segment radiates only a § field component, but the horizontal 
segment radiates both § and $ components. The electric field 


equations for the inverted L are [Ref 6: p. 169-170] 
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Equation Jc = ĝ component for the vertical segment, 
equation 9.2 is the - (vertically polarized) component for the 
horizontal segment, and equation 9.3 is the $ (horizontally 
polarized) component for the horizontal segment. The first 
two terms inside the brackets of each equation are the space 
wave components and any additional terms are the surface wave 
components. I, is taken to be one, since a Sinusoidal current 
input with maximum of unity is assumed. 

The requested inputs are used to calculate the following 
variables using a constant q of $-0 and ¢=m/2 for 632 discrete 


values of 0 which are equally incremented from -7/2 to m/2: 


Discs N half-length of the dipole (meters) 

ee wavelength of the operational frequency (meters) 

v — free space wavenumber for operational frequency 
EE I. distance from antenna to field point (meters) 
Ron dc distance from antenna image to field point (meters) 


EM o index of refraction 

Doc complex numerical distance (vertical polarization) 
complex numerical distance (horizontal polarization) 
vertical reflection coefficient 

ee ens horizontal reflection coefficient 

OS vertical surface wave attenuation factor 
horizontal surface wave attenuation factor 


o 


O 


The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 for $20 and $-7/2. The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding § and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of @ to obtain the total radiated electric field 
distribution for the $20 and $sT/2 vertical planes. The space 


wave, surface wave, and total electric field results are then 
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normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 0 to depict the radiation patterns. The inverted L 
Mathcad application computes the space wave, surface wave, and 
Poca lecture tela radiation patterns and radiation 
parameters in the $20 and $-7/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0, as $ varies from O to 2m in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the Ó and $ electric field 
components to the total average radiated power, Paa, equations 
3.8 and 3.9 are used to integrate equations 9.1, 9.2, and 9.3 
over the hemispherical Gaussian surface above the ground plane 
at a fixed radius (R) from the antenna. The sum of the 
integrals is the antenna’s P.,. With the discrete values of 
the electric field and total average radiated power 
determined, the Mathcad application predicts the following 


radiation parameters using the equations from Chapter 3: 


io ES. radiation resistance (Ohms) 

Became directivity 

ETRE M. effective isotropic radiated power (Watts) 

OR sonis maximum theoretical effective area (square meters) 
EE maximum theoretical effective length (meters) 

EA numerical distance (vertical polarization, 0-90?) 
DNE es numerical distance (horizontal polarization, 0-909?) 


.elevation angle of maximum directive gain (degrees) 
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The directivity (Dj, effective isotropic radiated power 
(EIRP), maximum theoretical effective area (A,,), maximum 
theoretical effective length (1,,,), and the elevation angle of 
maximum directive gain (Angle,,,) are all determined for both 


the $20 and $-7/2 vertical planes. 
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X. THE SLOPING LONG-WIRE ANTENNA 
UtUVmOgreitdpEet Of tie ssie@ping long-wire antenna is 
depicted in Figure 10.1, where 1 is length of the antenna, the 
feed is at ground level, o is the angle of the antenna with 
NNUS sa 1Sutne radial coordinate, 6 is the elevation 


SeOnrGtiace, ana @ 15 the azimuth coordinate. 


e FIELD POINT 
(R,0,0) 


X axis 


y 


SiS AROS Mie Orientation 





FIGURE 10.1: Spatial orientation of the sloping long-wire 
antenna for its corresponding Mathcad application. 
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The Mathcad application for the sloping long-wire antenna 


requires the following inputs: 


DE e a E length of the antenna (meters) 

Oon angle of antenna with horizontal (degrees) 
BET. MM operational frequency (Hertz) 

Y —P— distance from antenna (meters) 

E. ARCAS relative dielectric constant of ground plane 
E conductivity offground plane 

E EE EE elevation angle index (from Table 3.1) 


The user inputs the frequency (f,) for which the radiation 
parameters discussed in Chapter 3 are computed. The distance 
from the antenna (R) must meet the far-field requirements of 
Chapter 3.C. The elevation angle index (d) sets the elevation 
coordinate, 0, for which a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857° in 
increments of about 2°. Indices between those listed can be 
used to interpolate a better estimate of a desired elevation. 

The electric field equations for the sloping long-wire 


are given by [Ref 6: p. 180-181] 
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e IBE 





Erda cos(«)cos(o) diia 
E any oI 
1-cos? (y) 1-cos? (y!) i = 


where A,=cos[Bl cos(wp)]-cos(B1) 
B,=s1nT[p1 cos (y)]-sin(fP1)cos(w) 
cos (y) »cos(0) sin(a) *sin(0)cos(a)sin(o) 
A,-cos[B1 cos(w^)]-cos(Q1) 
B,-sin([Bl cos(y)]-sin(pl)cos (wy) 
cos(y^) »sin(0)cos(a)sin(o$)-cos(0)sin(a) 
The radiated electric field for the sloping long-wire is 22 
obtained in a manner analogous to the vertical dipole equation 
discussed in Chapter 4. The long-wire's radiated electric 
field is a varying combination of § and $ components, and its 
radiation patterns vary with changes in $. The total radiated 
electric field is the vector sum of the § and ¢ components. 
The first two terms inside the brackets of each equation are 
the space wave components and the additional terms are the 
surface wave components. I, is taken to be one since a 
Sinusoidal current input with a maximum of unity is assumed. 
Requested inputs are used to compute the variables on the 
next page using a constant $ of $20 and $2m7/2 for 632 discrete 


values of 0 which are equally incremented from --/2 to m/2. 


A wavelength of the operational frequency (meters) 
ees see free space wavenumber for operational frequency 
TS. index of refraction 


P complex numerical distance (vertical polarization) 
P complex numerical distance (horizontal polarization) 
ls vertical reflection coefficient 

Deae horizontal reflection coefficient 

ree cea s vertical surface wave attenuation factor 

F horizontal surface wave attenuation factor 
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The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 for $20 and $-7/2. The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding Ó and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of @ to obtain the total radiated electric field 
distribution for the $20 and $-7/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete @ to depict the radiation patterns. The sloping 
long-wire Mathcad application computes the space wave, surface 
wave, and total electric field radiation patterns and 
radiation parameters in the @=0 and $-7/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0, as $ varies from O to 2mw in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the Ó and electric dield 
components to the total average radiated power, Paa, equations 
3.8 and 3.9 are used to integrate equations 10.1 and 10.2 over 


the hemispherical Gaussian surface above the ground plane at 
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a fixed radius (R) from the antenna. The sum of the integrals 
is the antenna’s P.,. With the discrete values of the electric 
field and total average radiated power determined, the Mathcad 
application predicts the following radiation parameters using 


the equations from Chapter 3: 


Bes Wee. radiation resistance (Ohms) 

peces directivity 

EN. effective isotropic radiated power (Watts) 

D sos ET maximum theoretical effective area (square meters) 

LL maximum theoretical effective length (meters) 

O idol numerical distance (vertical polarization, 0-290?) 

geo ues numerical distance (horizontal polarization, 0-290?) 

Angle,, .elevation angle of maximum directive gain (degrees) 
The directivity (Dj, effective isotropic radiated power 
(EIRP), maximum theoretical effective area (Amx), maximum 


theoretical effective length (1,,), and the elevation angle of 
maximum directive gain (Angle,,,) are all determined for both 


the $20 and ¢=m/2 vertical planes. 
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XI. THE TERMINATED SLOPING V ANTENNA 

The orientation of the terminated sloping V antenna is 
depicted on the next page in Figure 11.1, where l is inet o 
the antenna wires (numbered 1 and 2), H, is the feed height 
above ground, H, is the height of the termination points above 
ground, y is the angle subtended by the wires, yn is the angle 
between the x axis and the projection of the antenna legs in 
the x-y plane, o is the angle between the plane which contains 
the antenna and the x-y plane, R is the radial coordinate, Q 
is the elevation coordinate, and $ is the azimuth coordinate. 
The Mathcad application for the sloping V antenna requires 


the following inputs: 


lA length of the antenna legs (meters) 

El rcu ceu feed height above ground (meters) 

He eee height of termination points above ground (meters) 
Vee eure angle subtended by the antenna legs (degrees) 

f.. TEE Operational frequency (Hertz) 

—€————— distance from antenna (meters) 

en nee relative dielectric constant of ground plane 

pr i conductivity of ground plane 

ec eed em elevation angle index (from Table 3.1) 


The user inputs the frequency (f,) for which the radiation 
parameters discussed in Chapter 3 are computed. The feed 
height (Hj and termination height (Hj) can be any value 
greater than zero. The angle o is found from Hy, H,, wire 
length (1), and the half-angle between the wires (y/2) by 
H-H, 


Icos(1] 2m 


OS ie 


12 
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FIGURE 11.1: Spatial orientation of the terminated 
sloping V antenna for its corresponding Mathcad 
application. 








The angle y is then found from o and y/2 by 


tan X) 
n=tanY 124) er 
cos (a) 


Distance from the antenna (R) must meet the far-field 
requirements of Chapter 3.C. The elevation angle index (d) 
sets the @ coordinate for which a horizontal radiation pattern 
is determined. Table 3.1 lists possible indices and their 


corresponding elevation angles from .285° to 88.8579? in 


T3 


increments of about 2°. Indices between those listed can be 
used to interpolate a better estimate of a desired elevation. 

The radiated electric field for the sloping V antenna is 
obtained in a manner analogous to that used for the vertical 
dipole discussed in Chapter 4. The sloping V electric field 


equations are given by [Ref 6: p. 185-1860] 
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The sloping V’s radiated electric field is a varying 
combination of Ó and 4 components, and its radiation pattern 
varies with changes in $. The total radiated electric field 
is the vector sum of the @ and $ components. The fields of 
each antenna leg are considered separately, and superposition 
is used to determine the total electric field. Equations 11.3 
and 11.4 are the radiated Ó and $ components, respectively for 
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number one, and equations 11.5 and 11.6 are the 6 and à 
components, respectively for wire number two. The first two 
terms inside each set of brackets represent the space wave 
components and any additional terms represent the surface wave 
components. I, is taken to be one, since a sinusoidal current 
input with a maximum of unity is assumed. 

The requested inputs are used to calculate the following 
variables using a constant ¢ of ¢=0 and ¢=7/2 for 632 discrete 
values of 0 which are equally incremented from -7/2 to m/2: 

Ko cmon wavelength of the operational frequency (meters) 
EN UEM free space wavenumber for operational frequency 
mc Lus index of refraction 
a e complex numerical distance (vertical polarization) 
esc complex numerical distance (horizontal polarization) 
21 vertical reflection coefficient 
pM E horizontal reflection coeficiente 


enn eee vertical surface wave attenuation actor 
horizontal surface wave attenua llantos 


< o 


a 


trj Nj 5g dd tU tgp oo 


The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 for $20 and $-7/2. The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding § and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of @ to obtain the total radiated electric field 
distribution for the $20 and $-7/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 


discrete 0 to depict the radiation patterns. The sloping V 
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Mathcad application computes the space wave, surface wave, and 
Beedle eGerte field | radiation patterns and radiation 
parameters in the ¢=0 and ¢=n/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (@,) as ¢ varies from 0 to 2m in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the 6 and @ electric freld 
components to the total average radiated power, P „a, equations 
3.8 and 3.9 are used to integrate equations 11.3, 11.4, 11.5, 
and 11.6 over the hemispherical Gaussian surface above the 
ground plane at a fixed radius (R) from the antenna. The sum 
of the integrals is the antenna’s P,,. With the discrete 
values of the electric field and total average radiated power 
determined, the Mathcad application predicts the following 


radiation parameters using the equations from Chapter 3: 


ae radiation resistance (Ohms) 

ph iu directivity 

¡A effective isotropic radiated power (Watts) 

D maximum theoretical effective area (square meters) 

IMS E. ee maximum theoretical effective length (meters) 

pr cun M numerical distance (vertical polarization, 0-909?) 

p numerical distance (horizontal polarization, 0-90?) 

Angle, .elevation angle of maximum directive gain (degrees) 
The directivity (Dj, effective isotropic radiated power 
(EIRP), maximum theoretical effective area (A,,,), maximum 


M 


theoretical effective length (lax), and the elevation angle of 
maximum directive gain (Angle»x) are all determined for both 


the $-0 and $-7/2 vertical planes. 
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XII. THE SIDE-LOADED VERTICAL HALF RHOMBIC ANTENNA 

The orientation of the side-loaded vertical half rhombic 
antenna is depicted in Figure 12.1, where 1 is length of the 
antenna wires (numbered 1 and 2), the feed is at ground level, 
y is the angle between the antenna wires and the horizontal 
plane, R is the radial coordinate, 0 is the elevation 


coordinate, and $ is the azimuth coordinate. 


Fee POINT 
(R, 8,4) * 


m xr 


Side-Loaded Vertical Half-Rhombic Orientation 





FIGURE 12.1: Spatial orientation of the side-loaded 
vertical half rhombic antenna for its Mathcad application. 


79 


The Mathcad application for the vertical half rhombic 


antenna requires the following inputs: 


EE length of the antenna legs (meters) 

AE angle of antenna wires with horizontal (degrees) 
P. ieee operational frequency (Hertz) 

Ra. 522. distance from antenna (meters) 

C eas coe relative dielectric constant of ground plane 
S.A conductivity of ground plane 

e TN elevation angle index (from Table 3.1) 


The user inputs the frequency (f, for which the radiation 
parameters discussed in Chapter 3 are computed. Distance from 
the antenna (R) must meet the far-field requirements of 
Chapter 3.cC. The elevation angle index (d) sets the @ 
coordinate for which a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857° in 
increments of about 2°. Indices between those listed can be 
used to interpolate a better approximation of a desired 
elevation. 

The vertical half rhombic’s radiated electric field is a 
combination of 6 and ¢ components, and its radiation pattern 
varies with changes in $. The total radiated electric field 
is the vector sum of the Ó and $ components. The radiated 
electric field for the vertical half rhombic antenna is 
obtained in a manner analogous to that used for the vertical 
dipole discussed in Chapter 4. The radiated fields of each 
antenna leg are considered separately, and superposition is 
used to determine the total electric field. The electric 


field equations for the half rhombic are [Ref 6: p. 198-200] 
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where U, and U, are given by 
DST COSA BASICOS (in) 
cos (¥,) =cos (0) sin(@) +san (8) ces (a) cos) 
cos (P,) =sin (8) cos (a) cos ($) -cos (8) sin(a@) 


Equations 1231" andel27 20are seme 0 (vertically polarized) and 
$ (horizontally polarized) components, respectively for leg 
number one, and equations 12.3 and 12.4 are the § and $ 
components, respectively for leg number two. The first two 
terms inside each set of brackets are the space wave 
components and any additional terms are the surface wave 
components. I, is taken to be one, since a sinusoidal current 
input with a maximum of unity is assumed. The requested 
inputs are used to calculate the following variables using a 
constant $ of $20 and ¢=m7/2 for 632 discrete values of 0 which 


are equally incremented from -7/2 to m/2: 


Noc M wavelength of the operational frequency (meters) 
DIEN CERE free space wavenumber for operational frequency 

NN e nen index of refraction 

p.ex complex numerical distance (vertical polarization) 
A complex numerical distance (horizontal polarization) 
A vertical reflection coefficient 

I RT as horizontal reflection coefficient 

RUNE o vertical surface wave attenuation factor 

ID... horizontal surface wave attenuation factor 


The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 for 4-20 and $--/2. The total space wave and surface wave 
distributions are determined from vector addition of 


corresponding Ó and $ components. The space wave and surface 
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wave results are then combined for corresponding discrete 
values of 10 to obtain the total radiated electric field 
distribution for the $¢=0 and ¿=1/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
Mtro co orde pre the radiation patterns. The vertical 
half rhombic Mathcad application computes the space wave, 
surface wave, and total electric field radiation patterns and 
radiation parameters in the $-0 and $-7/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0, as $ varies from O to 2m in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the @ and DEC CCS 
Eoussmu-nscascuPncapotalgaverage radiated power, P,,, equations 
3.8 and 3.9 are used to integrate equations 12.1, 12.2, 12.3, 
and 12.4 over the hemispherical Gaussian surface above the 
ground plane at a fixed radius (R) from the antenna. The sum 
of the integrals is the antenna's BP. With the discrete 
values of the electric field and total average radiated power 
determined, the Mathcad application predicts the following 


radiation parameters using the equations from Chapter 3: 


ag 


pe... radiation resistance (Ohms) 

Do cmm directivity 

Bak? . uu effective isotropic radiated power (Watts) 

Dl eee maximum theoretical effective area (square meters) 
ee. maximum theoretical effective length (meters) 

D: o ee numerical distance (vertical polarization, 602909?) 


Be. —À— numerical distance (horizontal polarization, 0-909?) 
Angle,, .elevation angle of maximum directive gain (degrees) 


The directivity (Dj, effective isotropic radiated power 
(EIRP), maximum theoretical effective area (A,,,), maximum 
theoretical effective length (lax), and the elevation angle of 
maximum directive gain (Angle,,,) are all determined for both 


the $20 and $-7/2 vertical planes. 
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XIII. THE TERMINATED SLOPING OR HORIZONTAL RHOMBIC ANTENNA 
The orientation of the terminated sloping or horizontal 
rhombic antenna is depicted on the next page in Figure 13.1, 
where 1 is the length of the antenna wires (numbered 1 to 4), 
Hy is the feed height above ground, H, is the height of the 
termination point above ground, y is the angle subtended by 
the feed wires (1 and 2), « is the angle between the x axis 
and the projection of the feed wires in the x-y plane, o is 
the angle between the plane which contains the antenna and the 
x-y plane, R is the radial coordinate, ÓÜ is the elevation 

coordinate, and $ is the azimuth coordinate. 
The Mathcad application for the rhombic antenna requires 


poc roll ewe s. 


B Nus length of the antenna legs (meters) 

I POS feed height above ground (meters) 

Hoc eMe height of termination point above ground (meters) 
a ecc m angle subtended by the antenna legs (degrees) 

b A. me Operational frequency (Hertz) 

RE e oes distance from antenna (meters) 

eye TN relative dielectric constant of ground plane 
DU nO conductivity of ground plane 

A elevation angle index (from Table 3.1) 


The feed height (Hj and termination height (H; can be any 
value greater than zero. The angle o is found from H,, H,, 
wire length (1), and the half-angle between the wires (y/2) by 


H-H, 


E CTS" 
21:c08| Y) 
2 


a=sin”? 
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Z axis 


FIELD POINT 
(R,0, 0) 
9 


m ais 
“Terminated Rhombic Antenna Orientation 


FIGURE 13.1: Spatial orientation of the terminated sloping 
or horizontal rhombic antenna for its Mathcad application. 





The angle y is then found from a and 2712 9:7 


tan Y) 
"n-2tan !| ——-—— 


SOS m 2) 





The user inputs the frequency (f,) for which the radiation 
parameters discussed in Chapter 3 are computed. Distance from 
the antenna (R) must meet the far-field requirements of 


Chapter 3 C. The elevation angle index (d) sets the @ 
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Sonda rome wilehe a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285? to 88.857? in 
increments of about 2?. Indices between those listed can be 
used to interpolate a better estimate of a desired elevation. 

The radiated electric field for the rhombic antenna is 
obtained in a manner analogous to that for the vertical dipole 
discussed in Chapter 4. The radiated field of each wire is 
considered separately, and superposition is used to determine 
Ehemtocal electric field. The rhombic's electric field 


equations are [Ref 6: p. 183-191] 
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ee 
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where U,, U,, U,, and U, are given by 


U,=1-cos(Y,) U,=1-cos(¥,) 
EOS a) A E 
where 
cos (V.) cos (0) sin(«) *sin(0)cos(a)cos(O-n) 
cos (Y,) cos (0) sin(«) +sin(0) cos (a) cos ($-n) 
cos (¥,) =sin(Q9) cos (a) cos (@-1n) -cos (8) sin (a) 


cos (¥,) =sin (89) cos (@) cos (@+n) -cos (6) sin (a) 


The rhombic’s radiated electric field is a varying combination 


of 6 and $ components, and its radiation pattern varies with 


changes in $. The total radiated electric field is the vector 


sum of the ÓÜ and $ components. Equations 13.3 and 13.4 are 


the $6 (vertically polarized) and ¢ (horizontally polarized) 


89 


components, respectively for wire one, and equations 13.5 and 
13.6 are the 6 and $ components, respectively for wire number 
two. Equations 13.7 and 13.8 are the analogous equations for 
wire three, and 13.9 and 13.10 are those for wire four. The 
first two terms inside each set of brackets represent are the 
space wave components, and any additional terms represent the 
surface wave components. In, is taken to be one, since a 
Sinusoidal current input with a maximum of unity is assumed. 
The requested inputs are used to calculate the following 
variables using a constant q of $-0 and $-7/2 for 632 discrete 
values of 0 which are equally incremented from -7/2 to m/2: 
Nese Gee wavelength of the operational frequency (meters) 
A free space wavenumber for operational frequency 
o cm index of refraction 
complex numerical distance (vertical polarization) 
complex numerical distance (horizontal polarization) 
vertical reflection coefficient 
E horizontal reflection coefficient 
vertical surface wave attenuation factor 
horizontal surface wave attenuation factor 
The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 for $20 and ¿=1/2. The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding § and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of @ to obtain the total radiated electric field 


distribution for the $20 and ¿=1/2 vertical planes. The space 


wave, surface wave, and total electric field results are then 
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normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete @ to depict the radiation patterns. The rhombic 
antenna Mathcad application computes the space wave, surface 
wave, and total electric field radiation patterns and 
radiation parameters in the $-0 and $-7/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
selected elevation angle (0, as $ varies from O to 27 in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the @ and tecla ena 
components to the total average radiated power, Pa, equations 
3.8 and 3.9 are used to integrate equations 13.3 to 13.10 over 
the hemispherical Gaussian surface above the ground plane at 
a fixed radius (R) from the antenna. The sum of the integrals 
is the antenna's P,,. With the discrete values of the electric 
field and total average radiated power determined, the Mathcad 
application predicts the following radiation parameters using 


the equations from Chapter 3: 


PME radiation resistance (Ohms) 

I a aci directivity 

E RP TENA effective isotropic radiated power (Watts) 
A E maximum theoretical effective area (square meters) 
E maximum theoretical effective length (meters) 

e Doti numerical distance (vertical polarization, 0290?) 
p INE numerical distance (horizontal polarization, 0-290?) 
Angle, .elevation angle of maximum directive gain (degrees) 


ON 


The directivity (Dj, effective isotropic radiated power 
(EIRP), maximum theoretical effective area (Am), maximum 
theoretical effective length (laa), and the elevation angle of 
maximum directive gain (Angle,,) are all determined for both 


the $20 and $sm/2 vertical planes. 
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XIV. THE TERMINATED SLOPING DOUBLE RHOMBOID ANTENNA 
The orientation of the terminated sloping double rhomboid 
antenna is depicted on the next page in Figure 14.1, where 1, 
is the length of the four short antenna wires (A, D, E, and 
H), 1, is the length of the four long antenna wires (B, C, F, 
and G), 1, is the length of the radials which connect the feed 
to the two termination points, Hj, is the feed height above 
ground, H, is the height of the termination points above 
ground, y is the angle between by the radials which connect 
the feed to the termination points, 6 is the angle between the 
termination radials and the short feed wires, y is the angle 
between the termination radials and the long feed wires, qa is 
the angle between the plane which contains the antenna and the 
x-y plane, R is the radial coordinate, @ is the elevation 
coordinate, and $ is the azimuth coordinate. 
The Mathcad application for the double rhomboid antenna 


requires the following inputs: 


"TT length of the short antenna wires (meters) 

l length of the long antenna wires (meters) 
Se a length of the termination radials (meters) 
ia feed height above ground (meters) 

A aE height of termination points above ground (meters) 
E E angle subtended by termination radials (degrees) 
Ee eee ne operational frequency (Hertz) 

A distance from antenna (meters) 

Es ee relative dielectric constant of ground plane 

Ur MeL conductivity of ground plane 

d. WERE elevation angle index (from Table 3.1) 


The feed height (Hj and termination height (Hj) can be any 


value greater than zero. The length of the radials to the 
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Sloping Double Rhomboid Orientation 


FIGURE 14.1: Spatial orientation of the terminated sloping 
double rhomboid antenna for its Mathcad application. 
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termination points must be less than the combined lengths of 
a long and short antenna wire (1, < 1,+1,). The angle y is 


found from the law of sines and cosines by 


1 
ST 
i, 


cele 


T) 
-1 1 3 2 4 ot 
cos eS MEAN Lo) 





The angle ð is also found from the law of cosines by 


(1,)2+(1,)?- (15)? 


21,1, M (142) 


d=\cos"+ 








The angle o is then found from Hy, H;, 1,, 1;, and the half- 
angle between the wires (y/2) by 


H, -H 


-— ——À (pras) 
1,cos|3 8] «1, cos(T «n 


ASI 
The user inputs the frequency (f, for which the radiation 
parameters discussed in Chapter 3 are computed. Distance from 
the antenna (R) must meet the far-field requirements of 
Chapter 3.C, and the elevation angle index (d) sets the @ 
coordinate for which a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857? in 
increments of about 2°. Indices between those listed can be 
used to interpolate a better estimate of a desired elevation. 

The radiated electric field for the double rhomboid is 
obtained in a manner analogous to that used for the vertical 
dipole discussed in Chapter 4. The electric field equations 


for the double rhomboid are given by [Ref 6: p. 238-244] 
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where U,, U, U,, Uy, Us, Us Un ana tarea” 
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The double rhomboid's radiated electric field is a combination 
of 0 and 4 components, and its radiation pattern varies with 
changes in $. The total electric field is the vector sum of 
the 6 and $ components. The radiated fields of each antenna 
wire are considered separately, and superposition is used to 
determine the total electric field. Equations 14.4 and 14.5 
are the 0 (Vertically polarized = ancdma Cnn (orem mE NE. 
polarized) components, respectively for wire A, and equations 
14.6 and 14.7 are the Y and $ components, respectively for 
wire B. Equations 14.8 and 14.9 are the analogous equations 
for wire C, 14.10 and 14.11 those for wire D, 14.12 and 14.13 


those for wire E, 14.14 and 14.15 those for wire F, 14.16 and 


14.17 those for wire G, and 14.18 and 14.19 those for wire H. 
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The first two terms inside each set of brackets represent the 
Space wave components and any additional terms represent the 
surface wave components. Ij is taken to be one since a 
sinusoidal current input with a maximum of unity is assumed. 
The requested inputs are used to calculate the following 
variables using a constant q of $-0 and $-7/2 for 632 discrete 
values of 0 which are equally incremented from -m/2 to m/2: 


D S te ae wavelength of the operational frequency (meters) 
Hp free space wavenumber for operational frequency 
A index of refraction 

i complex numerical distance (vertical polarization) 
complex numerical distance (horizontal polarization) 
E NM vertical reflection coefficient 

Mu a horizontal reflection coefficient 

vertical surface wave attenuation factor 
horizontal surface wave attenuation factor 


o 


< 


url UC) 


The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 with $20 and $-7/2. The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding Ó and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of to obtain the total radiated electric field 
distribution for the $20 and $-7/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete @ to depict the radiation patterns. The Mathcad 
double rhomboid antenna application computes the space wave, 
surface wave, and total electric field radiation patterns and 
radiation parameters in the $-0 and ¿=1/2 vertical planes. 
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The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0, as $ varies from O to 2m in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the 6 and $ electric field 
components to the total average radiated power, Paa, equations 
3.8 and 3.9 are used to integrate equations 14.4 to 14.19 over 
the hemispherical Gaussian surface above the ground plane at 
a fixed radius (R) from the antenna. The sum of the integrals 
is the antenna's P4. With the discrete values of the electric 
field and total average radiated power determined, the Mathcad 
application predicts the following radiation parameters using 


the equations from Chapter 3: 


ROME LIT radiation resistance (Ohms) 

jo, c 0t RE directivity 

BIRD 2 ee effective isotropic radiated power (Watts) 

a maximum theoretical effective area (square meters) 

UNT maximum theoretical effective length (meters) 

p Mm numerical distance (vertical polarization, 0-90?) 

Boo numerical distance (horizontal polarization, 0-90?) 

Angle, .elevation angle of maximum directive gain (degrees) 
The directivity (Dj, effective isotropic radiated power 
(EIRP), maximum theoretical effective area (Ama), maximum 


theoretical effective length (lax), and the elevation angle of 
maximum directive gain (Angle,,) are all determined for both 


the $20 and $-7/2 vertical planes. 
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XV. THE VERTICALLY POLARIZED LOG-PERIODIC DIPOLE ARRAY 
The orientation of the vertically polarized log-periodic 
dipole array is depicted in Figure 15.1, where N is the number 


of dipoles, n is an index equal to 0,1,2..N-1, 1, is the 


length of the n?«1 dipole, d, is the distance between the n*?^«1 
and n"«2 dipole elements, H, is the height of the antenna feed 
above ground, H, is the height of the n?«1 dipole center, Y is 
the angle between the vertical axis and the array's center 
axis, o, is the angle between the center axis and the line 
which connects the upper tips of the elements, o, is the angle 
between the center axis and the line which connects the lower 
tips of the elements, o, is the angle between the horizontal 
and the line connecting the lower tips of the elements, R is 
the radial coordinate, @ is the elevation coordinate, and ¢ is 
the azimuth coordinate. 


The Mathcad application for the vertically polarized log- 


periodic array requires the following inputs: 


IDE oa number of dipole elements 

DT TER length of the shortest dipole (meters) 

A length of the second shortest dipole (meters) 

c T M distance between first two elements (meters) 
rad ot. radius of the shortest dipole (meters) 

Ho cM, C feed height above ground plane (meters) 
MU ees angle from vertical axis to array axis (degrees) 
BE o operational frequency (Hertz) 

RS LL eus distance from array (meters) 

NEUE c. relative dielectric constant of ground plane 

eue d esci conductivity of ground plane 

WINE os elevation angle index (from Table 3.2) 

ADMa transmission line characteristic admittance 
INMB. termination impedance connected to longest dipole 
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FIGURE 15.1: Spatial orientation of the vertically polarized 
log-periodic dipole array for its corresponding Mathcad 
application. 
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The user inputs the frequency for which the radiation 
parameters discussed in Chapter 3 are computed. The lengths 
of the first two elements (1, and lj) are used to calculate the 


Tega. tedhesrelarcionshipgeia given by [Ref 10: p. 317] 


fo (n=0,1,2....N-2,N-1) MESSIS 





The remaining element lengths (l, to ly,) are then found from 


nel 


D 
1.,,== (duo) 
t 


As noted in equation 15.1, indices start at zero instead of 1, 
So an index simply refers to the index-plus-one position in 
terms of successive values for a given parameter. This is to 
maintain continuity with Mathcad which also indexes from zero. 
The radius of the first element (rad) is used to find the 
remaining radii by substituting rad, and rad,,, for 1, and l,,, 
in equation 15.2. The angle between the vertical axis and the 
array axis (Y) can be between 0? and 90°. The angles a, and 


a, are found by substituting 


1,4, 5 - Y (15.3) 
11260 


d z 
; 1-7 (15.4) 
Sia Sa.) Tar e.) Cosa + 040) 





which yields 


h,+d,sin{ > - m 


&,-tan'! 
e cos(5 -*| 
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The angle a, is then found from [Ref 10: pp. 351-354] 


Sin(a,) elie 


-— (155 859 
costa a ua COS 


Or 


Sin cascos ora) 
Cos (@,)'*sin(@,)cin (ance) 


&,-tan! 

The remaining separations (d, to d,,) are then found from 
equation 15.4. The dipole half-lengths (h,) are found from 
1,/2, the feed height (Hy) must be at least equal to the half- 
length of the longest dipole (h,,), and distance from the 
array (R) must meet the far-field requirements of Chapter 3.C. 
The elevation angle index (w) sets the @ coordinate for which 
a horizontal radiation pattern is determined. Table 3.2 lists 
possible indices and their corresponding elevation angles from 
.579 to 89.1? in increments of about 2.3 Indices dew 
those listed can be used to interpolate a better approximation 
of a desired elevation. The termination impedance (TIMP) is 
assumed to be connected to the array's center axis opposite 
the feed at a distance of h,,/2 from the longest element. 
Typical termination impedances are the transmission line’s 
matched impedance (complex conjugate) or a short circuit 
(zero). The characteristic admittance (ADM) is simply the 
inverse of the transmission line’s characteristic impedance. 
Tf computed lengths, separations and radii are not 
desired, they can be entered manually as described on the page 
two of the application. Consult a Mathcad manual to be 


certain the entries are made correctly. 
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The open Circuit impedance matrices must be calculated to 
determine the current distribution among the dipole elements. 
The matrix of open circuit mutual impedances between the 


dipole elements is represented as 


: 
Zoo Lo,1 Lo,2 EN 5 Zo, (N-1) 
1,0 211 21,2 Z1, (N-1) 
[Z] = 22,0 La, 1 La 2 eww t Z5, (N-1) (15 e 
|S (w-1),0 Z (N-1),1 Z (N-1),2 v T E Z (N-1) , (N-1) 


Each main diagonal term (2,) is the self-impedance of the i"«1 
dipole. The off-diagonal terms (Z) are the mutual impedances 
between the i%+1 and k®*+1 elements. Subscripts i and k are 
matrix indices. from 0,1,2... .N-1. 

The matrix of open circuit mutual impedances between the 


actual dipoles and their images is represented as 


a ZI A LlÍo 2 CELL Zl(w-1 
¿Lio Zla ¿liz o s Zl, (N-1) 
IZT] 3 ¿loo ZL, 1) ¿la 2 AN ZI; (-1) (15.7) 
Zlino Lleno Pera + 5 5 + SL (ery, (ner) 
The main diagonal terms (ZI) are the mutual impedances 


between the iW+1 element and its image, and the off-diagonal 
terms (ZI,,) are the mutual impedances between the i%+1 element 


and the image of the k?«1 element. 


Bo 


The main diagonal self-impedance terms of equation 15.6 


(Zu)lsare toundotrsomal bere p Oc 
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213° Sos (Phe | (Oe) -20(U,)] (15.8) 


Lel (0(V,) -20(V,) ]+2[0(U’,) -O(U,) -O(V,) ]+20(U',)[1+c0s (282) ] | 
y 
u-B[/2radi«ani-2n] ^ v^-/28rad, Y= [2radj+4ahj+2h,] 
U= | /2zad?+h?-h,| v, =P | /2rad?+h?+h,| 


For W=90°, the off-diagonal mutual impedance terms of 


where Q(x) =Ci(x) -jSi(x) =f ges ua ay jp dy Moon 
co Q 


equation 15-6 (2) aac son Vou eN 


mr, J ee O eee 
Lk "Ei cos[B(h;-h)]-cos[B (h;*h,) ] 


Z 
[e 7* ^? t0 (u,) -o(u) -o(0,) |ee" "n^" to(v,) -0(V) -9(,] 
«e TP Urtpp (gl -o(V,) -Q(U,) ee?" ^P to( v^.) -o(0,) -0(V;] 


*20(BDj ¿)(cos [P(h,-Aj)]+cos[P(A,+h,)])] 


where Q(x) =Ci(x) -j-Si(x) =f £95) ay - ¡EE gy 


y 0 
U,-p D (h*By)*- (hy+h,) | Voz | pes He (btp) | 
Up (DE + (2,2)? (2-5) | vs p VD} kt (hi-ho? (hy-h,) | 
p 


U,-p JD? ,«ni-hl v,=B |/D? .+h? +h, | 
U,=6 VD ,«hi-hj vV;-p | D 
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For Y less than 90°, the dipoles are in echelon. The 
expression for the mutual impedance terms of equation 15.6 


(2,,) for an echelon configuration is given by [Ref 12] 


Z, ,-15[e p OU -O(U,) pe Poseo (y) -Q(Vi)] mp5 SH 


te TPES? OCU) ZO pe Ec) ON Va) | 


+ eae poc OK U,) TO U,) +e? S ( V;) S OA | 


Il 


ee S 79, [9 (0, ) -O(U, pedis 1*2 hy* SD, & "Io (V) -O(V,)) 
+2cos(Bh,) e 7*5 EO Cw.) -Q(W,)]*2cos (Bh) eno (Y) -OCY,)] 


+2c08 a AO -O(W,) AO) ET XE | } 


where Q(x) =Ci(x)-j:(Si(x) = SOSY) ay-j-f RY ay 
0 


8 —— x 


U,-p D? ,« (SD; ,- hj) ?« (SD, ,- hj) | 
V;-B l/p x» (Sb; ,-h?- Gb; ,- hj | 
psp yD Dj k* (SD; ,* hj) ?- (SD; ,* hj J 
V/,- p [/D£ ,* (SD,, hy) 2+ (SD, ,+h,) | 
ND .+ (SD, ,-hy+h,) 2+ (SD, ,-hy+h,) 
[DE ¿+ (SD, ¿A+ hy) > (SD; ,-Ayt+h,) 


V,-B |/D2 ,« (SD, ,-h;*2h ?- (SD; ,-h;*2hy) | 


U,-p D at (SD; yc h;*2h,) ^- (SD; ,+h,+2h,| 
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V,=6 [yD k+ (SD, ti), +21) oe eee 
W,=B |/D? ,+SD? .-SD,, y, 
Y, =P DE +SDÍ g+SD,, y) 


W,=B |/D? e+ (SD; the) (SD; thy) | 


Y,-p p e (SD rh) 4 (SD; .+h,) | 
W,-p | Dj k* (SD ,*2h,)?- (SD, ,*2hj) | 
"n | (DE ,- (SD, ,*2h,) ?* (SD; ,*2h) | 
The D, terms in equations 15.10 and 15.11 are the perpen- 


dicular distances between the i%+1 and k"+1 elements given by 


k-1 
D, 5Dx, => dy: sin(®) GeT) 


msg 


The SD, terms fOr equate mS rece me eee? c 
from the center of the i"«1 dipole to the z coordinate of the 


nearest end of the k"+1 dipole and are given by 


k-1 
SD; ,5SD, ;2 9 discos yh; 


J 
n=l 


iS 59) 


For ground planes with high conductivities relative to 
the operating frequency, it is necessary to compute equation 
15.7, [21], for the mutual impedances between the elements and 
their images. The main diagonal terms (ZI) are the mutual 
impedances between the i%+1 dipole and its own image. These 
Z; terms are calculated from the mutual impedance equation for 


a collinear configuration given by Re cmi 


dub 


zr, nisle Ep (Aj Eis A" [o (U,) O CU ens 1n 





f CHI. E E Mq Bi. *2h, 
+e tb (hy = Oa BOL) Fe Jp (h; E zu i 
sso rol eo quee a in 
: AMET 2j ——À Ao h. 
A O e M n ra : 


+2c0s (Ph) e PO (Y,) -0(Y,)]+2c0s (Ph) e Es EET = | 


| 


Mieres OS CA SL) = [LED ay 7-f SAB) ay 
oo Q 


-2cos (Bh) Sg p 


#210(¥,) -O(Y) ] 
Mai. 
DINE EE 


m2 


*2cos (Bh;) e 7P Chithiu ^u 





U,- p [/D? c» (hi; A)? (ai, hp] 
vB |D} (hiz A- (his, erhi) | 
U' =P | YD? e+ (hii gth;)?- (hi; gth) | 
E +A) | 
U,-B |/D? + (Bi, y+, ?+ (hi; ,-Ay+h,) | 
v,-B [/D? thi, -hy+h,) 2- (hi; ,-y+h,) | 
U,-p DE &* (hi; jh; hy) ?- (Bi, ¿hh | 
v,-B YD? p+ (hi, ,+hj+h,) (m1, heh | 


U,-p | D a) (hi, ,-h,*2hy | 


V,-B |/D£ r*+ (hi, ,-h;*2hj) ?- (ii; ,-h;*2hj) | 


Ter: 


U,-p | Di 2 (hi, ,«h;*2hy | 
VID | Di, kt (hi, ,+hy+2h,)*+ (hi; ,+hy+2h,) | 


WB |D erhi e Bi  Y,=6 |/D? hi? +hi, , 
w,=B |yD? .+ (hi, ,+h,)?- (hi; +h) | 
Y,-p | Dj, «+ (hi; ,+h,)?+ (hi; «hy | 
W,-p |/D2 ,- (hi, ,*2h) ?- (hi; e 2h49 | 
Y,-p y DAS (hi, ,+2h,) | 
The symbol hi, in equation 15.14 is the height of the i%+1 


Gipole above the nearest end of its own image, and with i=k is 


obtained from 


The off-diagonal mutual impedances between the dipoles and the 
images (ZI,,) are found using equation 15.13, except the Q(x) 
arguments are those giveniifor eguacioni MIA QE UU 
replaced by hi, E&romweguabsong Ts 

The function Q(x) given by equation 15.9 is defined in 
terms of sine [Si(x)] and cosine [Ci(x)] integrals. Mathcad 
is incapable of evaluating Si(x) and Ci(x) directly, so a 
series expansion and polynomial approximation are written into 
the Mathcad code for both the sine and cosine integrals to 
evaluate the functions Q(x). The arguments for Q(x) are 
always real, so it is unnecessary to find Si(x) and Ci(x) 
expressions valid for complex numbers. For arguments less 
than one, Si(x) and Ci(x) are evaluated by series expansions 
given by [Ref 9: p. 232] 


T12 


= -1)7? xen+1 . 
E AE (2n+1)! E 


n=0 


n 2n 
O TE xy VC EE 7s) 


where y=0.5772156649 


For arguments greater than one, Si(x) and Ci(x) are 
defined in terms of two auxiliary functions, f(x) and g(x) 


given by [Ref 9: p. 233] 


pesi x9*ra x*ca,x**a4x^«a, E a TEN 
PE EDU acd 
where jenes» «uEitio-" 
a,738.027264 b,=40.021433 
a 206 187023 b,=322.624911 
a,=335.677320 ae O25 62 S10 
a,=38.102495 do, = 157 895423 
8+ + 44 24 
gun. B «e (x) QNS, 
where ¡MES LO 
a,=42.242855 b,=48.196927 
a,=302.757865 , . b,=482.485984 
= e AOT Jud 979885 
a,=21.821899 b,=449 .690326 
The sine and cosine integral functions are then defined in 
terms of these auxiliary functions by [Ref 9: p. 232] 
Si (x) =4-£(x) cos (x) -g(x) sin (x) MISION 
o) 


Ci(x)sf(x)ysan(x)-g(x)cos(x) 
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In addition to the mutual impedance matrices, it is 
necessary to calculate the short circuit admittance matrix to 
determine the elemental current distribution. The short 


circuit admittance matrix is Given by cs PONDUS 


VO 0 
Yo YQ Y 0 0 0 
O Y, E Ya Y, 3 O 
[Y] =| 0 O ae a (519751 
0 0 . -— Y(y2),(8-3). Yta-2), (N-2). Y (8-2) , (I1) 
0 Y ay taz) Ya) U-1 | 


where the diagonal terms are given by 
Y,o--JY.cot(Bd) and Y, rm "Y 7 JY,Cot (Bdy.,) 
and forQ<n< N-1: 
Y uns ES Des SII 


and the off diagonal terms are given by 


A y =-JY¿csc(Pad,..,) 


n, (n-1 


The term Y, is the characteristic admittance of the trans- 


mission line given by the input (ADM), and Y’; is given by 








cos| p > 
Y' Y, 


v,zscos[B Pa), jsin(B VT 


where Z, 1s the termination impedance input (TIMP). 








The ground conductivity relative to the operating 


frequency (Ca) is then computed by [Ref 6: p. 640] 


a) allo. Teen 


rel 
WEE, f. 5 


(15.24) 
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to determine which equation to use for the array’s current 
distribution. A surface with a relative conductivity greater 
than 20 1S considered to be highly conductive relative to the 
transmitted frequency and the current distribution is given by 


FRED 599427 


ME LtUIREF| (15.25) 


where [U] is the N X N identity matrix. The entries of the 
[IB] matrix are the base currents of the dipole elements. The 
matrix [IREF] represents the input current to which the [IB] 
distribution is referenced. Since the current input is only 
to the feed at the center of the first element, and since the 
input current is sinusoidal with a maximum of unity, [IREF] is 
alXN matrix with a one as the first entry and N-1 zeroes in 
the remaining positions. 

For surfaces with a relative conductivity of less than 
20, the mutual impedances due to the images are neglected and 


the current distribution is given by 


[IB] -([U] «[Y] LZ] y! [IREF] (159559) 


The electric field equations are referenced to the 
current maxima values for the dipole elements, [I], which are 


obtained from the [IB] matrix using [Ref 6: p. 209] 


pa ll (T5 27) 
i sin (BA,) | 


TES 


The mutual impedance calculations are not valid when any 
element length is an exact integer multiple of the wavelength. 
When this occurs, there is a Singularity error in the mutual 
impedance calculations, and it is necessary to vary the 
frequency such that no element is an exact integer multiple of 
the wavelength. The change required is only one or two 
percent, and the predicted radiation parameters are still a 
good estimate of those for the original frequency. 

The electric field for the vertical log-periodic array is 
obtained for each individual dipole in a manner analogous to 
that used for the vertical dipole discussed in Chapter 4. The 
equations for the individual elements are combined into a 
Single expression for the array’s total electric field by the 
array factor. The array factor for the vertical log-periodic 


array is [Ref 6 p. 222] 


N-1 
A+jB=Y Fete cos [Bh,cos (8) ]-cos(BA;)] - CON 


1=0 


A 4 (1-T.) F, [sin (0) " p orn (0) COS (8) )] 
n 





eJBY;sec («42*44) [cos (00 sin (a2*a,) *sin(8)cos (a2*a,) sin($)] 


where Y, is the y coordinate of the i,+1 dipole given by 
Ic 
Y,=| Y d,[sin(Y) (15-29) 
n=0 
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The first two terms inside the brackets of equation 15.28 
represent the space wave, the third term represents the 
surface wave, and I, is the current term for each element based 
on a sinusoidal current input with a maximum of unity. The 
array factor is combined with the element factor to find the 
expression for the total radiated electric field given by 


-JBR 7 
E,- 7J60e A+ JB 


R sin (8) cil 


The requested inputs are used to calculate the following 
variables using a constant $ of $20 and $-7/2 for 312 discrete 


values of @ which are equally incremented from -m/2 to m/2: 


NU vs half-length of the i%+1 dipole (meters) 

AA i height of the i^"«1 dipole (meters) 

pe y coordinate of the i"«.1 dipole 

— current for the i"«1 dipole 

D m wavelength of the operational frequency (meters) 

e free space wavenumber for operational frequency 
IT index of refraction 

CCS complex numerical distance for vertical polarization 
OA vertical reflection coefficient 

F vertical surface wave attenuation factor 


The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 with $20 and $-7/2. The space wave and surface wave results 
are combined for corresponding values of @ to obtain the total 
radiated electric field distribution for the ¢=0 and q=1/2 
vertical planes. The space wave, surface wave, and total 
electric field results are then normalized with respect to the 
maximum field intensity of each, and the normalized magnitudes 
are plotted for each discrete 0 to depict the radiation 


patterns. The Mathcad vertical log-periodic dipole array 


im 


application computes the space wave, surface wave, and total 
electric field radiation patterns and radiation parameters in 
the $20 and $-7/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (w) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0,) as $ varies from 0 to 21 in 312 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

Equations 3.8 and 3.9 are used to integrate equation 
15.30 over the hemispherical Gaussian surface above the ground 
plane at a fixed radius (R) from the array to find total 
average radiated power (Paa). With the discrete values of the 
electric field and total average radiated power determined, 
the Mathcad application predicts the following radiation 


characteristics from the equations in Chapter 3. 


rn radiation resistance (Ohms) 

Do. 08. directivity 

EBENE. . ie effective radiated isotropic power (Watts) 

Ae e maximum theoretical effective area (square meters) 

ee eee maximum theoretical effective length (meters) 

De e numerical distance (vertical polarization, 602909?) 

BI. ee numerical distance (horizontal polarization, 0-90?) 

Angle, .elevation angle of maximum directive gain (degrees) 
The directivity (D), effective isotropic radiated power 
(EIRP), maximum theoretical effective area (A,,,), maximum 


theoretical effective length (1,,), and the elevation angle of 
maximum directive gain (Angle,,) are all determined for both 


the $20 and $-7/2 vertical planes. 


PES 


As an example, the Mathcad vertical log-periodic array 


application was executed with the following inputs: 


number of elements 12 
lengenson first dipole 3.246 meters 
length of second dipole 2.732 meters 
first separation distance 1.096 meters 
radius of first element 0.00325 meters 
height of antenna feed 8.0 meters 
frequency 10-10? Hertz 
distance from the antenna 3000 meters 
relative dielectric constant 4 

ground conductivity fe ena 
characteristic admittance 1/450 Mhos 
termination impedance O Ohms 
angle from vertical to array axis 90° 
elevation angle index 285 (47°) 


Figures 15.2 through 15.7 are the space wave and surface wave 


radiation patterns in the $-0 and $-7/2 vertical planes, and 


che designated horizontal plane for this exampie. 
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FIGURE Vertical log- 
periodic space wave radiation 
pattern for d$-m/2 vertical 
plane. 
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periodic space wave radiation 
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FIGURE 15.5: Vertical log- 
periodic surface wave 
radiation pattern for q=1/2 
vertical plane. 
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FIGURE 15.7: 


The following radiation parameters were predicted by the 


Mathcad application for a sinusoidal current input of one Amp. 


Total power radiated (Watts) nu go 
Radiation resistance (Ohms) 75.934 
Numerical distance (vertical) ENG 
pirectivity SES 12. 196 
EIRP (Watts) 129.656 462.272 
Max eff area (sq meters) 244.578 Sze One 2 
Max eff length (meters) 14.038 AS DE 
Angle, (degrees) Dal 21.14 


These results are consistent with expectations for this 
pare cular cont 1guration. Appendix D contains computer 
hardcopies of additional example calculations for the vertical 
log-periodic dipole array and compares predicted radiation 


parameters to those expected based on previous calculations. 
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XVI. THE HORIZONTALLY POLARIZED LOG-PERIODIC DIPOLE ARRAY 





The Orientation of the dhorsvzonpod EC HEN 
periodic dipole array is depicted in Figure 16.1, where N is 
the number o£ dipoles, n is an index equal ron 1L 
is the length of the n"«1 dipole, d, is the distance between 
the n"«.1 and n"«2 dipole elements, H, is the height of the 
antenna feed above ground, H, is the height of the n?«1 dipole, 
Y is the angle between the vertical axis and the array's 
center axis, o is the angle between the array's center axis 
and the lines which connect the tips of the elements, R is the 
radial coordinate, Ü is the elevation coordinate, and $ is the 
azimuth coordinate. 


The Mathcad application for the horizontally polarized 


log-periodic array requires the following inputs: 


Nec cM number of dipole elements 

TUM A length of the shortest dipole (meters) 

Ls e n length of the second shortest dipole (meters) 

Go sfr distance between first two elements (meters) 
rol cm radius of the shortest dipole (meters) 

Ho ELS feed height above ground plane (meters) 

pre ce angle from vertical axis to array axis (degrees) 
[oe c Operational frequency (Hertz) 

Ries cee distance from array (meters) 

p EU relative dielectric constant of ground plane 

Qu co RE conductivity of ground plane 

Wi E elevation angle index (from Table 3.2) 

ADM... ae transmission line characteristic admittance 

T IMP... eae termination impedance connected to longest dipole 


The user inputs the frequency for which the radiation 


parameters discussed in Chapter 3 are computed. The lengths 
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FIGURE 16.1: Spatial orientation of the horizontally 
polarized log-periodic dipole array for its corresponding 
Mathcad application. 
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of the first two elements (l, and 1,) are used to calculate the 


log-periodic relationship (7) given Dyke cR ID s 





T= —— (n=0,1,2.. "NS (5.1) 


The remaining element lengths (1, to l,,) are then found from 


E 


TIS 


1 
= (16.2) 
= 


As noted in equation 16.1, indices start at zero instead of 1, 
so the index simply refers to the index-plus-one position in 
terms of successive values for a given parameter. This is to 
maintain continuity with Mathcad which also indexes from zero. 
The dipole half-lengths (h,) are found from 1,/2. The radius 
of the first element (rad) is used to find remaining radii by 
substituting rad, and rad,,, for 1l, and l,,, into equation 16.2. 
The angle between the vertical axis and array axis (Y) can be 


between 0° and 90°, and œ is found from [Ref 10: p. 319] 


Sa = (1-1) cot (a) (16.3) 


n+1 





which can be written as 


(iS DNE. 
d 


a=tan”?* 
n 


| (16.4) 


The remaining separations (d, to dy) are found from equation 
16.3. Feed height (Hj) can be greater than or equal to zero. 
Distance from the array (R) must meet the far-field require- 
ments of Chapter 3.C. The elevation angle index (w) sets the 


0 coordinate for which a horizontal radiation pattern is 
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determined. Table 3.2 lists possible indices and their cor- 
responding elevation angles from .57° to 89.1° in increments 
Orta 2.37; Indices between those listed can be used to 
interpolate a better estimate of a desired elevation. The 
termination impedance (TIMP) is assumed to be connected to the 
array's center axis opposite the feed at a distance of h,,/2 
from the last element. Typical termination impedances are the 
transmission line's matched impedance (complex conjugate) or 
a short circuit (zero). The characteristic admittance (ADM) 
is just the inverse of the transmission line's characteristic 
impedance. If computed lengths, separations and radii are not 
desired, they can be entered manually on page two of the 
application. Consult a Mathcad manual to be certain the 
entries are made correctly. 

The open circuit impedance matrices are calculated to 
determine the current distribution among the dipole elements. 
The matrix of mutual impedances between dipole elements is 


represented as 


Zo, o Zo,1 Gora + + + + £o, (NA) 
210 241 41,2 mr l, (N-1) 
[Z] = ; , . m i (ES!) 
Liar o Lao Lar 2 * * c * naa), (0-10) 

Each main diagonal term (2Z,) is the self- impedance of the 
i^-«1 dipole. The off-diagonal terms (Z,) are the mutual 
impedances between the i%+1 and k"+1 elements. Subscripts i 
dgugEIEamesmabrIx indices (0,1,2....N-1). 
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The matrix of open circuit mutual impedances between the 


actual dipoles and their images is represented as 


Ze DH Z2. a M 
O arm ZI, 3 E E ; : Zl, (Ni) 
r] = i S OMNES (16.6) 
LL rayo Limo )iqeiya coc c + SL en), (1)! 


The main diagonal terms (Z21,) are the mutual impedances 
between the i%+1 element and its image, and the off-diagonal 
terms (ZI,,) are the mutual impedances between the i%+1 element 
and the image of the k"+1 element. 

The main diagonal self-impedance terms of equation 16.5 


are found from [Ref 6; pp- 205-206] 


2 60 -j2ph,; E 
^ii 1-cos2Dn») [e [9 (05) -20 (01) ] (16.7) 


PEAT [O V.) O OA a sos (2BA,)] | 


where 0(x)=Ci(x)-3j:Si(x) =f Sosy y) dy- zfs dy Giera 


Up | /2rad: con il Nu Vo=B | /2rad?+4h;+2h,| 
U,-p |/2rad?«ni-nj| v,-B | /2rad?+h?+h,| 


The off-diagonal mutual impedance terms of equation 16.5 


are given by [Ref 11] 


60 


DO cos [B (h;-h) ] -cos [B (4, +h,) ] 7 (16.9) 


[e 79 ^? to (9) -o(u) -9Q(0,) «e? P i^? t0(v.) -o0(v.) -O( Va) j 
e 7PUcP TO (U^. -9(V,) -Q(U,) ee? P ^77? to (v^) -o(0) -QCV2) | 
*20(BD; Q(cos(B(h;-h) ] «cosIB (A;*h,) 1)] 
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where Q(x) =Ci (x) -3+Si (x) - [ 9833 ay E jf RM ay 


=P [Diet (hith)? ith) vap [VD et (hyh) + (hith) 


u'=B VD} k*t (hiha ?- (hi-h,) | v! -p 
p 


The D., terms are the distances between the i%+1 and k%+1 


dipole elements and are given by 
k-1 


D; 47D, 17, d, (16.10) 
Fir 


For ground planes with high conductivities relative to 
the operating frequency, it is necessary to compute equation 
16.6, [ZI], for the mutual impedances between the elements and 
their images. The mutual impedance terms (Z1,,) are given by 


equation 16.9 with the Q(x) arguments 


=B Yaiza (hiha? (hith) voB [ydi e+ (hith) + (hith) | 
U',-B [/did.-(h;-h)?-(h-hg| — v^-p |/diZ ,- (h;- ho ?* (3;- 54) | 
U,-p [jai nia, v,=p [Jai nia, 
U,-p [/dif ,«h£-h,] v,-p [/ai£ .«nz«h, 


The di,, entries are the distance from the i%+1 dipole to the 


image of the k"+1 dipole and are given by 


di, ,=/ (H,+H,) *+ (D, ,8in(®) )° (16.11) 
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The function Q(x) given by equation 16.8 is defined in 
terms of sine [Si(x)] and cosine [Ci(x)] integrals. Mathcad 
is incapable of evaluating Si(x) and Ci(x) directly, so a 
Series expansion and polynomial approximation are written into 
the Mathcad code for both the sine and cosine integrals to 
evaluate the functions Q(x). The arguments for Q(x) are 


always real, so it is unnecessary to find Si(x) and Ci(x) 


expressions valid for complex numbers. 


For arguments greater than one, two auxiliary functions, 


f(x) and g(x) are determined by [Ref 9: p. 233] 


8 6 4 2 
imei ot x tduxX tad omite 
foy P (er (16.12) 
X\ 308 +) cS occus 


where lene 0-7 
a,=38.027264 b,=40.021433 
a,=265.187033 b,=322.624911 
a,=335.677320 b,=570.236280 
a,=38.102495 b,=157.105423 

x9*«a,x*«a,x'«a,x?«a 
g emr c dcc ci o (16.13) 
Xx*| xBUpxSIbooci pp GEM. 
where le EROS 
a,=42.242855 b,=48.196927 
a,=302./57865 b,-482.485984 
a,=352.018498 D:=11 M08895 
a,221.821899 b,=449 .690326 
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and Si(x) and Ci(x) are defined in terms of these auxiliary 


Eumeeions by WWRef 9:"p. 232] 


EDS s cos(x)-g(x)sin(x) 016014) 


Buc um )-cox)cos (x) uo s) 


For arguments less than one, Si(x) and Ci(x) are 


evaluated by series expansions given by [Ref 9: p. 232] 


fa) n x^n*l 


SS) m: (2n*1) 2n*1)! Ede 
E g (-1)2 x?n 
C1(x) =y+1n(x) + 2 Ea (EC) 


n 


where y=0.5772156649 


It is also necessary to calculate the short circuit 
admittance matrix to find the elemental current distribution. 


The admittance matrix is given by [Ref 10: p. 321] 


Y (N-2) 20N-3J Y (N-2) 7428-2] Y (N-2) , (N-1) | 


Co Co. 
© 


0 Yon-1),(N-2) Yin-1) (x-1) | 
where the diagonal terms are given by 
Saad Y 1 mem =Y 7 nct (Bady) 
and for 0 < n < N-1: 

e) J Yo Corpa mcos (pa) 


and the off diagonal terms are given by 


Y (0-1) n^ Y" 


nee oe Sel di 


129 


The term Y, is the transmission line's characteristic admit- 


tance obtained from the input (ADM), and Y’, is given by 


h h 
COS jn; Yr NN pm 


[e 
DA 


(16.19) 
h h 
v, zscos[a Pes J. sion = 


The ground conductivity relative to the operating 


frequency (Ca) 1s then computed by [Ref 6: p. 640] 


c _ 18000 o 


(e220) 
(coc TM 


O tel” 
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A surface with a relative conductivity greater than 20 is 
considered to be highly conductive relative to the frequency 
and the current distribution is given by (Ref 10: p. 342] 

[IB]-([U]« [Y] [Z] « [Y] [ZI] ' [IREF] (16.21) 

where (U] is the N X N identity matrix. The entries of the 
[IB] matrix are the base currents of the dipole elements. The 
matrix [IREF] represents the input current to which the [IB] 
distribution is referenced. Since the current is input only 
to the feed at the center of the first element, and since the 
input current is sinusoidal with a maximum of unity, [IREF] is 
a 1 X N matrix with a one as the first entry and N-1 zeroes in 
the remaining positions. 

For surfaces with a relative conductivity of less than 
20, the mutual impedances due to the image elements are 
neglected and the current distribution is given by 


[IB] ={ [U] + [Y] [Z] }' [IREF] (16.22) 


ins 


The electric field equations are referenced to the 
current maxima values for the dipole elements, [I], which are 


obtained rrom the [IB] matrix using. [Ref 6 p.209] 


(a (16.23) 
+ sin(ßh;) 

The mutual impedance calculations are not valid when any 
element length is an exact integer multiple of the wavelength. 
When this occurs, there is a singularity error in the mutual 
impedance calculations, and it is necessary to vary the 
frequency such that no element is exactly an integer multiple 
of the wavelength. The change required is only a percent or 
two, and the predicted radiation parameters are still a good 
estimate of those for the original frequency. 

The electric field for the horizontal log-periodic array 
is obtained for each individual dipole in a manner analogous 
to that for the vertical dipole discussed in Chapter 4. The 
equations for the individual elements are combined into a 
Single expression for the array’s total electric field by the 
array factors. The array factors for the horizontal log- 


periodic array are [Ref 6: p. 214] 


N-1 
Se= Y T¡¡cos[Bh,(cos (0) cos (¥) +sin (8) sin(P) sin (g) )]-cos (BA;) ] 


1=0 


ue i 4 i 2E d 
TS DI OSA e 2BAicos (0), 5i.n* (9). (16.24) 
n*cos(8) 





CU D (cos (0) cos (Y) -sin(0) sin(Y)sin($)] 


[sint (e) - 8-81 (GT cos (oy 
n? 





et, 


N-1 


SoD I ,{cos[Bh,(cos (8) cos (P) +sin(@) sin(¥) sin())]-cos (BA)! 


EX uc (1-D,) F e aaa (ic 25) 


ele Y,csc(V) [cos(0) cos (Y) -sin(0) sin( Y) sin($)] 


where Y, is the y coordinate of the i,+1 dipole given by 


rl asin? (EZ 6) 
n=0 

The first two terms inside the brackets of equations 
16.24 and 16.25 represent the space wave, the third term 
represents the surface wave, and I, is the current term for 
each element based on a sinusoidal current input with a 
maximum of unity. The array factors are combined with the 
element factors to find the total radiated electric field 
equation, and the expression for the radiated electric field 
distribution of the horizontal log-periodic dipole array is 


given by [Ref 6: p. 213] 


_-j60e7FF cos() cos (6) S 


E, : (Pc 
R 1-sin? (8) cos? (4) 

pe A (16.28) 
R 1-sin? (8) cos? (8) 


The requested inputs are used to calculate the following 
variables using a constant $ of $-0 and $-7/2 for 312 discrete 


values of 0 which are equally incremented from -7/2 to m/2: 


T» 


ee half-length of the i"+1 dipole (meters) 

Ho E. . height of the i"«1 dipole (meters) 

poe Ver. y coordinate of the i^"«1 dipole 

DAE E current for the i".1 dipole 

A C NN wavelength of the operational frequency (meters) 
DNE. free space wavenumber for operational frequency 
Tues index of refraction 

EN eos complex numerical distance (vertical polarization) 
P complex numerical distance (horizontal polarization) 
o vertical reflection coefficient 

[o c horizontal reflection coefficient 

Borm. vertical surface wave attenuation factor 

Pe este horizontal surface wave attenuation factor 


The calculated variables are used to evaluate the 
radiated far-field space wave and surface wave for the 
discrete values of @ with $¢=0 and ¢=7/2. The space wave and 
surface wave results are combined for corresponding values of 
0 to obtain the total radiated electric field distribution for 
the $20 and $-2-7/2 vertical planes. The space wave, surface 
wave, and total electric field results are then normalized 
with respect to the maximum field intensity of each, and the 
normalized magnitudes are plotted for each discrete @ to 
depict the radiation patterns. The Mathcad horizontal log- 
periodic dipole array application computes the space wave, 
surface wave, and total electric field radiation patterns and 
radiation parameters in m $-0 and $-7/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (w) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0,) as $ varies from O to 2m in 312 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 


electric field just as those for the vertical planes. 


s 


Equations 3.8 and 3.9 are used to integrate equations 
16.27 and 16.28 over the hemispherical Gaussian surface above 
the ground plane at a fixed radius (hereon: me lat pay se om minal 
total average radiated power (Puu). With the discrete values 
of the electric field and total average radiated power 
determined, the Mathcad application predicts the following 


radiation characteristics from the equations in Chapter 3: 


Ro A radiation resistance (Ohms) 

DS directivity 

EPR effective radiated isotropic power (Watts) 

io E maximum theoretical effective area (square meters) 

PR maximum theoretical effective length (meters) 

P on numerical distance (vertical polarization, 602909?) 

PENES ed numerical distance (horizontal polarization, 0-909?) 

Angle,, .elevation angle of maximum directive gain (degrees) 
The directivity (Dj, effective isotropic radiated power 
(EIRP), maximum theoretical effective area (A,,,), maximum 


theoretical effective length (lw), and elevation angle of the 
directivity (Angle,,,) are determined for both vertical planes. 
As an example, the Mathcad horizontal log-periodic array 


application was executed with the following inputs: 


number of elements ee 

length of first dipole 3.246 meters 
length of second dipole 3.732 meters 
first separation distance 1.096 meters 
radius of first element 0.00325 meters 
height of antenna feed 8.0 meters 
frequency 20-10? Hertz 
distance from the antenna 3000 meters 
relative dielectric constant 4 

ground conductivity 5-10? 
characteristic admittance 1/450 
termination impedance O Ohms 
angle from vertical to array axis SOS 
elevation angle index 295 iid Lo 
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HEENnDIessIoxrhrougueo v are thefspace wave and surface wave 


radiation patterns in the $-0 and $-7/2 vertical planes, and 


the designated horizontal plane for this example. 


EN 


FIGURE 16.2: Horizontal log- 
periodic space wave radiation 


pattern for $20 vertical 
plane. 
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The following radiation parameters were predicted by 


the Mathcad application for a sinusoidal current input of one 


Amp. The predicted radiation parameters are: 

Total power radiated (Watts) i5 546 
Radiation resistance (Ohms) roi oo 
Numerical distance (vertical) 93.742 
Numerical distance (horizontal) 3399. 15 

=0 ZNA 
Directivity 0.838 18.384 
EIRP (Watts) | 63-301 138.882 
Max eff area (sq meters) 152003 OS 
Max eff length (meters) 4.904 227971 
Ang le e degi CES) 89.7 24.0 


These results are consistent with expectations for this 
particular Foon igoa aron Appendix E contains computer 
hardcopies of additional example calculations for the 
horizontal log-periodic dipole array and compares predicted 
radiation parameters to those expected based on previous 


calculatqyons 
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XVII. THE HORIZONTAL YAGI-UDA ARRAY 

The orientation of the horizontal Yagi-Uda array is 
G@epmeeea im Figure 17:1, where l,, is the length of the 
current driven element, NR is the number of parasitic 
reflector elements, ND is the number of parasitic director 
elements, n is an index equal to 0,1,2..NR+ND, and 1, is the 
length of the n?«1 element starting with the outermost 
reflector and counting consecutively toward the opposite end. 
The d, terms are the separation distances between the n"«1 and 
n?-«2 elements, Hy is the height of all antenna elements above 
ground, R is the radial coordinate, @ is the elevation 
coordinate, and $ is the azimuth coordinate. 

The Mathcad application for the horizontal Yagi-Uda array 


requires the following inputs: 


DE. oos number of reflector elements 
NDS. number of director elements 
MA. T length of the n*+1 element (meters) 
o Ss. distance between n%+1 and n%+2 elements (meters) 
radi.. radius of the n*+1 element (meters) 
HUNE oan antenna height above ground plane (meters) 
[E sous operational frequency (Hertz) 
FEE. uns distance from array (meters) 
A relative dielectric constant of ground plane 
E 1s conductivity of ground plane 
(m PS elevation angle index (from Table 3.2) 
Unlike the  log-periodic applications, all element 


lengths, radii, and separation distances must be input 
manually on page two of the application. Consult a Mathcad 


manual to be certain the entries are made correctly. The user 


237 


MO Area dor 10 Hey Sbig-Teen 1EJUOZTIQH 


C - CIN-- HN 


Ex s. 





'Sjuouo[e ÁAeiie TTE 
JO au5reu eya st 





MALTA CVHHHAAO 


Spatial orientation of the horizontal Yagi-Uda 


array for its corresponding Mathcad application. 
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138 


inputs the frequency for which the radiation parameters 
discussed in Chapter 3 are calculated. As noted in the first 
paragraph, indices start at zero instead of 1, so the index 
Simply refers to the index-plus-one position in terms of 
successive values for a given parameter. This is to maintain 
continuity with Mathcad which also indexes from Zero. The 
dipole half-lengths (h,) are found from 1,/2. Antenna height 
(H) can be any value greater than or equal to zero, and the 
distance from the array (R) must meet the far-field 
requirements of Chapter 3.C. The elevation angle index (w) 
sets the @ coordinate for which a horizontal radiation pattern 
is determined. Table 3.2 lists possible indices and their 
corresponding elevation angles from .57? to 89.1? in 
increments of about 2.3°. Indices between those listed can be 
used to interpolate a better approximation of a desired 
elevation. 

The open circuit impedance matrices are calculated to 
determine the current distribution among the dipole elements. 
The matrix of mutual impedances between the dipole elements is 


represented as 


Zo o Lo,1 Lo, 2 mt Zo, (N-1) 
Zio La 1 L1 2 GENET T Zi, (N-1) 
42,0 22,1 22,2 s 22, (N-1) 
[Z] = ! —— 7 
Z (N-1),0 Z (N-1),1 Z (N-1),2 ee 4 oœ zr 
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Each main diagonal term (Z,,) is the self impedance of the i%+1 


elementi The off-diagonal terms (Z) are the mutual 
impedances between the i®+1 and k®*+1 elements. Subscripts i 
and k are matrix indices (0,1,2....NR+ND). 


The matrix of open circuit mutual impedances between the 


actual elements and their images is represented as 


ZL Zio, ¿lo 2 IZ o ZIg.(N-1) 
Zl... Gling ZE A AR 
[zr] E Llao ZI» ZI, > oe 8 8 og ZI, (w-1) (17.2) 
Ll gays £2i(q-13)5,3 €i (w-30,2 * * 5 -* É£Íqwin,(Qü-1) 
The main diagonal terms ZI). are che mucus ame dances 


between the i%+1 element and its image, and the off-diagonal 
terms (Z1,,) are the mutual impedances between the i%+1 element 
and the image of the k"+1 element. 

The main diagonal self-impedance terms of equation 17.1 


are found from [Ref 6: pp. 205-206] 


60 


z O SS se 
26^ deseos DSL) 


I [O(U,) -20(U,)] (a 99) 
+e? (O( Vy) -20( V4) }+2[0(U",) -0(9) -9(V.) ]|-20 (07) [50s (285 ] ] 


where Q(x) =Ci(x) -j-Si (x) = SOSU) ay- j.f 3380. ay (17.4) 
0 


8 — X 


U,=B |/2rad?+4h?-2h,| U',=/2Brad, Vo=B |/2rad?+4h?+2h,] 
u,= | /2rad?+h?-h,| V,-B |/2rad?+h?+h,| 
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The off-diagonal mutual impedance terms of equation 17.1 


are given by [Ref 11] 


= PRENNE " 
Es kr cos[B(h,-h,)]-cos[B(h,+h,) ] er! 


[e 7i *9to (0) -o(0) -Q(05) j«e7 P n^? t0(v.) -O(V,) -O( V5) ] 
re PERO (0-9 (V,) -Q(U,) LN OD -Q(U.) -QU) | 


+20(BD; ,) (cos [B (h,-A,) ]+cos [B (A;+h,) J) 


where ME EEES (Xx) = [ S28 ay - q E ay 
0 


UB |D + (2 +A (a AD] — vp |/D£ (a +A) (A A) | 
0.8 [0 (a ay aja] vB |/DZ (hyp) 2+ (h;-h,) | 
0,28 [Di vrhi-a)] v,7B [/pz ibit] 


The D, terms are the distances between the i,*1 and k,-«1 


dipole elements and are given by 
k-1 
D, ,7 Dy, ;7 9. d, (EOD 
n=1 


The mutual impedances between the actual elements and 
their images must also be computed to determine the elemental 
current distribution. The mutual impedance terms of equation 


17.2 (21;,) are given by equation 17.5, except the D, terms in 
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the Q(x) arguments are replaced by di, terms which represent 
the distance from the 1"+1 element to the image of the k%+1 


element and are given by [Ref 6: pp. 205-208] 
di; ,=VHo+Di, x Giga) 


The function Q(x) given by equations 17.4 is defined in 
terms of sine [Si(x)] and cosine [Ci(x)] integrals. Mathcad 
is incapable of evaluating Si(x) and Ci(x) directly, so a 
series expansion and polynomial approximation are written into 
the Mathcad code for both the sine and cosine integrals to 
evaluate the functions Q(x). The arguments for Q(x) are 
always real, so it is unnecessary to find Si(x) and Ci (x) 
expressions valid for complex numbers. 

For arguments less than one, Si(x) and Ci(x) are 


evaluated by series expansions given by [Ref 9: p. 232] 


P ES (-1)2 x20* 

c da (2n+1) (2n+1)! SN 
- ET (Pe yen 

GIC OE mec) DZ Te (17-9) 


where y=0.5772156649 


For arguments greater than one, Si(x) and Ci(x) are 
expressed in terms of auxiliary functions, f(x) and g(x), 


given by Ref 9: p.232] 


Si (x) its cos (x) -g (x) sin(x) (172910) 


Ciox)fx)saintciecicos aA (1.7 Salil) 
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The auxiliary functions f(x) and g(x) are evaluated by the 


polynomial approximations [Ref 9: p. 233] 
y Y AA ei) ae 
lear XS x b x “+b, 
where Leoa 
a,=38.027264 b,=40.021433 
a 20959957035 p.=322 62491 
E O O b,=570.236280 
a,=38.102495 b,=157.105423 
EN ES e. NES (17.13) 
E DAA AD, AD, 
where [E Co ME EID! 
a,=42.242855 y= 485196927 
aa oer) SOS 5b, =482 .485984 
a,=352.018498 b,=1114.978885 
a coded og b,=449.690326 


When the mutual impedances are found, the array's base 


current distribution is computed from [Ref 10 p. 258-259] 


[TB] =([2] +P, ¿[211y* [VREF] Ee 
, where 
1- 
Mo Ts | n 


and n is the complex index of refraction. The entries of the 


[IB] matrix are the base currents of the dipole elements, and 
[VREF] is the voltage matrix to which the [IB] distribution is 
referenced. Since the Yagi-Uda has only one driven element 
(corresponding to the NR” index), and since a sinusoidal 


voltage response with a maximum of unity is assumed across the 
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input terminals, [VREF] is a 1 X NR+ND matrix with a 1 as the 
NR” entry and NR+ND-1 zeroes in the remaining positions. The 
electric field equations are referenced to the current maxima 
values for the dipole elements, [I], which are determined from 


the [IB] matrix using [Ref 6 p. 209] 


y. (1D) 
sin(Ph,) 

The mutual impedance calculations are not valid when any 
element length is an exact integer multiple of the wavelength. 
When this occurs, there is a singularity error in the mutual 
impedance calculations, and it iS necessary to vary the 
frequency such that no element is exactly an integer multiple 
of the wavelength. The change required is only a percent or 
two, and the predicted radiation parameters are still a good 
estimate of those for the original frequency. 

The electric field for the horizontal Yagi-Uda array is 
obtained for each element in a manner analogous to that used 
for the vertical dipole in Chapter 4. The equations for the 
individual elements are simplified into two expressions for 
the array’s electric field components by the array factor. 


The Yagi-Uda’s array factor is given by [Ref 6: p. 204] 


NR+ND-1 
s= Y rIjcos(Bh(sin(0)sin($))-cos(ph,)] e? 6n 9*9 (17.16: 
1=0 


where Y, is the y coordinate for the i*"«1 element given by 


E 
Y=% d, (E) 
n=0 
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The array factor is multiplied by the element factors to 
find the equations for the total radiated electric field 


distribution of the horizontal Yagi-Uda array, 


! -JBR 5 
E=-j60 e COS (0) e JP Pocos (0) fe Cigars) 


R 1-sin? (6) cos? (o) 


[cos (0) (1-T Cd. c 909 4 (1-T )F @ ~72BH cos (8) , 
v v e 


127-8127 07 (5155 (o) YB Si cos (0) | 
de ras 
o-JPR sin (d) e Picos (8) (S) x 


(17.19) 
R 1-sin? (6) cos? (6) 





E,-j60 
E De -j2pHscos (8) T (1-T,) Eu cy PE 2j 


The first two terms inside the brackets of each equation are 
the space wave, the third terms are the surface wave, and I, 
is the current term for each element based on a sinusoidal 
voltage across the input terminals with a maximum of unity. 
The requested inputs are used to calculate the following 
variables using a constant $ of $20 and $-T7/2 for 312 discrete 


values of 0 which are equally incremented from -T7T/2 to m/2: 


p ds: half-length of the i%+1 element (meters) 

O y Coordinate of the i"«1 element 

TIMES. . current for the i?«1 element 

Agua . wavelength of the operational frequency (meters) 
D... free space wavenumber for operational frequency 

n... WX index of refraction 

DP. E. complex numerical distance (vertical polarization) 
Pao. —-! complex numerical distance (horizontal polarization) 
IDE a vertical reflection coefficient 

y . ues horizontal reflection coefficient 

ID ln vertical surface wave attenuation factor 

Foe. horizontal surface wave attenuation factor 
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The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 with $20 and $-7/2. The space wave and surface wave results 
are combined for corresponding values of @ to obtain the total 
radiated electric field distribution for the ¢=0 and d=1/2 
vertical planes. The space wave, surface wave, and total 
electric field results are then normalized with respect to the 
maximum field intensity of each, and the normalized magnitudes 
are plotted for each discrete @ to depict the radiation 
patterns. The Mathcad horizontal Yagi-Uda array application 
computes the space wave, surface wave, and total electric 
field radiation patterns and radiation parameters in the ¢=0 
and ¢=7/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (w) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0,) as ¢ varies from 0 to 2m in 312 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

Equations 3.8 and 3.9 are used to integrate equations 
17.18 and 17.19 over the hemispherical Gaussian surface above 
the ground plane at a fixed radius (R) from the array to find 
total average radiated power (Paa). Since the Yagi-Uda current 
distribution is referenced to a sinusoidal voltage input with 


a maximum of unity, the radiation resistance can not be 
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calculated from equation 3.10 as in previous applications. 
Instead, radiation resistance must be found with respect to 


the unity voltage input using 


NIS gut 


P OMe aie. = 
rad DU rad E 


O 
With the discrete values of the electric field and total 
average radiated power, the Mathcad application predicts the 


following radiation characteristics from Chapter 3: 


5r cee radiation resistance (Ohms) 

Ip Cn coss directivity 

A effective radiated isotropic power (Watts) 
AAA maximum theoretical effective area (square meters) 
lee ic maximum theoretical effective length (meters) 
Do ee numerical distance (vertical polarization, 02909?) 
Se numerical distance (horizontal polarization, 0-90?) 


Angle,, .elevation angle of maximum directive gain (degrees) 


The directivity (Dj, effective isotropic radiated power 
(EIRP), maximum theoretical effective area (A,,,), maximum 
theoretical effective length (1,,), and elevation angle of the 
directivity (Angle,,,) are determined for both vertical planes. 

As an example, the Mathcad horizontal Yagi-Uda array 


application was executed with the following inputs: 


number of reflectors 1 
number of directors 1 
element lengths 15.6, 15, and 14 meters 
element separations 7.5 and 6.5 meters 
element radii 0.001, 0.001, and 0.001 meters 
height of antenna 8.0 meters 
frequency 10:10? Hertz 
distance from the antenna 3000 meters 
relative dielectric constant 72 
gcound conductivity 4 
elevation angle index 285 (17°) 
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Figures 17.2 through 17.7 are the space wave and surface wave 


radiation patterns in the $20 and $-7/2 vertical planes, and 


the designated horizontal plane for this example. 
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The predicted radiation parameters based on a one volt 


response across the input terminals are: 


Total power radiated (Watts) 0.00482 

Radiation resistance (Ohms) ISO 

Numerical distance (vertical) 0.04363 

Numerical distance (horizontal) O 

=0 =1/2 

Directivity meo 14.631 

EIRP (Watts) 02008 00705 

Max eff area (sq meters) 120.. 26 1047.9 

Max eff length (meters) 14955506 33. 96 

Angle, (degrees) EI 41.714 

These results are consistent with expectations for this 

Pparticulam™ configuratioB*A Appendix F contains computer 
hardcopies of additional example calculations for the 
horizontal Yagi-Uda array and compares predicted radiation 
parameters to those expected based on previous calculations. 
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XVIII. REMARKS AND CONCLUSION 

This thesis is the culmination of nine months of research 
and computer programming. No new electric field equations 
were derived exclusively for this project. It was directed by 
NAVMARINTCEN prior to commencing that existing electric field 
equations would provide the basis for the radiation parameter 
predictions. The expressions used in this thesis are all 
previously derived by pioneers in antenna radiation theory 
such as Sommerfeld, Norton, Cox, King, Ma, and Walters. The 
equations are compiled here solely for the purpose of 
crediting the source references and to describe how the 
Mathcad code was assembled. 

The accuracy of the predicted radiation characteristics 
is totally dependent upon the extent to which the equations 
used realistically model the actual radiated electric fields 
of the antennas. It is near impossible to obtain analytic 
results which accurately model an  antenna's radiation 
characteristics in all cases. There are too many operational 
and environmental variables to obtain one general, all-purpose 
expression. Even if all possible variables could be accounted 
for in a single expression, computer processing times would be 
unacceptably long. The equations used herein are simple 
enough to be evaluated by Mathcad in a reasonable length of 
time, but the radiation parameter predictions are accurate 
enough to be useful for antenna analysis. 


150 


As addressed in appendices A-F, the equations presented 
in previous chapters and evaluated by the Mathcad applications 
provide radiation parameters consistent with expectations for 
the inputs which have been executed to date. All antenna 
configurations which have been computed exhibit radiation 
characteristics consistent with other computational programs 
and empirically obtained patterns and parameters. There is no 
way to test the accuracy of the Mathcad predictions other than 
empirical measurements for each antenna configuration, clearly 
a task beyond the scope and purpose of this report. However, 
adequate analysis of the Mathcad results has been provided to 
demonstrate that the applications provide very good estimates 
of antenna radiation parameters as a base level analysis tool. 

Most of the applications have computations times under 
ten minutes. However, even with the simplified equations used 
by Mathcad, some of the applications can take upwards of two 
hours to compute on a 33 MHz 80386 PC. The applications for 
BucNrHuonbrc  deuble FRembold” and array (with more than four 
to five elements) antennas can take an exorbitantly long time 
for power calculations. For this reason, the applications 
Should be run on the fastest PC available; a 50 MHz 80486 is 
preferable. Computation times can also be reduced by 
computing only the radiation patterns until it is determined 
which frequencies are of greatest interest. Then the average 


power calculations for those frequencies can be executed. 
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The Mathcad applications were written for NAVMARINTCEN to 
provide antenna analysts with the capability to predict 
radiation parameters based solely on antenna physical 
dimensions and ground properties. The Mathcad results 
computed thus far are consistent with expectations and are 
most likely providing accurate predictions of actual antenna 
radiation characteristics. Only extended use of the 
applications and empirical confirmation of results will prove 
out the accuracy of the Mathcad code. Results obtained thus 
far certainly justify continued use of the Mathcad code as an 


antenna analysis tool. 


152 


APPENDIX A: 
VERTICAL DIPOLE ARRAY COMPUTER OUTPUT 


This appendix contains computer hardcopies from the Mathcad 
vertical dipole application which show the input values and 
predicted radiation characteristics for two sample calculations. 
The configuration is a half-wave dipole at one quarter wavelength 
above soil (e€,=10 and 0=10%) for the first example and the same 
configuration above seawater (€,=72 and o=4) for the second. 
Reference 6 [p. 91] provides the radiation patterns and gain 
predictions from several sources for the configuration in the 
first example. 

The elevation of maximum directive gain is slightly higher 
for the Mathcad output than for those given in reference 6, but 
the overall radiation patterns are very Similar. The 4.17 (6.2 
dB) value of directivity for Mathcad is quite a bit higher than 
the maximum gain values of about zero dB in reference 6. 
However, Since a half-wave dipole has a free-space directivity of 
1.64 (2.1 @B), one would expect the actual directivity to be 
closer to the Mathcad prediction because of the effect of the 
reflected wave (constructive and destructive interference) on the 
space wave, and the lower total average radiated power resulting 
from ground plane losses. The seawater example yields results 
consistent with expectations. With respect to the soil example, 
the seawater example’s directivity is slightly higher due to a 
Stronger reflected wave, and the surface wave iS stronger at 
grazing angles (#=90°) due to higher conductivity. 
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VERTICAL DIPOLE 


This application calculates far field radiation patterns and parameters associated 
with vertical thin-wire dipole antennas (diameter << wavelength). The antenna is 
mounted vertically along the z-axis in a rectangular coordinate system with the feed 
at the antenna's center at a set height above the surface. Required inputs are the 
antenna length, feed height above the surface, transmitted frequency, distance from 
the antenna, the conductivity and dielectric constant of the surface below the 
antenna. The planar earth model is assumed tn predicting radiation patterns. 
Predicted operating frequencies assume that the antenna is a quarter- wavelength. 
half- wavelength, three-quarter-wavelength, or full-wavelength dipole. A sinusoidal 
current input with a maximum of unity is assumed. All radiation patterns are 
normalized with respect to the maximum electric field intensity transmitted by the 
antenna. Plotted radiation patterns are valid for any vertical plane passing through 
the antenna axis. The electric field magnitudes to which the Radiation Patterns are 
normalized are displayed below their respective plots. Polarization is vertical for all 
vertical dipoles. 


Input the Dipole li 5 p 1 
length in meters 2 


Wavelengths and Frequencies 


4-1 
Ay .= 4-1 A» = 2-1 Ag = 3 ^. = l 
c c c c 
nei. ees (gc s — 
1 à 3 4 
i pe A3 Ay 
Possible Operating Frequencies (Hertz) 
z Dipole: EOS 10" : Dipole: f = 3- 10" 
Y Dipole: [4 7 4.5- 107 A Dipole: O 10" 


i Input the Distance 
Input the operating An 1n$ P 
Frequency (Hertz) Ig = 30:10 from the Antenna R -- 3000 


(meters) 
Input the ground zl Input the ground 5 " 10° 2 
Dielectric Constant Conductivity 


Input the Height of 5 
the Antenna Feed "Ome 
(meters) 


Radiation Patterns valid for any Vertical Plane passing through the Antenna 


Space Wave Radiation Pattern 





Maximum Space Wave Electric Field 
Intensity (Volts per meter) max( MagE1) - 0.01733 


"OS 


Surface Wave Radiation Pattern 





Heb, 6 Ke, 


Maximum Surface Wave Electric Field 4 
Intensity (Volts per meter) max (MagE2) = 2.69275: 10 
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Combined Space and Surface Wave Radiation Pattern 





Maximum Radiated Electric Field 
Intensity (Volts per meter) max ( MagE3) = 0.01732 


iw 


2 : 
"AREE 
cos(B:h-cos(3)) - cos(B:h) | e e 9 M n. 
sin( $) R - Hg: cos( $) R + Hy 


Ig E 
power := 30- 


2 
E3 -x RÊ. 
ap, . MED Directivity - 47 R max ( sqE3) 
2:( 120-2) power 
Radres = 2-power 
Total Power Radiated (Watts) power = 10.79029 


Radiation Resistance (Ohms) Radres = 21.58059 


Directivily (maximum Power Gain pea 
assuming 100% Antenna efficiency) Directivity — 4.16784 


Effective Isotropic Radiated 


Power (EIRP) (Watts) | Directivity: power = 44.97225 
Maximum Effective Área 2 
As ) ^: Directivit 
(Along Radial of. Directivity) (x) CR. Laid 
(square meters) 4. x 


Maximum Effective Length 
(Along Radial of Directivity) 
(meters) 


= 2.75579 





Numerical Distance for 
P = 74.95924 
Vertical Polarization | Peg | 


Elevation Angle of Maximum Angle3 = ( 


19.39778 
Power Gain (Degrees) 


19.39778 


VERTICAL DIPOLE 


This application calculates far field radiation patterns and parameters associated 
with vertical thin-wire dipole antennas (diameter << wavelength). The antenna ts 
mounted vertically along the z-axis in a rectangular coordinate system with the feed 
at the antenna's center at a sel height above the surface. Required inputs are the 
antenna length, feed height above the surface, transmitted frequency, distance from 
the antenna, the conductivity and dielectric constant of the surface below the 
antenna. The planar earth model is assumed in predicting radiation patterns. 
Predicted operating frequencies assume that the antenna is a quarter- wavelength, 
half-wavelength, three-quarter- wavelength, or full-wavelength dipole. A sinusoidal 
current input with a maximum of unity is assumed. All radiation patterns are 
normalized with respect to the maximum electric field intensity transmitted by the 
antenna. Plotted radiation patterns are valid for any vertical plane passing through 
the antenna axis. The electric field magnitudes to which the Radiation Patterns are 
normalized are displayed below their respective plots. Polarization is vertical for all 


vertical dipoles. 


Input the Dipole esas E 
length in meters 2 
Wavelengths and Frequencies 
4-1] 
Ay = 4. Ag = 2-1 Aq im 3 ^, =f 
c c c c 
f - — f = — f I — f I 
1 2 3 4 
Ay A, Az A, 
Possible Operating Frequencies (Hertz) 
z Dipole: f = 1.5: 10* > Dipole: p = 3: 10? 
^u Dipole: [4 7» 4.5: 10? A Dipole: [76 10? 


Input the Distance 


Input the operating f.. 30.108 
Frequency (Hertz) — ? ES Antenna R -- 3000 


Input the ground m-7 Input the ground 5 d 4 
Dielectric Constant Conductivity 


Input the Height of 225 
the Antenna Feed Ho S 


(meters) 
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Radiation Patterns valid for any Vertical Plane passing through the Antenna 


Space Wave Radiation Pattern 





Maximum Space Wave Electric Field 
Intensity (Volts per meter) max ( MagE1) = 0.03394 


Tol 


Surface Wave Radiation Pattern 
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Heb, 4 Ha, 


Maximum Surface Wave Electric Field 
Intensity (Volts per meter) max ( MagE2) = 0.03331 
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Combined Space and Surface Wave Radiation Pattern 





0.5 


0 
Hreb, H3, 


max(MagE3) = 0.03398 
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P 
-IB R- 8. AIBR 
cos( 8 )b-'costs@))) i acasies any mire Ie yet) e ve 


= 30: A A 
as sin( ¢) R - Hg cos(£) R + Ho 


a 
~ 


2 
E3 a Re. 
ol (| E3]) Directiviy - 57 RÀ max(sqE3) 
2- (120: x) power 


Radres := 2- power 
Total Power Radiated (Watts) power = 37.64371 
Radiation Resistance (Ohms) Radres = 75.28741 


Directivity (maximum Power Gain VEM 
assuming 10096 Antenna efficiency) Drrectivity = 4.59967 


Effective Isotropic Radiated 


Power (EIRP) (Watts) Directivity: power = 173.14876 

Maximum Effective Area | (As) 2 Directivity 

(Along Radial of Directivity) AA == 36.60304 
(square meters) 4% 





Maximum Effective Length Padre: (As) 2 Directivity 
(Along Radial of Directivity) — 2. | 77 ^2 5.49735 
(meters) 480. x? 


Numerical Distance for E 
Vertical Polarization | P eg | 0.392352 


Elevation Angle of Maximum Angle3 = ( 


8.55784 
Power Gain (Degrees) 


8.55784 


APPENDIX B: 
VERTICAL MONOPOLE ARRAY COMPUTER OUTPUT 


This appendix contains computer hardcopies from the Mathcad 
vertical monopole application which show the input values and 
predicted radiation characteristics for two sample calculations. 
The configuration is a quarter-wave monopole above soil (€,=10 and 
0-107) for the first example and the same configuration above 
seawater (€=72 and o=4) for the second. Reference 6 [p. 89] 
provides the radiation patterns and gain predictions from several 
sources for the configuration in the first example. 

The elevation of maximum directive gain is slightly higher 
for the Mathcad output than for those given in reference 6, but 
the overall radiation patterns are very similar. The 3.24 (5.1 
dB) value of directivity for Mathcad is quite a bit higher than 
the maximum gain values of about zero dB in reference 6. 
However, Since a half-wave dipole has a free-space directivity of 
1.64 (2.1 dB), the actual directivity of a quarter-wave monopole 
above a ground plane should be closer to the Mathcad prediction 
because of the effect of the reflected wave (constructive and 
destructive interference) on the space wave and the lower total 
average radiated power resulting from ground plane losses. The 
seawater example yields results consistent with expectations. 
With respect to the soil example, the seawater example’s 
directivity is slightly higher due to a stronger reflected wave, 
and the surface wave is stronger at grazing angles (#=90°) due to 
higher conductivity. 
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VERTICAL MONOPOLE 


This application calculates far field radiation patterns and parameters associated 
with vertical thin-wire monopole antennas (diameter << wavelength). The antenna 
is mounted vertically along the z-axis in a rectangular coordinate system with the 
feed at the origin. Required inputs are the antenna length, transmitted frequency, 
distance from the antenna, the conductivily and dielectric constant of the surface 
below the antenna. The planar earth model is assumed in predicting radiation 
patterns. Predicted operating frequencies assume that the antenna is an eighth- 
wavelength, quarter- wavelength, three-eighths-wavelength, or half- wavelength 
monopole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electic field intensity 
transmitted by the antenna. The electric field magnitudes to which the patterns are 
normalized are displayed below their respective plots. Radiation patterns are valid 
for any vertical plane passing through the z-axis, because the radiation pattern is 
symmetrical with respect to phi. Polarization is vertical for all vertical monopoles. 


Input the Monopole DEUS jh 
length in meters 


Wavelengths and Frequencies 


4-1 
= 4.] Aa = 2-1 Az i 3 A ist 
c c c c 
f = = fs = — fa := — fa = — 
1 2 3 4 
Ay Ao ha ^4 
Possible Operating Frequencies (Hertz) 
; Monopole: f, - 5: 10° z Monopole: f, = 1: 10? 
3. À 
T Monopole: f} = 1.5- 10" > Monopole: fy = 2-10" 
Input the operating fs = 10: 102 Input the Distance from - 3000 
Frequency (Hertz) the Antenna (meters) 
Input the ground O Input the ground o = 107° 
Dielectric Constant Conductivity: 


Radiation Patterns valid for any Vertical Plane passing through the Antenna 


Space Wave Radiation Pattern 
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Maximum Space Wave Electric Field 


Intensity (Volts per meter) max( MagE!) - 0.01115 
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Surface Wave Pattern 





Hreb, Hé; 


-4 


aximum 


M 


Surface Wave Field Intensity 


(Volts per meter) 


max(MagE2) = 8.23677: 10 
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Combined Space and Surface Wave Pattern 
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Maximum Radiated Electric Field 


Intensity (Volts per meter) max( MagE3) - 0.01114 
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Nia 


Spa ho - cos( B: h) ~ | 


+j sin( B-h-cos( ¢) ) T. 
ere laca + (-(cos(3):sin(B:h))) S 
R2 sin( ¢) 
cos( d: h:cos( $)) - cos(B.h) ... 





un sin(B-h-cos({)) 4) ` 
sin( 5) 





n2-cos(t) - {a2 - sine)? 


n?.cos( 0) ince sin( t) * 


2 
E3 a- RÊ 
sqE3, - (EAD 1) Directivity := E 
2-(120-1) power 


Radres = 2: power 


Total Power Radiated (Watts) power = 5.74454 
Radiation Resistance (Ohms) Radres = 11.48908 


Directivity (maximum Power Gain Bon 
assuming 100% Antenna Efficiency) Directivity — 3.23898 


Effective Isotropic Radiated 


Power (EIRP) (Watts) Directivity: power = 18.60648 
Maximum Effective Area ZI A e 
| Ag)” Directivit 
(Along Radial of Directivity) ( 5) di y = 231.97519 
(square meters) 4-1 


| ^ 
Maximum Effective Length E (Ag) *- Directivity 
M ME ^ WE 4 


(Along Radial of Directivity) 2 SS 
(meters) 480- x 

Numerical Distance for E 

Vertical Polarization | Peg] S 

Elevation Angle of Maximum Rock | = 25.67353 

Power Gain (Degrees) 25.67353 
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VERTICAL MONOPOLE 


This application calculates far field radiation patterns and parameters associated 
with vertical thin-wire monopole antennas (diameter << wavelength). The antenna 
is mounted vertically along the z-axis in a rectangular coordinate system with the 
feed at the origin. Required inputs are the antenna length, transmitted frequency, 
distance from the antenna, the conductivity and dielectric constant of the surface 
below the antenna. The planar earth model is assumed in predicting radiation 
patterns. Predicted operating frequencies assume that the antenna is an eighth- 
wavelength, quarter-wavelength, three-eighths-wavelength, or half-wavelength 
monopole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electic field intensity 
transmitted by the antenna. The electric field magnitudes to which the patterns are 
normalized are displayed below their respective plots. Radiation patterns are valid 
for any vertical plane passing through the z-axis, because the radiation pattern is 
symmetrical with respect to phi. Polarization is vertical for all vertical monopoles. 


Input the Monopole h = 75 o 
length in meters 


Wavelengths and Frequencies 


4-1 

Ay = 4.1 Xo =e Az = E Ng = | 
c c c c 
[-— [5 :- — fa := — fa = — 

1 2 3 4 
M i 23 ^a 

Possible Operating Frequencies (Hertz) 

> Monopole: fq = 5.10 4 Monopole: Í, = 1-10? 
7 


| 


- Monopole: f3 = 1.5- 107 Monopole: PA = 2-10 


Input the operating fs = 10: 6 Input the Distance from - 3000 
Frequency (Hertz) ie the Antenna (meters) 
Input the ground a= 7 Input the ground 2a 


Dielectric Constant o Conductivity: 


173 


Radiation Patterns valid for any Vertical Plane passing through the Antenna 


Space Wave Radiation Pattern 
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Maximum Space Wave Electric Field 
Intensity (Volts per meter) max( MagE1) - 0.01855 
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Surface Wave Pattern 
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Heb. Hé. 


aximum 


M 


Surface Wave Field Intensity 


(Volts per meter) 


max( MagE2) - 0.01976 


yes 


Combined Space and Surface Wave Pattern 
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Maximum Radiated Electric Field 
Intensity (Volts per meter) max( MagE3) = 0.01976 
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+ j-| sin( B-h-cos(¢)) ... 
NET L (t C (cos(S)'sin(B-h))) j J 
R? sin( £) 


" n?-cos(£) - |n? - sin(¢)* 


n?.cos( () Poné sin( ¿y 


X 
- VP sio - cos(B- h) " | 





cos( B: h:cos($)) - cos(B:h) ... 
+ Ee» a 
L_ Ut 1+ (-(cos(§)-sin(s-b)) 
sin(¢) 











ney» (ED Dedi x formano) 
Radres - 2-power 
Total Power Radiated (Watts) power = 16.19454 
Radiation Resistance (Ohms) Radres = 32.38909 


Directivity (maximum Power Gain NM. 
assuming 10096 Antenna Efficiency) Directivity = 3.61503 


Effective Isotropic Radiated 


Power (EIRP) (Watts) Directivity: power = 58.54372 
Maximum Effective Area 2 nisiut 
Ac) Directivit 
(Along Radial of Directivity) Mom Lin: 
(square meters) 4-1 
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(Along Radial of Directivity) = 9.4327 


| ^ 
Maximum Effective Length : = (As) ^- Directivity 
(meters) 


480- x * 


Numerical Distance for _ 
Vertical Polarization | P € | 0.04363 


Flevation Angle of Maximum a 
Power Gain (Degrees) Angli? | 7 


APPENDIX C: 
HORIZONTAL DIPOLE ARRAY COMPUTER OUTPUT 


This appendix contains computer hardcopies from the Mathcad 
horizontal dipole application which show the input values and 
predicted radiation characteristics for four sample calculations. 
The configuration is a half-wave dipole at one- quarter, one-half, 
and three-quarters wavelengths above soil (¢€,=10 and o=107) for the 
first three examples, respectively, and a half-wave dipole at one- 
quarter wavelength above seawater (€,=72 and 0=4) for the fourth. 
Reference 6 [pp. 86-88] provides the radiation patterns and gain 
predictions from several sources for the configurations in the 
first three examples. 

The radiation patterns and maximum directive gain computed by 
the first three Mathcad examples are almost identical to the those 
given in reference 6. The directivity predictions are slightly 
higher for Mathcad, but the overall similarity between the 
predicted radiation characteristics 1S noteworthy. The seawater 
example also yields results consistent with expectations. With 
respect to the soil example, the seawater example’s directivity is 
slightly higher due to a stronger reflected wave, and the surface 
wave is slightly stronger at grazing angles (090?) due to higher 
conde ME —s FOn worl zontal polarization, the higher conductivity 
surface below the antenna does not result in a greatly enhanced 
surface wave and increased directivity as it does for vertical 


polarization, 


t79 


HORIZONTAL DIPOLE 


This application calculates far field radiation patterns and parameters associated with 
horizontal thin-wire dipole antennas (diameter << wavelength). The antenna is 
mounted above and parallel! to the x-axis in a rectangular coordinate system. The 
feed is at the center of the antenna at a set height directly above the origin. 
Required inputs are the antenna length, feed height above the surface, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. Predicted operating frequencies assume that the antenna ts a 
quarter- wavelength, half- wavelength, three-quarter- wavelength, or full wave- 
length dipole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect (o the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. The radiation patterns are plotted for the phi-0 and phi-pi/2 
vertical planes parallel and perpendicular to the x-axis respectively. A horizontal 
radiation pattern is plotted at an elevation selected by the index from the elevation 
angle index table. Polarization is a combination of vertical and horizontal 
depending upon spatial orientation relative to the dipole. 


Input the dipole f= 15 h- I 
length in meters 2 


Input the index of the elevation angle for 
which to calculate the horizontal radiation d = 535 
pattern (from the angle index table) 


Wavelengths and Frequencies 


4-1 
Ay a 4-1 Ao = 2-1 Àa = 3 Ng = l 
C C C C 
[,- [A e me ra 
3 4 
M $0 A9 M 


Possible Operating Frequencies (Hertz) 


Dipole: io = 1: 10? 


M 6 
—  Dipole: f, = 5:10 
4 ipo 1 

X Dipole: [4 7 15-10" A Dipole: [72 10? 
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Input the operating - 10-109 esu ELR 
Frequency (Hertz) ie 3 f - p As 
Input the Height of He Input the ground - 30 
Antenna Feed (meters) 0 Dielectric Constant 
Input the Distance H Input the ground 2. 3.10 ? 
from Antenna (meters) Wo song Conductivit y: B 
Index of Refraction 
Rd; = R - Hg:cos(6j) — ]B000.o 
"n ses) 
fs 107 $) 
R = R+ Hg cos(6;) S 
2 
; n? 
Complex Numerical Distance Pet Ne a os ( 0. ) + m - 
[or Vertical Polarization ). 
sin E 
Complex Numerical Distance E 0 y 
for Horizontal Polarization ran s Bes (4) * (sin )^) 
(n?.cos( 4) | - ( n? - sin (6j) 2) 

Vertical Reflection Coefficient RE E ELEM uM 


Horizontal Reflection Coefficient 








(n2. cos(8) ) s ( n? - sin(6j) 2 


cos (6j) + (h? - "— 
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Radiation Patterns in Phi=Pi/2 Plane (Perpendicular to Dipole) 
Space Wave Radiation Pattern (Phi=Pi/2) 





-1 -0.5 0 0.5 1 
Heb HIP, 


Max E-Field Intensity (Volts per meter) max ( MagE1P) = 0.03632 
Surface Wave Radiation Pattern (Phi=pi/2) 





net, ; HEP, 


Max E-Field Intensity (Volts per meter) — max( MagE2P) - 1.03856: 10 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 





Heb, xSP, 


Max E-Field Intensity (Volts per meter) max( MagE3P) - 0.03632 
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Radiation Patterns in Phi=O Plane (Parallel to Dipole) 
Space Wave Radiation Pattern (Phi=0) 





Rreh.» #10, 


Max E-Field Intensity (Volts per meter) max( MagE10) = 0.03594 
Surface Wave Radiation Pattern (Phi=0) 





Max E-Field Intensity (Volts per meter) — max( MagE20) - 3.84453: io * 


184 


Combined Space and Surface Wave Radiation Pattern (Phi=0) 





utet, i 430, 


Max E-Field Intensity (Volts per meter) — max( MagE30) = 0.03594 


roe 


Radiation Patterns in Horizontal Plane 


Space Wave Radiation Pattern 


Phi=0 





Heb, RIE, 


Max E-Field Intensity (Volts per meter) ^ max( MagEIH) - 0.02634 


Surface Wave Radiation Pattern (Horizontal Plane) 





Max E-Field Intensity (Volts per meter) max( MagE2H) - 2.62147- 10° 
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Combined Space and Surface Wave Radiation Pattern (Horizontal Plane) 


yeb 


Phi=0 





Max E-Field Intensity (Volts per meter) — max( MagE3H) - 0.02634 
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Total Power Radiated (Watts) power = 37.80286 


Radiation Resistance (versus 
Maximum Antenna Current) 
(Ohms) 


Radres 


75.60572 


Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP = 5.23417 DirectivityO = 5.12535 


Effective Isotropic Radiated Power (EIRP) (Watts) 


Phi = pi/2 Plane Phi = 0 Plane 


DirectivityP- power = 197.86674 DirectivityO- power = 193.75275 


Maximum Effective Area (Along Radial of Directivity) 
(square meters) 


Phi = pi/2 Plane Phi = 0 Plane 
2 ? rct 2 . S 
he ) ^ Directivit yP he)“: Directivity0 
ie = 374.87006 E = 367.07586 
«X X 
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Maximum Effective Length (Along Radial of Directivity) 
(meters) 


Phi = pi/2 Plane Phi = 0 Plane 


Radres: (As) . DirectivityP Radres- ( As) ê. DirectivityO 


= 17.34132 2- = 17.1601 





480- x* 480. 22 
Numerical Distances 
Vertical Polarization Horizontal Polarization 


| Peg| = 5.04617 | Pmg| = 1.92562- 10% 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi 7 pi/2 Plane Phi = 0 Plane 
65.8954 89.85737 
AngleP = Angled = ( | 
65.8954 89.85737 
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HORIZONTAL DIPOLE 


This application calculates far field radiation patterns and parameters associated with 
horizontal thin-wire dipole antennas (diameter << wavelength). The antenna is 
mounted above and parallel to the x-axis in a rectangular coordinate system. The 
feed is at the center of the antenna at a set height directly above the origin. 
Required inputs are the antenna length, feed height above the surface, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. Predicted operating frequencies assume that the antenna is a 
quarter- wavelength, half-wavelength, three-quarter- wavelength, or full wave- 
length dipole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. The radiation patterns are plotted for the phi=0 and phi=pi/2 
vertical planes parallel and perpendicular to the x-axis respectively. A horizontal 
radiation pattern is plotted at an elevation selected by the index from the elevation 
angle index table. Polarization is a combination of vertical and horizontal 
depending upon spatial orientation relative to the dipole. 


Input the dipole I- 75 - 1 
length in meters 2 


Input the index of the elevation angle for 
which to calculate the horizontal radiation d - 535 
pattern (rom the angle index table) 


Wavelengths and Frequencies 


4-1 
Ay - 4.] Az 22] hy = 3 Ng = | 
C C C C 
f = — f um —— f M > f m => 
1 2 3 4 
ry Az Az Ay 


Possible Operating Frequencies (Hertz) 


Dipole: 155-2: 10° 


À ? 
— Dipole: f, = 1:10 
4 Dipole i 


= Dipole: fy = 3-10" A Dipole fy = 4-10" 


Input the operating 


N 
= 






6 c 
= 20-10 Ae = = = 
Frequency (Hertz) 3 mf 5 E As 
Input the Height of H. = 7.6 Input the ground - 30 
Antenna Feed (meters) 0 — Dielectric Constant 
Input the Distance - 3000  inputthe ground . 13.1 2 
from Antenna (meters) P. Conductivity: | AL 
Index of Refraction 
Rd, = R- Hg: cos ( 6;) | E 
n: [tr - j:[z2—— —— 
-6 
fz- 10 
Rr, = R + Hg cos( 6j) 5 ) 
: n* 2 
Complex Numerical Distance pame -p B: Rr ness (6, ) + - sin € 
for Vertical Polarization l 2-sin a? ' 


Complex Numerical Distance p 
for Horizontal Polarization 


Vertical Reflection Coefficient 


Horizontal Reflection Coefficient 





s || eos( 64) « ln? - (sin (,) 7) jl 
" (n? cos(8) ) - ( n? - sin (5) 

| 

(a2.cos(0;) ) " i -— 


P cos (6,) - (fa? sin (4)*) 
cos ( 6) + (h? - sin (02) 


Radiation Patterns in Phi=Pi/2 Plane (Perpendicular to Dipole) 
Space Wave Radiation Pattern (Phi=Pi/2) 





-1 -0.5 0 0.5 1 
Hel, HIP, 


Max E-Field Intensity (Volts per meter) max( MagE1P) - 0.03737 
Surface Wave Radiation Pattern (Phi-pi/2) 





Max E-Field Intensity (Volts per meter)  max( MagE2P) - 8.03309: io? 
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Radiation Patterns in Phi=O Plane (Parallel to Dipole) 


Space Wave Radiation Pattern (Phi=0) 





-] -0.5 0 0.5 1 
Heb. #10, 


Max E-Field Intensity (Volts per meter) max( MagE10) = 0.01561 
Surface Wave Radiation Pattern (Phi=0) 





Max E-Field Intensity (Volts per meter)  max( MagE20) = 1.54902: TE 


49:3 


Radiation Patterns in Horizontal Plane 


Space Wave Radiation Pattern 


Phi-0 





Hrebo : HIE, 


Max E-Field Intensity (Volts per meter) max( MagE!H) = 0.03734 
Surface Wave Radiation Pattern (Horizontal Plane) 





-1 -0.5 0 05 1 
wet, : 2H, 


Max E-Field Inensity (Volts per meter) — max( MagE2H) - 1.62046- 10 ^ 
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Total Power Radiated (Watts) power = 27.68975 


Radiation Resistance (versus 
Maximum Antenna Current) 


(Ohms) 


Radres = 55.3795 


Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP = 7.56408 DirectivityO = 1.32021 
Effective Isotropic Radiated Power (EIRP) (Watts) 
Phi = pi/2 Plane Phi = 0 Plane 


DirectivityP- power = 209.44744 DirectivityO: power = 36.55637 


Maximum Effective Area (Along Radial of Directivity) 
(square meters) 


Phi = pi/2 Plane Phi = 0 Plane 
2 > Means 2 . © s 
Ae ) ~: DirectivityP Ae ) ^ DirectivityO 
ae = 135.43431 p" = 23.63833 


Maximum Effective Length (Along Radial of Directivity) 
(meters) 


Phi = pi/2 Plane Phi = 0 Plane 


Radres: (Ag) * DirectivityP Radres: (Ag) “- Directivityo 





= 8.92079 2- = 3.7269 





480-2 480. x * 
Numerical Distances 


Vertical Polarization Horizontal Polarization 
= _ 4 
| Peg| = 15.283 | Pm | = 2.4896- 10 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi = pi/2 Plane Phi = 0 Plane 
.24089 41.64818 
AngleP = ( E AngleO - 
28.24089 41.64818 


HORIZONTAL DIPOLE 


This application calculates far field radiation patterns and parameters associated with 
horizontal thin-wire dipole antennas (diameter << wavelength). The antenna is 
mounted above and parallel to the x-axis in a rectangular coordinate system. The 
feed is at the center of the antenna at a set height directly above the origin. 
Required inputs are the antenna length, feed height above the surface, transmitted 
frequency. distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. Predicted operating frequencies assume that the antenna is a 
quarter- wavelength, half-wavelength, three-quarter- wavelength, or full wave- 
length dipole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. The radiation patterns are plotted for the phi=0 and phi=pi/2 
vertical planes parallel and perpendicular to the x-axis respectively. A horizontal 
radiation pattern is plotted at an elevation selected by the index from the elevation 
angle index table. Polarization is a combination of vertical and horizontal 
depending upon spatial orientation relative to the dipole. 


Input the dipole 1-5 s 


i 
length in meters 2 


Input the index of the elevation angle for 
which to calculate the horizontal radiation d - 535 
pattern (from the angle index table) 


Wavelengths and Frequencies 


4-1 
^ E 4-] Ao = 2-1 A = = Ng = | 
C C C C 
f ez — f M — f = =— f = — 
1 2 3 4 
A ^2 Az Ay 


Possible Operating Frequencies (Hertz) 


> Dipole: fy = 1.5: 10? Dipole: D-3 10° 
3. À 


4 Dipole lud» 10? A Dipole: i= 6-10 








Input the operating 2 30.109 m: _ 28 
Frequency (Hertz) Isis DO s fs d E 
Input the Height of Ha = 7,6 Input the ground er := 30 
Antenna Feed (meters) 0 ~~ Dielectric Constant | 
Input the Distance - 3000 Input the ground mU o2 
from Antenna (meters) ES Conductivity: a 
Index of Refraction 
Rd = R - Hg: cos ( 6j) 
ii = 
Rr, :- R + Hg cos( 6) 
! 2 tC à 
Complex Numerical Distance p, . _ a "— ( 6 f ian sia (9) 
for Vertical Polarization i 2. sin (6j) 2 | n* 
Complex Numerical Distance Pm. - i mss ES | cos( ) + ES (sin Ow Ji 
for Horizontal Polarization ID ay i 
(n?.cos(6,)) - ( n* - sin(6}) 2) 
Vertical Reflection Coefficient Ivy SSS eng 


(n?. cos (4) | " n? - sin (6,)2) 


2 , 2 
Horizontal Reflection Coefficient Thy := cos( 6j) y y - sin(6,) | 
cos(6,) + (h a "— 


es 


Radiation Patterns in Phi=Pi/2 Plane (Perpendicular to Dipole) 
Space Wave Radiation Pattern (Phi=Pv2) 





Max E-Field Intensity (Volts per meter) max( MagE!P) = 0.03799 
Surface Wave Radiation Pattern (Phi=pi/2) 





Hel, } H2P, 


Max E-Field Intensity (Volts per meter)  max( MagE2P) - 6.21706: 10° 
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Radiation Patterns in Phi=0 Plane (Parallel to Dipole) 
Space Wave Radiation Pattern (Phi=0) 





Max E-Field Intensity (Volts per meter) max( MagE10) - 0.03422 
Surface Wave Radiation Pattern (Phi=0) 





Max E-Field Intensity (Volts per meter) — max( MagE20) — 9.38091-10 ^ 
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Radiation Patterns in Horizontal Plane 


Space Wave Radiation Pattern 





Max E-Field Intensity (Volts per meter) max( MagE1H) - 0.02972 


Surface Wave Radiation Pattern (Horizontal Plane) 


Phi-0 





net, : 428, 


Max E-Field Intensity (Volts per meter) max(MagE2H) = 1.14409- 106 
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Total Power Radiated (Watts) power = 31.69054 


Radiation Resistance (versus 
Maximum Antenna Current) 
(Ohms) 


Radres 


63.38107 


Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP = 6.8304 DirectivityO = 5.54301 
Effective Isotropic Radialed Power (EIRP) (Watts) 
Phi = pi/2 Plane Phi = 0 Plane 


DirecivilyP-power — 216.45893 DirectivilyO: power — 175.66087 


Maximum Effective Area (Along Radial of Directivity) 
(square melers) 


Phi = pi/2 Plane Phi = 0 Plane 
À nune 2 ec 
Ae) ^ DirectivityP Ae ) ^ DirectivityO 
Me a = 54.35457 Ep x = 44.10985 


Maximum Effective Length (Along Radial of Directivity) 
(meters) 


Phi = pi/2 Plane Phi = O Plane 


Radres: ( As) 2 DirectivityP Radres- ( As) à DirectivityO 





2- = 6.04592 2- = 5.44643 
480 X 2 480 X 2 
Numerical Distances 
Vertical Polarization Horizontal Polarization 
| Peg] = 26.28165 | Pmg| = 3.21687- iof 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi = pi/2 Plane Phi = 0 Plane 
18.542 89.85737 
AngleP = AngleO = 
18.542 89.85737 


HORIZONTAL DIPOLE 


This application calculates far field radiation patterns and parameters associated with 
horizontal thin-wire dipole antennas (diameter << wavelength). The antenna is 
mounted above and parallel to the x-axis in a rectangular coordinate system. The 
feed is at the center of the antenna at a set height directly above the origin. 

Required inputs are the antenna length, feed height above the surface, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. Predicted operating frequenctes assume that the antenna is a 
quarter- wavelength, half-wavelength, three-quarter-wavelength, or full wave- 
length dipole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. The radiation patterns are plotted for the phi=0 and phi-pi/2 
vertical planes parallel and perpendicular to the x-axis respectively. A horizontal 
radiation pattern is plotted at an elevation selected by the index from the elevation 
angle index table. Polarization is a combination of vertical and horizontal 
depending upon spatial orientation relative to the dipole. 


Input the dipole "E D 
length in meters 


No | — 


Input the index of the elevation angle for 
which to calculate the horizontal radiation d = 535 
pattern (from the angle index table) 


Wavelengths and Frequencies 


4-1 
^ jm 4-1 Az = 2.1 Ag se 3 As = | 
c c c c 
f = — f ‘=o — f = f Jd e 
1 à 3 4 
A D? ^3 M 
Possible Operating Frequencies (Hertz) 
KEN 7 Age | 7 
Dipole: f = 2.5.10 5 Dipole: [5 - 5:10 
s Dipole: f} = 7.5.10" A Dipol: {4 = 1-105 


Input the operating - 89.198 -— € 2 
Frequency (Hertz) D 9 fs s As 
Input the Height of Ha - 15 Input the ground - 7 
Antenna Feed (meters) 0 — Dielectric Constant 

Input the Distance R - 3000 _ Input the ground E 4 

from Antenna (meters) Conductivity: 


Rd; =R- Hg'cos(0;) 
Rr, = R+ Hp cos( 6i) 


Complex Numerical Distance 


P 
for Vertical Polarization i 


Complex Numerical Distance 
for Horizontal Polarization 


Vertical Reflection Coefficient 


Horizontal Reflection Coefficient 


Index of Refraction 








—— 
Bi eds i cos (6i) * OS i sin (9) 
2: sin (8,) | n? 
- j B: Rr | — 12 
E cos( 6) P |»? = (sin (6,) 2) 
sin ( 6j 


- eoo) - Ue - sinc?) 


I'v, = Ac 


(02.cos(a)) + (fa? - inca)? 


rh - cos ( 44) : h? - (ny?) 


cos(6,) 4 IE - "o 


Radiation Patterns in Phi=Pi/2 Plane (Perpendicular to Dipole) 
Space Wave Radiation Pattern (Phi=Pi/2) 





Max E-Field Intensity (Volts per meter) max ( MagE!IP) = 0.03926 
Surface Wave Radiation Pattern (Phi=pi/2) 





-] -0.5 0 0.9 ] 
net, HeP. 


Max E-Field Intensity (Volts per meter) = max( MagE2P) = 8.83121: 10? 
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Combined Space and Surface Wave Radiation Pattern (Phi-pi/2) 





Heb, A H3E, 


Max E-Field Intensity (Volts per meter) max( MagE3P) - 0.03926 
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Radiation Patterns in Phi-O Plane (Parallel to Dipole) 


Space Wave Radiation Pattern (Phi=0) 





Heo HIG, 


Max E-Field Intensity (Volts per meter) max( MagE10) = 0.03925 
Surface Wave Radiation Pattern (Phi=0) 
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neb, ) 20, 


Max E-Field Intensity (Volts per meter) max(MagE20) = 3.2224: 10* 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 








w30, 


net, 


Max E-Field Intensity (Volts per meter) 


max( MagE30) - 0.03925 
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Radiation Patterns in Horizontal Plane 


Space Wave Radiation Pattern 


Phi=0 





net, : alf 


Max E-Field Intensity (Volts per meter)  max( MagE1H) - 0.02647 


Surface Wave Radiation Pattern (Horizontal Plane) 


Phi-0 





Max E-Field Intensity (Volts per meter) max( MagE2H) = 4.07737: TS 
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Combined Space and Surface Wave Radiation Pattern (Horizontal Plane) 


Phi=0 





Heb, «3B. 


Max E-Field Intensity (Volts per meter)  max( MagE3H) - 0.02647 


2 1 


Total Power Radiated (Watts) power = 41.98012 


Radiation Resistance (versus 
Maximum Antenna Current) 
(Ohms) 


Radres 


83.96025 


Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP — 5.50819 DirectivityO = 5.50321 


Effective Isotropic Radiated Power (EIRP) (Watts) 


Phi = pi/2 Plane Phi = 0 Plane 


DirectivilyP: power = 231.23457 DirectivityO:- power = 231.02563 


Maximum Effective Area (Along Radial of Directivity) 
(square meters) 


Phi = pi/2 Plane Phi = 0 Plane 
(Ag) *: DirectivityP (^5) Directivityo 
Eu 0 457780 = 1876555 
. . A 


Maximum Effective Length (Along Radial of Directivity) 
(meters) 


Phi = pi/2 Plane Phi = 0 Plane 





Radres: (Ag) *: DirectivityP 





2: = 3.74932 2- = 3.74762 
480. x * 
Numerical Distances 
Vertical Polarization Horizontal Polarization 
| Peg| = 1.08943 | Pmg| = 2.26469- 109 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi = pi/2 Plane Phi = 0 Plane 
78.73217 89.85737 
AngleP = AngleO = ( 
78.73217 89.85737 


APPENDIX D: 
VERTICAL LOG-PERIODIC DIPOLE ARRAY COMPUTER OUTPUT 


This appendix contains computer hardcopies from the 
Mathcad vertical log-periodic dipole array application which 
show the input values and predicted radiation characteristics 
for two sample calculations. The configuration is a given by 
the inputs on the first two pages of each printout. The first 
antenna is mounted above soil (e,-10 and c-10?) and the second 
above seawater (ec,272 and 0-4). Reference 6 [p. 114] provides 
the radiation patterns and gain predictions from several 
sources for the configuration in the first example. 

The radiation patterns and maximum directive gain 
computed by the first Mathcad example are almost identical to 
the those given in reference 6. The directivity prediction is 
slightly higher for Mathcad, and the elevation of the maximum 
directive gain is also slightly higher, but the overall 
Similarity between the predicted radiation characteristics 1s 
Still quite good. The seawater example also yields results 
consistent with expectations. With respect to the soil 
example, the seawater example’s directivity is higher due to 
a stronger reflected wave, and the surface wave is stronger at 


grazing angles (090?) due to higher conductivity. 
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VERTICAL LOG PERIODIC DIPOLE ARRAY 


This application calculates far field radiation patterns and parameters associated with 
vertical log periodic dipole arrays. The antenna is oriented such that the projection 
of its center axis lies on the positive y-axis of a rectangular coordinate system. The 
antenna axis can be oriented with respect to the vertical at any angle between zero 
and ninety degrees. The feed is at the center of the shortest element directly above 
the origin. The dipole elements are bisected by the antenna axis and are parallel to 
the z-axis. Required inputs are the number of dipole elements, shortest and second 
shortest element lengths, separation between the shortest and second shortest 
elements, radius of the shortest element, height of the feed above the surface, 

angle of the antenna axis with the vertical, characteristic admittance of the line 
feeding the antenna, termination impedance of the antenna, transmitted frequency, 
distance from the antenna, the conductivity and dielectric constant of the surface 
below the antenna. The planar earth model is assumed in predicting radiation 
patterns. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. Radiation patterns are plotted for the phi-0 and phi-pi/2 vertical 
planes perpendicular and parallel to the y-axis respectively. A horizontal radiation 
pattern is plotted at an elevation selected by the index from the elevation angle 
index table. Polarization is vertical for all vertical log-periodic antennas. 


Input the number 
of elements 


Input the length 
of the shortest 
element (meters) 


Input the distance 
from shortest to 
second shortest 
element (meters) 


Input the operating 
Frequency (Hertz) 


Input the ground 
Conductivity 


Input the Characteristic 
Admittance of the Line ADM = 
Feeding the Antenna 


(Ohms) 


N 


[5 


O - 


= 12 


= 3.073 


= 1.0766 


= 18-10° 


1:107 3 


Es 
450 


Input the shortest 
element's Radius 
(meters) 


Input the length 
of the second 
shortest element 
(meters) 


Input the height of 
the shortest element 
above the surface 
(meters) 


Input the Distance 
from the Antenna 
(meters) 


Input the ground 
Dielectric Constant 


Input the Termination 


radg - .00918 
ly = 4.373 
Hg `= 6.87 

R = 3000 

er = 10 


Impedance connected TIMP = 450 


to longest element 
(Ohms) 
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Input the angle between Input the elevation index 
the vertical (z axis) and y - 78 for which to calculate the w - 285 
the antenna axis (degrees) horizontal radiation pattern 

(from angle index table) 


NOTE: The Vertical Log-Periodic Dipole mutual impedance calculations are nol 

valid if there is an antenna element whose length is an exact integer multiple of the 
wavelength. If this occurs, there will be a singularity error in the mutual impedance 
calculations. If this problem arises, it will be necessary to vary the operating frequency 
such that no element is exactly an integer multiple of the wavelength. 


If the log-periodic calculations of element length and spacing do not represent the 
desired antenna configuration, you can enter the values directly by selecting the 
variables | and d, using the define key (shift, colon), and entering the lengths and 
spacings by separating each successive entry by a comma. 


Calculated distance 


Calculated length between successive Calculated radii 
of elements from elements-shortest of elements from 
shortest to longest to longest shortest to longest 


rad 





eo 


Radiation Patterns in Phi=P1/2 Plane 
(Perpendicular to X -Axis) 


Space Wave Radiation Pattern (Phi=Pi/2) 





=| -0.5 0 0.5 1 
nre , HIP, 
Max E-Field Intensity (Volts per meter) max( MagE1P) - 0.05599 


Surface Wave Radiation Pattern (Phi-pi/2) 





El -0.5 


- 0 0.5 1 
nref, Mf 


Max E-Field Intensity (Volts per meter) max( MagE2P) = 0.00124 


eu 


Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 





nef, RF; 


Max E-Field Intensity (Volts per meter) max( MagE3P) - 0.05583 
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Radiation Patterns in Phi=O Plane 
(Perpendicular to Phi=pi/2 Plane) 


Space Wave Radiation Pattern (Phi=0) 





-1 -0.5 0 0.5 1 
ure, x10, 


Max E-Field Intensity (Volts per meter) max( MagE!0) = 0.03221 
Surface Wave Radiation Pattern (Phi=0) 
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-1 -0.5 0 0.5 1 
Reb, : #20, 


Max E-Field Intensity (Volts per meter) ^ max( MagE20) - 7.64667. ro * 


ZI 


Combined Space and Surface Wave Radiation Pattern (Phi=0) 





Reb, : 430. 


Max E-Field Intensity (Volts per meter) — max( MagE30) - 0.0321 


Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 


Space Wave Radiation Pattern (Horizontal) 
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-1 -0.5 0 0.5 1 
net , «tH, 


Max E-Field Intensity (Volts per meter) max( MagE1H) - 0.05438 


Surface Wave Radiation Pattern (Horizontal) 


Phi-0 





Max E-Field Intensity (Volts per meter) — max( MagE2H) - 1.09861: 10 $ 
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Combined Space and Surface Wave Radiation Pattern (Horizontal) 


Phi=0 





Max E-Field Intensity (Volts per meter) max( MagE3H) - 0.05427 


ae 





ane 30 y pas + a3)-(cos(5)-sm(al + 103) + sm(Z) 
2 Y 
TR 
X | Y 
-— 40 
J 2 
power = poweré Radres - 2: power 
2 
MagE 3P, us. 
aram . VM) pep ~- £7 Rmax ( sqE3P) 
2-(120: 7) power 
(| MagE30,|) * 4-1: R2. max ( sqE30 ) 
sqE30, = ————— — —— DirectivityO = —————————————— 
| 2:(120- x) power 
Total Power Radiated (Watts) power = 34.29572 


Radiation Resistance (versus 
Maximum [nput Current) 
(Ohms) 


Radres = 68.59144 


Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP = 13.63346 DirectivityO = 4.50794 
Effective Isotropic Radiated Power (EIRP) (Watts) - 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP- power = 467.56917 DirectivityO- power = 154.6031 


Maximum Effective Area (Along Radial of Directivity) 
(square meters) 


Phi = pi/2 Plane Phi = 0 Plane 
ON IO A E. 
Ac): DirectivityP Az ) ^ DirectivityO 
CS ae = 301.36553 Ma = 99.64739 
X E 


Maximum Effective Length (Along Radial of Directivity) 





(meters) 
Phi = pi/2 Plane Phi 7 0 Plane 
Radres- ( As) 2. DirectivilyP Radres (As) E DirectivityO 
2: [——————————————— = 14.8097 2: |——————————————— = 8.515394 
480. x* 480. * 


Numerical Distance for Vertical Polarization 


| Peg | = 50.7 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi = py2 Plane Phi = O Plane 
13.71429 
AngleP = 14.28571 Angle0 = 
13.71429 
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VERTICAL LOG PERIODIC DIPOLE ARRAY 


This application calculates far field radiation patterns and parameters associated with 
vertical log periodic dipole arrays. The antenna is oriented such that the projection 
of its center axis lies on the positive y-axis of a rectangular coordinate system. The 
antenna axis can be oriented with respect to the vertical at any angle between zero 
and ninety degrees. The feed is at the center of the shortest element directly above 
the origin. The dipole elements are bisected by the antenna axis and are parallel to 
the z-axis. Required inputs are the number of dipole elements, shortest and second 
shortest element lengths, separation between the shortest and second shortest 
elements, radius of the shortest element, height of the feed above the surface, 

angle of the antenna axis with the vertical, characteristic admittance of the line 
feeding the antenna, termination impedance of the antenna, transmitted frequency, 
distance from the antenna, the conductivity and dielectric constant of the surface 
below the antenna. The planar earth model is assumed in predicting radiation 
patterns. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. Radiation patterns are plotted for the phi-0 and phi=pi/2 vertical 
planes perpendicular and parallel to the y-axis respectively. A horizontal radiation 
pattern is plotted at an elevation selected by the index from the elevation angle 
index table. Polarization is vertical for all vertical log-periodic antennas. 


Input the number N - 12 Input the shortest radg :- .00918 
of elements element's Radius 

(meters) 
Input the length Input the length 
of the shortest lg * 3.673 of the second I, = 4.373 
element (meters) shortest element 

(meters) 
Input the distance Input the height of 
from shortest to dg = 1.0766 the shortest element Hg = 6.87 
second shortest above the surface 
element (meters) (meters) 
Input the operating — ( _ ¡8.199 Input the Distance R - 3000 
Frequency (Hertz) S from the Antenna 

(meters) 
Input the ground Bu Input the ground a 7 
Conductivity Dielectric Constant 
Input the Characteristic 1 Input the Termination 
Admittance of the Line ADM = 450 [mpedance connected TIMP = 450 
Feeding the Antenna to longest element 
(Ohms) (Ohms) 
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Input the angle between Input the elevation index 
the vertical (z axis) and V - 78 for which to calculate the w - 285 
the antenna axis (degrees) horizontal radiation pattern 

(from angle index table) 


NOTE: The Vertical Log-Periodic Dipole mutual impedance calculations are not 

valid if there is an antenna element whose length is an exact integer multiple of the 
wavelength. If this occurs, there will be a singularity error in the mutual impedance 
calculations. If this problem arises, it will be necessary to vary the operating frequency 
such that no element is exactly an integer multiple of the wavelength. 


If the log-periodic calculations of element length and spacing do not represent the 
desired antenna configuration, you can enter the values directly by selecting the 
variables ] and d, using the define key (shift, colon), and entering the lengths and 
spacings by separating each successive entry by a comma. 


Calculated distance 
Calculated length between successive Calculated radii 
of elements from elements-shortest of elements from 
shortest to longest to longest shortest to longest 
rad 
Y 





Radiation Patterns in Phi=P1/2 Plane 
(Perpendicular to X -Axis) 


Space Wave Radiation Pattern (Phi=Pi/2) 





-1 -0.5 0 
utet, HIP. 
Max E-Field Intensity (Volts per meter) 


0.5 1 


max( MagEIP) - 0.10756 
Surface Wave Radiation Pattern (Phi=pi/2) 


raro rr 


-05 [X 





-1 


Max E-Field Intensity (Volts per meter) max( MagE2P) - 0.12392 


Cai 


Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 
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max( MagE3P) - 0.12392 


Max E-Field Intensity (Volts per meter) 


Radiation Patterns in Phi=O Plane 
(Perpendicular to Phi=pi/2 Plane) 


Space Wave Radiation Pattern (Phi=0) 


yret, 





nel, , x10. 


Max E-Field Intensity (Volts per meter) max( MagE10) - 0.06687 
Surface Wave Radiation Pattern (Phi=0) 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 
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max ( MagE30) = 0.07827 


Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 


Space Wave Radiation Pattern (Horizontal) 


Phi=0 
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net, HIE, 


Max E-Field Intensity (Volts per meter) max( MagE1H) - 0.06821 


Surface Wave Radiation Pattern (Horizontal) 
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-0.3 


0.5 
Heb 42H. 
Max E-Field Intensity (Volts per meter) 


max( MagE2H) = 4.93203- 10 
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Combined Space and Surface Wave Radiation Pattern (Horizontal) 


Phi=0 





Heh, i H3H. 


Max E-Field Intensity (Volts per meter) max( MagE3H) - 0.0682 


eee 





X aen 
2 |2 
j-p-Y -sec(a2 + a3)-(cos(C)-sm(a2 + a3) + smt) 
power - 2e y e ' 
2 Y 
xR 
1 Y 
a 
4 2 
power := power Radres =- 2- power 
2 
MagE3P. a- RÊ. 
sqE3P. = qa)" DrrectivityP = 4: z: R max( sgE3P) 
| 2. (120- x) power 
2 
MagE30, a RÊ 
Bii, E Directivityo - S RÀ max (sqE30) 
2. (120: x) power 
Total Power Radiated (Watts) power = 99.48703 


Radiation Resistance (versus 
Maximum Input Current) 
(Ohms) 


Radres = 198.97405 


Directivity (or Maximum Power Gain assuming 10076 Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP = 23.15454 DirectivityO = 9.2377 
Effective Isotropic Radiated Power (EIRP) (Watts) | 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP: power = 2.30358: 10° DirectivityO: power = 919.03096 


Maximum Effective Area (Along Radial of Directivity) 
(square meters) 


Phi = p1/2 Plane Phi = 0 Plane 
2 r oe i à . f 
Ae ) ^: DirectivityP Ae): DirectivityO 
Cs ma = 511.82763 Sl ee = 204.19793 
R A 


Maximum Effective Length (Along Radial of Directivity) 
(meters) 


Phi 7 pi/2 Plane Phi 7 0 Plane 


Radres: ( As) e DirectivityP Radres- (As) A DirectivityO 


= 32.87188 2. = 20.76291 





480-27 480-22 


Numerical Distance for Vertical Polarization 


| Peg| = 0.14135 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi = pi/2 Plane Phi = 0 Plane 


0 
AngleP = 0 AngleO = | | 
0 


APPENDIX E: 
HORIZONTAL LOG-PERIODIC DIPOLE ARRAY COMPUTER OUTPUT 


This appendix contains computer hardcopies from the 
Mathcad horizontal log-periodic dipole array application which 
show the input values and predicted radiation characteristics 
for two sample calculations. The configuration is a given by 
the inputs on the first two pages of each printout. The first 
antenna is mounted above soil (e,-4 and 0-107) and the second 
above seawater (c,-72 and 0-4). Reference 6 [p. 108] provides 
the radiation patterns and gain predictions from several 
sources for the configuration in the first example. 

The radiation patterns and maximum directive gain 
computed by the first Mathcad example are almost identical to 
the those given in reference 6. The directivity prediction is 
slightly higher for Mathcad, but the overall similarity among 
the predicted radiation patterns is remarkably good. The 
seawater example has results very much like the soil example. 
For horizontal polarization, the higher conductivity surface 
below the antenna does not result in a greatly enhanced 
surface wave and increased directivity as it does for vertical 


polarization. 
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HORIZONTAL LOG PERIODIC DIPOLE ARRAY 


This application calculates far field radiation patterns and parameters associated with 
horizontal log periodic dipole arrays. The antenna is oriented such that the 
projection of its center axis lies on the positive y-axis of a rectangular coordinate 
system. The antenna axis can be oriented with respect to the vertical at any angle 
between zero and ninety degrees. The feed is at the center of the shortest element 
directly above the origin. The dipole elements are bisected by the antenna axis and 
are parallel to the x-axis. Required inputs are the number of dipole elements, 
shortest and second shortest element lengths, separation between the shortest and 
second shortest elements, radius of the shortest element, height of the feed above 
the surface, angle of the antenna axis with the vertical, characteristic admittance of 
the line feeding the antenna, termination impedance of the antenna, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. À sinusoidal current input with a maximum of unity is assumed . 
All radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. Radiation patterns are plotted for the phi=pi/2 and phi=0 vertical 
planes parallel and perpendicular to the y-axis respectively. A horizontal radiation 
pattern is plotted at an elevation selected by the index from the elevation angle 
index table. Polarization is a varying combination of horizontal and vertical 
depending on spatial orientation with respect to the antenna. 


Input the number N - 12 
of elements 

Input the length 

of the shortest lg ^ 
element (meters) 


7.249 


Input the distance 
from shortest to 
second shortest 
element (meters) 


dg - 1.659 


Input the operating | 05. 108 
Frequency (Hertz) is E 
Input the ground 


2 4.1973 
Conductivity at ae 


Input the Characteristic 

Admittance of the Line i 
Feeding the Antenna “DM = 300 
(Ohms) 


Input the shortest 


element's Radius rado = -00144 
(meters) 

Input the length 

of the second l ="8:297 


shortest element 
(meters) 


Input the height of 

the shortest element H g^ 94 
above the surface 

(meters) 


Input the Distance 


from the Antenna R = 3000 
(meters) 
Input the ground - 4 


Dielectric Constant 


Input the Termination 
Impedance connected 
to longest element 
(Ohms) 


TIMP - 300 


EID 


Input the angle between Elevation index for which 


the vertical (z axis) and Y - 90 to calculate horizontal pattern W = 285 
the antenna axis (degrees) (from angle index table) 


NOTE: The Horizontal Log-Periodic Dipole mutual impedance calculations are not 
valid if there is an antenna element whose length is an exact integer multiple of the 
wavelength. If this occurs, there will be a singularity error in the mutual impedance 
calculations. If this problem arises, it will be necessary to vary the operating frequency 
such that no element is exactly an integer multiple of the wavelength. 


If the log-periodic calculations of element length and spacing do not represent the 
desired antenna configuration, you can enter the values directly by selecting the 
variables | and d, using the define key (shift, colon), and entering the successive 
lengths and spacings by separating each successive entry by a comma. 


Calculated distance 


Calculated length between successive Calculated radii 
of elements from elements-shortest of elements from 
shortest to longest to longest shortest to longest 





Radiation Patterns in Phi=Pi/2 Plane 
(Perpendicular to X-Axis and Dipole Elements) 
Space Wave Radiation Pattern (Phi=Pi/2) 





-0.5 


0.5 1 
Het. HIP, 
Max E-Field Intensity (Volts per meter) 


max ( MagF iP) 


= 0.08512 
Surface Wave Radiation Pattern (Phi=pi/2) 
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Reb. x2P. 
Max E-Field Intensity (Volts per meter) 


max( MagE2P) = 4.54242- 10? 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 





neb, , "3P, 


Max E-Field Intensity (Volts per meter) max( MagE3P) - 0.08511 
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Radiation Patterns in Phi-0 Plane 
(Perpendicular to Phi=pi/2 Plane) 


Space Wave Radiation Pattern (Phi-O) 
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-1 -0.5 0 0.5 1 
reto. #10. 


Max E-Field Intensity (Volts per meter) max( MagE10) - 0.03495 
Surface Wave Radiation Pattern (Phi=0) 
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utet, «20, 


Max E-Field Intensity (Volts per meter)  max( MagE20) - 1.45261: (ole 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 
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Max E-Field Intensity (Volts per meter)  max( MagE30) - 0.03495 
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Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 


Space Wave Radiation Pattern (Horizontal) 


Phi-0 





-1 -05 0 0.5 1 
utet, . ATH 


Max E-Field Intensity (Volts per meter) max( MagE1H) - 0.0676 
Surface Wave Radiation Pattern (Horizontal) 


Phi=0 





ure, 1 42H. 


Max E-Field Intensity (Volts per meter) — max( MagE2H) — 3.44012: 10 
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Combined Space and Surface Wave Radiation Pattern (Horizontal) 


Phi=0 





Max E-Field Intensity (Volts per meter) max ( MagE3H) = 0.06758 
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1.5707963 1.5707963 
COs econ Poca 








power? = AC eee 
x: R? y m (sin( £):cos( £) )* 
- - 1.5707963 -0 
E 
2 2 i 
0030 SE | -P-Y, -ese(') C eos( 7 ) cost 
poweró - ———MÀÀ— NM 
xR? >| 1 - (sin(£):cos( £))? 
a 
hit 
+ 2 
power = poweré + power¢ Radres = 2- power 
(| MagE3P; |) $ 4: x R^ max( sqE3P) 
SQE3P, s — ee DirectivityP ::. ————————— ——— 
2: (120m ) power 
MagE30 a Re. 
sq 30, ‘= (MED. Directivit yO = 4- z: R max( sqE30) 
2: (120: x) power 
Total Power Radiated (Watts) power = 83.98604 
Radiation Resistance (versus Radres = 167.97209 
Maximum Input Current) 
(Ohms) 


Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi 7 pV/2 Plane Phi = 0 Plane 


DirectivityP = 12.93749 DirectivityO = 2.18175 
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a EULUPTPPÍ- IL PL» LL LLLIPLLLUL'LLLV]OCLUL,EOGO 


Effective Isotropic Radiated Power (EIRP) (Watts) 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP- power = 1.08657- 10° DirectivityO: power = 183.23657 


Maximum Effective Area (Along Radial of Directivity) 
(square meters) 


Phi = pi/2 Plane Phi = 0 Plane 
Qu. e Peg cas 
Ag)” DirectivityP Ag)” DirectivityO 
(Ns) "Direcüviy _ 643.45795 IS On 1096112: 

4-1 4-1 

Maximum Effective Length (Along Radial of Directivity) 

(meters) 
Phi = pi/2 Plane Phi = O Plane 


Radres: (As) a DirectivityP 





2 
480. x? 
Numerical Distances 
Vertical Polarization Horizontal Polarization 
| Peg| = 69.28584 | Pmg| 7 1.26447- 10? 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi = pi/2 Plane Phi = 0 Plane 
89.71429 
AngleP = 31.42857 Angle0 = 
89.71429 
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HORIZONTAL LOG PERIODIC DIPOLE ARRAY 


This application calculates far field radiation patterns and parameters associated with 
horizontal log periodic dipole arrays. The antenna is oriented such that the 
projection of its center axis lies on the positive y-axis of a rectangular coordinate 
system. The antenna axis can be oriented with respect to the vertical at any angle 
between zero and ninety degrees. The feed is at the center of the shortest element 
directly above the origin. The dipole elements are bisected by the antenna axis and 
are parallel to the x-axis. Required inputs are the number of dipole elements, 
shortest and second shortest element lengths, separation between the shortest and 
second shortest elements, radius of the shortest element, height of the feed above 
the surface, angle of the antenna axis with the vertical, characteristic admittance of 
the line feeding the antenna, termination impedance of the antenna, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radialion patterns. A sinusoidal current input with a maximum of unity is assumed . 
All radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. Radiation patterns are plotted for the phi=pi/2 and phi=0 vertical 
planes parallel and perpendicular to the y-axis respectively. A horizontal radiation 
pattern is plotted at an elevation selected by the index from the elevation angle 
index fable. Polarization is a varying combination of horizontal and vertical 
depending on spatial orientation with respect to the antenna. 


Input the number N = 12 
of elements 


Input the length 
of the shortest Io = 
element (meters) 


7.219 


Input the distance 
from shortest to 
second shortest 
element (meters) 


‘= 1.659 


Input the operating 
Frequency (Hertz) 5 


Input the ground 
Conductivity 


Input the Characteristic 
Admittance of the Line 
Feeding the Antenna 
(Ohms) 


ADM - 


g|- 


Input the shortest 
element's Radius 
(meters) 


rado = .00144 


Input the length 
of the second 
shortest element 
(melers) 


Input the height of 

the shortest element H 0^ 
above the surface 

(meters) 


9.4 


Input the Distance 
from the Antenna 
(meters) 


R - 3000 


Input the ground 
Dielectric Constant 


Input the Termination 
Impedance connected 
to longest element 
(Ohms) 


TIMP - 300 
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Input the angle between Elevation index for which 
the vertical (z axis) and V - 90 tocalculate horizontal pattern W = 289 


the antenna axis (degrees) (from angle index table) 


NOTE: The Horizontal Log-Periodic Dipole mutual impedance calculations are not 
valid if there is an antenna element whose length is an exact integer multiple of the 
wavelength. If this occurs, there will be a singularity error in the mutual impedance 
calculations. If this problem arises, it will be necessary to vary the operating frequency 
such that no element is exacily an integer multiple of the wavelength. 


If the log-periodic calculations of element length and spacing do not represent the 
desired antenna configuration, you can enter the values directly by selecting the 
variables | and d, using the define key (shift, colon), and entering the successive 
lengths and spacings by separating each successive entry by a comma. 


Calculated distance 


Calculated length between successive Calculated radii 

of elements from elements-shortest of elements from 

shortest to longest to longest shortest to longest 
! d rad 


V 
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Radiation Patterns in Phi=Pi/2 Plane 
(Perpendicular to X-Axis and Dipole Elements) 
Space Wave Radiation Pattern (Phi=Pi/2) 





Heb HIP, 


Max E-Field Intensity (Volts per meter) max ( MagE1P) = 0.13724 


Surface Wave Radiation Pattern (Phi=pi/2) 





net, : R2P. 


Max E-Field Intensity (Volts per meter) max( MagE2P) - 3.19958. 10 ^ 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 





neb, : HIP, 


Max E-Field Intensity (Volts per meter) max( MagE3P) - 0.13724 


Radiation Patterns in Phi-0 Plane 
(Perpendicular to Phi-pi/2 Plane) 


Space Wave Radiation Pattern (Phi=0) 





-1 -0.5 0 
nel , #10; 


Max E-Field Intensity (Volts per meter) max( MagE10) - 0.05843 


0.5 1 


Surface Wave Radiation Pattern (Phi=0) 
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-1 -0.5 0 
utet, : #20; 


Max E-Field Intensity (Volts per meter) 


max( MagE20) - 4.46943. 10 * 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 





Kieb, i #30, 


Max E-Field Intensity (Volts per meter)  max( MagE30) - 0.05843 
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Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 


Space Wave Radiation Pattern (Horizontal) 





fme 
p 
o5 p 
n j 
yet, : la 
tee 0 p | - Phi=0 
PS ME 
EM X uu t 
-05 [-. 
-1 -0.5 0 0.5 1 


Reb, ; HIE, 


Max E-Field Intensity (Volts per meter) max( MagE1H) = 0.09414 
Surface Wave Radiation Pattern (Horizontal) 


Phi-0 





Het, 2H, 


Max E-Field Intensity (Volts per meter) |= max( MagE2H) = 5.57296- (DES 


Ae 


Combined Space and Surface Wave Radiation Pattern (Horizontal) 


Phi-0 





nre, à 43H, 


Max E-Field Intensity (Volts per meter) max ( MagE3H) = 0.09414 


29:3 


1.5707963 1.5707963 


cos( £ ):cos( S hoec 
"e 





xR? y 1- (sin(£)-cos(£))* 
| --1.5707963 -'0 


sin( £)-1 Jg: Y -esc( P )-.( cos(6 )-cos( 4 
poweró - y 2r ol — ONT 





xR? y | 1 - Cin(£):cos(£))* 
x 
-> -0 
/ 2 
power = power + power¢ Radres = 2-power 
2 
MagE3P q RE. 
sqE3P; 'z (Map DirectivityP = 4: x: R max( sqE3P) 
2°( 120-2) power 
2 
MagE30, ARE 
qe3o > IAB irectivityo LT RÓmax(59E30) 
2-(120-1x) power 
Total Power Radiated (Watts) power = 261.928 


Radiation Resistance (versus 
Maximum Input Current) 
(Ohms) 


Radres = 523.856 


Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 


DirectivityP = 10.78579 DirectivityO = 1.95488 
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Effective Isotropic Radiated Power (EIRP) (Watts) 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP: power = 2.8251 - 10° DirectivityO: power = 512.03765 


Maximum Effective Area (Along Radial of Directivity) 
(square meters) 





Phi = pi/2 Plane Phi = 0 Plane 
A« ) * DirectivityP Ac) *: Directivit 

(As) ivityP (As) irectivityO 

—————————— z $36.4412 —————————— 29722772 
4. x 4. x 
Maximum Effective Length (Along Radial of Directivity) 

(meters) 
Phi = pi/2 Plane Phi = 0 Plane 
Radres: (As) a DirectivityP Radres: (As) 2 DirectivityO 
2: [——————————————— = 54.60488 2. |[——— ——— — ——————— —— z 23.24692 

480: x? 480. x^ 


Numerical Distances 
Vertical Polarization Horizontal Polarization 


| Peg| = 0.06283 | Pmg| = 2.26211-10° 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi = pi/2 Plane Phi = 0 Plane 
6.5714 
AngleP = 34.85714 Angled = E à 
56.57143 
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APPENDIX F: 
HORIZONTAL YAGI-UDA ARRAY COMPUTER OUTPUT 


This appendix contains computer hardcopies from the 
Mathcad horizontal Yagi-Uda array application which show the 
input values and predicted radiation characteristics for two 
sample calculations. The configuration is a given by the 
inputs on the first three pages of each printout. The first 
antenna is mounted above soil (e,24 and 0-10?) and the second 
above seawater (ec,272 and 0-4). Reference 6 [p. 107] provides 
the radiation patterns and gain predictions from several 
sources for the configuration in the first example. 

The radiation patterns and maximum directive gain 
computed by the first Mathcad example are almost identical to 
the those given in reference 6. The directivity prediction is 
slightly higher for Mathcad, but the overall similarity among 
the predicted radiation patterns is very good. The seawater 
example has results very much like the soil example. For 
horizontal polarization, the higher conductivity surface below 
the antenna does not result in a greatly enhanced surface wave 
and increased  directivity as it does for vertical 


polarization. 
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HORIZONTAL YAGI-UDA ARRAY 


This application calculates far field radiation patterns and parameters associated with 
horizontal Yagi-Uda arrays. The antenna is oriented such that the projection of its 
center axis lies on the posilive y-axis of a rectangular coordinate system. The feed 
is at the center of a single driven element between a set number of longer parasitic 
reflector elements and shorter parasitic director elements. Antenna elements are 
bisected by the antenna axis and are parallel to the x-axis. The center of the first 
reflector is directly above the origin at a set height on the z-axis. Required inputs 
are the number of reflector and director elements, individual element lengths, 
individual element radii, separation between adjacent elements, height of the feed 
above the surface, transmitted frequency, distance from the antenna, conductivity 
and dielectric constant of the surface below the antenna. The planar earth model is 
assumed in predicting radiation patterns. À sinusoidal voltage with a one volt 
maximum across the input terminals is assumed. All radiation patterns are 
normalized with respect to the maximum electric field intensity transmitted by the 
antenna in the plane of interest. The electric field magnitudes to which the patterns 
are normalized are displayed below their respective plots. Radiation patterns are 
plotted for the phi=pi/2 and phi-O vertical planes parallel and perpendicular to the 
y-axis respectively. A horizontal radiation pattern is plotted at an elevation 
selected by the index from the elevation angle index table. Polarization is a 
varying combination of horizontal and vertical depending on spatial orientation with 
respect to the antenna. 


NOTE: The Yagi-Uda mutual impedance calculations are not valid if there is an 
antenna element whose length is an exact integer multiple of the wavelength. If 

this occurs, there will be a singularity error in the mutual impedance calculations. 
This should not be a problem since Yagi-Uda Antennas are designed to operate at 
frequencies such that the driven element is approximately one-half- wavelength long, 
the reflectors are slightly longer, and the directors are slightly shorter. 1f the problem 
does arise, it will be necessary to vary the operating frequency such that no element 
is exactly an integer multiple of the wavelength. 


Input the number of wp. y; Input the numberof yr ... , 
reflector elements director elements 
y := O.. (NR + ND) 0 := O..(NR + ND - 1) 


This entire page is alloted for these entries to allow for as many elements as necessary. 


Input the lengths and radii of the Yagi 
elements, starting with the outermost 
reflector and proceeding inward to the 
driven element. After entering the driven 
element length and radius, enter those of 
the director elements, starting with the 
director next to the driven element and 
proceeding outward to the final element. 


l = rad - 
Y Y 





‚001 | 


Input the distances between 
successiv elements in the same 
order as the length and radius 
inputs. There should be one less 
entry in the separation distance 
array than there is in the length 
and radius array. 


Input the operating fs =- 10: 10° Input the height of Hy = 7.5 
Frequency (Hertz) the array above the 
surface (meters) 


Input the ground o - 1-10 : Input the ground f - d 
Conductivity Dielectric Constant 

Input the Distance = 3000 Index of Refraction 
from the Antenna 

(meters) 





Input the index of the elevation 

angle for which to calculate the W = 285 
horizontal radiation pattern 

(from the angle index table) 


Complex Numerical Distance for Vertical Polarization 


2 
In? - sin (6) ^ 


Sip s + Hg cos(8,)). cos(0) " 


2 - sin (9;) 2 


Pe; = 


Complex Numerical Distance for Horizontal Polarization 


Pm, - PA A] cos 0; x n* - m 6)? : 
2-( sin (6;) 2) | ( ) ( ( ) )| 


(n2cos(0;)) = M" a sin ( 6;) 2 
Vertical Reflection Coefficient EA AAA AA 


© (1%cos(a)) + (fa? - inca)? 


on 2) 
Horizontal Reflection Coefficient Fh; = cos (6;) j TN sin (6) 


cos ( 6j) + (h? - — 


Radiation Patterns in Phi=Pi/2 Plane 
(Perpendicular to X-Axis and Dipole Elements) 
Space Wave Radiation Pattern (Phi-Pi/2) 





net, , HIP. 


Max E-Field Intensity. (Volts per meter) max ( MagE1P) = 5.13224: io f 


Surface Wave Radiation Pattern (Phi=pi/2) 





Max E-Field Intensity (Volts per meter) max( MagE2P) = 4.08957: TUR 
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Combined Space and Surface Wave Radiation Pattern (Phi=p1/2) 
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Heb, ' H3P, 


Max E-Field Intensity (Volts per meter) max( MagE3P) = 5.13132: (0 * 


IET 


Radiation Patterns in Phi-O Plane 
(Perpendicular to Phi=pi/2 Plane) 


Space Wave Radiation Pattern (Phi=0) 
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Max E-Field Intensity (Volts per meter) max(MagE10) = 1.35264- T 


Surface Wave Radiation Pattern (Phi=0) 
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-1 -0.5 0 0.5 1 
Refi i x20, 
Max E-Field Intensity (Volts per meter)  max( MagE20) = 5.32341- 10? 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 
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Heb, : 430. 


Max E-Field Intensity (Volts per meter) 
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max( MagE30) = 1.35264- 10 


4 


Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 


Space Wave Radiation Pattern (Horizontal) 


Phi=0 





neb, ! RIH, 


Max E-Field Intensity (Volts per meter) max( MagE1H) = 3.79643- Ww 


Surface Wave Radiation Pattern (Horizontal) 





-1 -0.5 0 05 1 
Heb «2H. 


Max E-Field Intensity (Volts per meter) — max( MagE2H) - 2.73697- 10 © 
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Combined Space and Surface Wave Radiation Pattern (Horizontal) 





-4 
Max E-Field Intensity (Volts per meter) max(MagE3H) = 3.7942-10 ” 


2/55 








1.5707963 
power - La 
nm 
| --1.5707963 -'0 
M 
2 2 
powerg = 
x-R? 
-3 Jo 
“ 2 


power =- power + power¢g 


^W MagE3P; |) * 


SES; = 120.1) 
(| MagE30, | ) 


Total Power Radiated (Watts) 


Radiation Resistance (versus 
Maximum Input Current) 
(Ohms) 


| - (sin(£):cos( £))? 


1.5707963 
cos(€)-cos(2)-L  jp-D smc: 
ae 
| - (sin(£):cos(£))* 
sin(£)«1 |. 9-D -si(C) si(£) 


—Ó (cosd 


(VREF yp)? 


2: power 


Radres - 


4: x- R^ max ( sqE3P) 
power 


DirectivityP :=- 


4. x- R*. max ( sqE30) 
power 


DirectivityO = 


power = 0.00266 


Radres = 187.94973 


Directivity (or Maximum Power Gain assuming 10096 Antenna Efficiency) 


Phi = pi/2 Plane 


DirectivityP = 14.84639 


Phi = O Plane 


DirectivityO = 1.03163 
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Effective Isotropic Radiated Power (EIRP) (Watts) 
Phi = pi/2 Plane Phi = 0 Plane 
DirectivityP- power — 0.0395 DirectivilyO: power = 0.00274 


Maximum Effective Area (Along Radial of Directivity) 
(square meters) 








Phi = pi/2 Plane Phi = 0 Plane 
(Ag) * DirectivityP 3... (s)? Direcivityo 
A—— —— —————— cz 1.06329. 10 ` = 73.88536 
4. X 4-1 
Maximum Effective Length (Along Radial of Directivity) 
(meters) 
Phi = pi/2 Plane Phi = 0 Plane 
Radres- (As) : DirectivityP Radres: (Ag) 2. Directivity0 
2. [————————————— 2 46.04814 2. |————————————— —— - 12.13848 
480. x? 480- x? 
Numerical Distances 
Vertical Polarization Horizontal Polarization 
pa n 3 
| Peg | = 57.12622 | Pmg | = 1.09911-10 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi = pi/2 Plane Phi = 0 Plane 
89.71429 
AngleP = 36 AngleO = | | 
89.71429 
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HORIZONTAL YAGI-UDA ARRAY 


This application calculates far field radiation patterns and parameters associated with 
horizontal Yagi-Uda arrays. The antenna is oriented such that the projection of its 
center axis lies on the positive y-axis of a rectangular coordinate system. The feed 
is at the center of a single driven element between a set number of longer parasitic 
reflector elements and shorter parasitic director elements. Antenna elements are 
bisected by the antenna axis and are parallel to the x-axis. The center of the first 
reflector is directly above the origin at a set height on the z-axis. Required inputs 
are the number of reflector and director elements, individual element lengths, 
individual element radii, separation between adjacent elements, height of the feed 
above the surface, transmitted frequency, distance from the antenna, conductivity 
and dielectric constant of the surface below the antenna. The planar earth model is 
assumed in predicting radiation patterns. A sinusoidal voltage with a one volt 
maximum across the input terminals is assumed. All radiation patterns are 
normalized with respect to the maximum electric field intensity transmitted by the 
antenna in the plane of interest. The electric field magnitudes to which the patterns 
are normalized are displayed below their respective plots. Radiation patterns are 
plotted for the phi=pi/2 and phi=0 vertical planes parallel and perpendicular to the 
y-axis respectively. A horizontal radiation pattern is plotted at an elevation 
selected by the index from the elevation angle index table. Polarization is a 
varying combination of horizontal and vertical depending on spatial orientation with 
respect to the antenna. 


NOTE: The Yagi-Uda mutual impedance calculations are not valid if there is an 
antenna element whose length is an exact integer multiple of the wavelength. If 

(his occurs, there will be a singularity error in the mutual impedance calculations. 
This should not be a problem since Yagi-Uda Antennas are designed to operate at 
frequencies such that the driven element is approximately one-half-wavelength long. 
the reflectors are slightly longer, and the directors are slightly shorter. If the problem 
does arise, it will be necessary to vary the operating frequency such that no element 
is exacily an integer multiple of the wavelength. 


Input the number of NR. į Input the numberof — wj . 
reflector elements director elements 
y - 0. (NR + ND) $ = O..(NR + ND - 1) 


This entire page is alloted for these entries to allow for as many elements as necessary. 


Input the lengths and radii of the Yagi 
elements, starting with the outermost 
reflector and proceeding inward to the 
driven element. After entering the driven 
element length and radius, enter those of 
the director elements, starting with the 
director next to the driven element and 
proceeding outward to the final element. 


l- rad = 
Y 





Zo 


Input the distances between 
successiv elements in the same 
order as the length and radius 
inputs. There should be one less 
entry in the separation distance 
array than there is in the length 
and radius array. 


Input the operating fe - 10- 10° Input the heightof Hg = 7.5 
Frequency (Hertz) the array above the 
surface (meters) 


Input the ground o - 4 Input the ground m o -7 
Conductivity Dielectric Constant 

Input the Distance p _ 3000 Index of Refraction 
from the Antenna 

(meters) 





Input the index of the elevation 

angle for which to calculate the W = 285 
horizontal radiation pattern 

(from the angle index table) 


Complex Numerical Distance for Vertical Polarization 


-jB (R + Hg cos(0;)) — n? - sin (8) f 
2 


2- sin (0;) n 


Pe; = 


Complex Numerical Distance for Horizontal Polarization 


Po -j-B-(R + Le (Ras Ho Salt) s fe? - (since?) | 


2-( sin i 2) 


(n?-cos(&) ) - n^ - in (4)? 


0%cos(0y)) + (fas : ! 


2 : 2) 
Horizontal Reflection Coefficient Fh; = cos(6) J a sin (6j) 


cos(6)) + (h? . -— 


Vertical Reflection Coefficient Py; 'z 


Radiation Patterns in Phi=Py2 Plane 
(Perpendicular to X-Axis and Dipole Elements) 
Space Wave Radiation Pattern (Phi=Pi/2) 





-1 -0.5 0 0.5 1 
nek, HIP, 


Max E-Field Intensity (Volts per meter) max( MagE!P) = 6.85595: [UR 
Surface Wave Radiation Pattern (Phi-pi/2) 





net, l Hep, 


Max E-Field Intensity (Volts per meter) max( MagE2P) = 2.16533- 19 19 
209] 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 





-1 
Het, : «JP. 


Max E-Field Intensity (Volts per meter) max( MagE3P) = 6.85595: 10° 


27 
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Radiation Patterns in Phi=0 Plane 
(Perpendicular to Phi=pi/2 Plane) 


Space Wave Radiation Pattern (Phi=0) 





-1 -0.5 0 0.5 1 
Heeb. #10, 


Max E-Field Intensity (Volts per meter) max( MagE10) - 2.32256: 10 * 
Surface Wave Radiation Pattern (Phi=0) 
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Max E-Field Intensity (Volts per meter) ^ max( MagE20) - 9.86965: io? 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 
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Heb : H30; 


Max E-Field Intensity (Volts per meter) — max( MagE30) - 2.32256: io 
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Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 


Space Wave Radiation Pattern (Horizontal) 
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nre£, HIB, 


Max E-Field Intensity (Volts per meter) max( MagE1H) - 425018- 10 * 


Surface Wave Radiation Pattern (Horizontal) 





Heb, 1 «2H. 


Max E-Field Intensity (Volts per meter) max(MagE2H) = 3.10467- 10° 


21S 


Combined Space and Surface Wave Radiation Pattern (Horizontal) 
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Max E-Field Intensity (Volts per meter) max( MagE3H) - 4.25018: 10 * 
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* 1.5707963 
power0 = Ta 
TR 
| --1.5707963 -'0 
LE ls 
2 2 
30 
powerg - 5 
x.R 
5 ta 
2 


power = power + powerø 


— (| MagE3P;|) f 
a. + 0-3) 

— (| MagE30, |)? 
A OO) 


Total Power Radiated (Watts) 


Radiation Resistance (versus 
Maximum Input Current) 
(Ohms) 


1 - (sin(¢)-cos(t))? 


1.5707963 
Ñ E Dale eD gai: 
na $$" 2 
y | 1 - (sin( $)-cos(€) )? 
SUR SU us MA ents | 


p 
\ 


e (eos(8: 


( VREFyg) e 


2: power 


Radres - 


4. x- R2: max( sqE3P) 
power 


DirectivityP - 


AMARE max ( sqE30) 
power 


DirectivityO -= 


power = 0.00482 


Radres = 103.76004 


Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 


Phi = pi/2 Plane 


DirectivityP = 14.63143 


Phi = 0 Plane 


DirectivityO = 1.67913 


Cy 


Effective Isotropic Radiated Power (EIRP) (Watts) 
Phi = py/2 Plane Phi = 0 Plane 
DirectivityP: power — 0.07051 DirectivityO- power = 0.00809 


Maximum Effective Area (Along Radial of Directivity) 
(square meters) 





Phi = pi/2 Plane Phi = 0 Plane 
(As) : Directivit yP 4 (As) : DirectivityO 
~—_______— = 1.0479. 10 >~ = 120.2591! 
4. x 4. 
Maximum Effective Length (Along Radial of Directivity) 
(meters) 
Phi = pi/2 Plane Phi = 0 Plane 
Zu etu À Tieentivi 
Radres: ( Ae ) ^: DirectivityP Radres: (Ag) ^ DirectivityO 
2° — = 33.9559 2. s ii = {1.50637 
480- 2? 480. x? 
Numerical Distances 
Vertical Polarization Horizontal Polarization 
| Peg| - 0.04363 | Pmg| = 2.26206- 10° 


Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 


Phi = pi/2 Plane Phi = 0 Plane 
89.71429 
AngleP = 41.71429 AngleO = 
89.71429 


Oe 
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